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ay I. Allgemeine Ubersicht Sage 


Die ferromagnetischen Stoffe verdanken ihre physikalisch wichtigste Eigen- | 

- schaft, die spontane Magnetisierung, der Tatsache, dass zwischen benachbarten : 
¥ Elementarmagneten ihres Kristallgitters Austauschkrafte wirken, welche diese — 
- Dipole parallel zu stellen suchen. Infolge dieser Krafte bildet sich unterhalb — 
einer gewissen kritischen Temperatur, der Curie-Temperatur, ein durch eine 
_magnetische Polarisation charakterisierter Ordnungszustand aus. 
In der Natur treten jedoch ausser den die Parallelstellung benachbarter 
Dipole begiinstigenden Austauschkraften mindestens ebensohaufig solche auf, 
_ welche die Antiparallelstellung benachbarter Dipole bevorzugen. Genau wie im _ 
Falle des Ferromagnetismus bildet sich daher unterhalb einer bestimmten kriti- _ 


-schen Temperatur (aus Analogiegriinden ebenfalls Curie-Temperatur genannt) A 


-spontan ein Ordnungszustand in bezug auf die Richtungsverteilung derDipoleaus. _ | 
Wenn dieser Ordnungszustand durch eine resultierende Magnetisierung . 
_gekennzeichnet ist, spricht man von Ferrimagnetismus, bei verschwindender — 
7 4 . . . . a 2 = 
- resultierender Magnetisierung dagegen von Antiferromagnetismus’). 
1) AFIF, Eidgenéssische Technische Hochschule. ' 
2) Diese Bezeichnungen stammen von N&éex, der diese Erscheinungskomplese als erster 
pearbeitet hat. 
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Ferrimagnetismus tritt haufig bei komplizierteren Verbindungen auf, bei 
welchen die wechselwirkenden Dipole auf verschiedenartigen Gitterplatzen 
sitzen. Die Eigenschaften dieser Substanzen sind ahnlich denen der Ferro- 
magnetika. 

Antiferromagnetismus tritt dagegen vor allem bei einfachen anorganischen 
Verbindungen, eventuell auch bei Elementen auf. Die magnetische Suszeptibi- 
litat dieser Stoffe liegt mit derjenigen gewdhnlicher Paramagnetika in der 
gleichen Grdéssenordnung, zeigt jedoch unterhalb des Curie-Punktes einen 
anderen Temperaturverlauf sowie ausgesprochene Feldabhangigkeit. Am Curie- 
Punkt findet man eine durch eine Gitterverzerrung bedingte Anomalie der 
thermischen Ausdehnung wie auch die beim Ubergang von Order zu Disorder 
zu erwartende Anomalie der spezifischen Warme. 

Um einen Uberblick iiber die bis heute gewonnenen Erkenntnisse auf dem 
Gebiete des Antiferromagnetismus zu vermitteln, sollen im folgenden die Fra- 
gen besprochen werden, die mit der Wechselwirkung zwischen benachbarten 
Dipolen, mit der Art des Ordnungszustandes, mit den magnetischen Eigen- 
schaften und mit den Erscheinungen beim Ubergang von Order zu Disorder 
zusammenhangen. 


II. Die Wechselwirkung der Elementarmagnete 


Im Falle der ferromagnetischen Stoffe steht fest, dass fiir die Wechselwir- 
kung der Elementarmagnete quantenmechanische Austauschkrafte verantwort- 
lich gemacht werden miissen. Alle andern Méglichkeiten, zum Beispiel die 
magnetostatische Dipol-Dipol-Wechselwirkung, liefern viel zu kleine Energien. 
Man muss daher im Falle elementarer Antiferromagnetika, deren Curie-Punkte 
in vielen Fallen wesentlich iiber Zimmertemperatur liegen, einen dhnlichen 
Wechselwirkungsmechanismus erwarten. 

Die meisten der bekannten Antiferromagnetika sind jedoch Verbindungen, 
in deren Kristallgitter die Dipole durch unmagnetische Ionen separiert sind. 
Diese experimentelle Tatsache beweist, dass auch indirekte Austauschkrafte 
bestehen miissen (Superexchange). 


1. Direkte Wechselwirkung 


Die Energie zweier benachbarter Elementarmagnete betragt nach HEISEN- 
BERG [27]}) 


W=-2A(S, 5), (1) 


wo S,und S, die Spinvektoren der beiden Dipole in Einheiten //(2z) und A das 


AAMAS a : 3 3 : 
) Die Ziffern in eckigen Klammern, beziehen sich auf das Literaturverzeichnis auf Seite 22. 


— 
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Austauschintegral 

A =| Yat, Va %1) Yo(%e, Vo, 22) HZ Pal%s, Vo, 2) Wola, Vip %) dv, dv, (2) 
bedeuten. y, und y, sind die Eigenfunktionen, die den beiden beteiligten Elek- 


tronen zur Verfiigung stehen, deren Koordinaten die Indizes 1 bzw. 2 tragen. 
.H bedeutet den Hamilton-Operator. 


6 7 
Atomeatstang 
Ochalenradius 


Fig. 1 
Abhangigkeit des Austauschintegrals vom Verhaltnis Atomabstand : Schalenradius. 


Aus Gleichung (1) geht hervor, dass fiir positives Austauschintegral Parallel- 
stellung der magnetischen Momente, fiir negatives Austauschintegral dagegen 
Antiparallelstellung die energetisch giinstigere Lage bildet. 

Das Vorzeichen von A wurde von SLATER [59] diskutiert. In Figur 1 ist die 
za erwartende Abhangigkeit des Austauschintegrals vom Verhaltnis zwischen 
Atomabstand und Radius der d-Schale dargestellt. Antiferromagnetismus ware 
demnach fiir Elemente mit grosser d-Schale oder kleinem Atomabstand zu er- 
warten. Darunter fallen Mn, Cr, V, Ti, Sc in der ersten Ubergangsreihe, Ru, 
Rh, Pd in der zweiten und Os, Ir und Pt in der dritten Ubergangsreihe. NEEL 
[49] berechnete die Grésse des Austauschintegrals aus dem paramagnetischen 
Verhalten dieser Metalle sowie ihrer binéren Legierungen. Die Slatersche 
Kurve fand sich weitgehend bestatigt. Pd, Mn, Cr, Rh, Pt, V, Ru, Ti, Ir, Os 
zeigen in dieser Reihenfolge zunehmend negative Austauschwechselwirkung. 

ZENER [84] bis [87] vertritt die Auffassung, dass das Austauschintegral 
zwischen den d-Elektronen immer negativ sei und dass die starke, zum Ferro- 
magnetismus fiihrende positive Wechselwirkung zwischen den Dipolen durch 
Vermittlung der Leitungselektronen zustande komme. Auf dieser Basis gelingt 
die Erklarung einer Reihe von Erfahrungstatsachen, zum Beispiel der raum- 
zentrierten Kristallstruktur der Ubergangsmetalle. 
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2. Indirekte Wechselwirkung 


Das Problem der indirekten Wechselwirkung fand seine erste quanten- 
theoretische Bearbeitung durch KRAMERs [35]. ANDERSON [1] legte einer spate- 
ren Behandlung des Effektes das folgende Modell zugrunde: Man betrachtet 
von zwei Metallionen Met++ mit magnetischem Moment je ein Elektron im 
Zustand d, bzw. do. 


Zwischen diesen Metallionen sitzt ein unmagnetisches Kopplungsion, zum 


Beispiel O-~, von welchem zwei p-Elektronen in den Austauschmechanismus 
einbezogen werden. Schematisch : 
| vg ly ved? 
Met+ O-- Mett 


Aus diesem Grundzustand lasst sich ein angeregter Zustand ableiten durch 
Ubergang eines #-Elektrons des O-~ zu einem der Me**. 


da, > ds 
' Met O- Mett 


Zu diesem Vorgang werde die Energie. Ey bzw. Ey bendétigt, je nachdem die 
Spin der d,- und dj-Elektronen in Me* parallel oder antiparallel stehen. 

Die quantenmechanische Berechnung zeigt bei st6rungstheoretischer Be- 
riicksichtigung des angeregten Zustandes, dass die Energie des Grundzustandes 
tatsachlich von der gegenseitigen Orientierung der Spin der d,- und d,-Elek- 
tronen abhangt, sofern die #- und d,- bzw. f- und d,-Eigenfunktionen sich ge- 
niigend tiberlappen und Ey, verschieden von Ey ist. 

Die Winkelabhangigkeit der Wechselwirkung ist wie im Falle der direkten 
Wechselwirkung durch ein Kosinusgesetz gegeben : 


| W = —2 Aing (Si Se) (3) 
mit 
1 1 
A . => “j;5 — Ta 2 , 
ind (2 zr) Q y (4) 
worln 
0 = | wa(x, ¥, 2) H pp(x, y, 2) dv (5) 


mit dem Hamilton-Operator H das 1-Elektron-Resonanzintegral bedeutet und 
fiir J das 2-Elektronen-Austauschintegral 


a aa Ya, (Xr Vas 21) Wo(%e, Vo, 22) H Ya, (%2, Var 22) Wo(%1, Va» 21) dv, dv, (6) 


Zu setzen ist. 


— 
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Der Ausdruck (4) lasst einige physikalisch direkt kontrollierbare Schliisse zu. 


a) Das Vorzeichen von A ing 


Da / einer normalen chemischen Bindung entspricht, bei welcher die Spin 
sich absattigen, ist es negativ. Damit A,,q negativ wird, wodurch Antiferro- 
magnetismus bedingt ist, muss daher | Ey|<|E,| gelten, das heisst, die 
antiparallele Anlagerung eines weiteren Elektrons im Me++ muss energetisch 
giinstiger sein. Wie aus magnetischen und spektroskopischen Untersuchungen 
bekannt ist, trifft dies bei den zweiwertigen Ionen der ersten Ubergangsreihe 
Sce++, Tit+, Vt+, Cr++, Mnt+, Fet+, Co++, Ni++ nur fiir die Tonen zwischen 
Mn++ und Nit+ zu, weshalb zum Beispiel MnTe antiferromagnetisch, CrTe 
dagegen ferromagnetisch gefunden wird. 


b) Die absolute Grisse von A jnq 


Sie hangt besonders stark von der Grosse @ ab. Es ist daher verstandlich, 
dass stark elektronegative Kopplungsionen, bei denen die Uberlappung der 
w-Funktion gering und damit o@ klein ist, zu kleinerer Wechselwirkung Anlass 
geben als weniger stark elektronegative. Man erwartet daher, dass zum Beispiel 
in der Reihe MnO, MnS, MnSe, in welcher die Elektronegativitat des Kopplungs- 
ions abnimmt, die Wechselwirkung zunimmt, was durch die Reihe der entspre- 
chenden Curie-Punkte (122° K fiir MnO, 165° K fiir MnS und 247° K fiir MnSe) 
experimentell bestatigt wird. 


c) Winkelabhangigkeit 


Die den #-Elektronen im Kopplungsion entsprechende Ladungsverteilung 
ist hantelformig, was bedingt, dass 9 am gréssten und damit die Wechselwirkung 
am starksten wird, wenn die drei betrachteten Ionen in einer Geraden liegen. 


III. Curie-Temperatur und Ordnungszustand 


Die Frage nach dem sich unterhalb der Curie-Temperatur infolge der Wech- 
selwirkung zwischen den Dipolen einstellenden Ordnungszustand lasst sich bei 
den Antiferromagneten im wesentlichen nach denselben Methoden behandeln 
wie bei den Ferromagneten. 

Im Vordergrund steht die bei grosser Einfachheit durch thre erstaunliche 
Leistungsfahigkeit sich auszeichnende Molekularfeldtheorie, welche sich aus 
der alten Langevin-Weiss-Theorie des Ferromagnetismus ableitet. Die Schwache 
dieser Methode besteht vor allem in der Vernachlassigung der Riickwirkung 
der Dipole auf ihre Umgebung. 

Das Problem ist deshalb auch mit verschiedenen exakteren Methoden an- 
gegriffen worden, die sich aber haufig bei tragbarem mathematischem Aufwand 
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nur auf die lineare Kette anwenden lassen. Eine Aufzaéhlung dieser Arbeiten 
folgt Im zweiten Abschnitt dieses Kapitels. 

Fiir die theoretisch abgeleiteten Ordnungszustainde besteht eine experi- 
mentelle Priifungsméglichkeit durch Interferenz von Neutronen am Gitter, 
iiber welche wir im dritten Abschnitt berichten. 


1. Molekularfeldtheorien 


Da nach Definition des Antiferromagnetismus der unterhalb des Curie- 
Punktes sich einstellende Ordnungszustand ohne dusseres Feld keine resultie- 
rende Magnetisierung besitzen soll, miissen wir annehmen, dass sich die Elemen- 
tarmagnete innerhalb sehr kleiner Bereiche absattigen. Aus Energiegriinden 
werden hiebei zum Beispiel nach links gerichtete Dipole im Gitter periodisch 
mit nach rechts gerichteten Dipolen abwechseln. Diese Uberlegung fiihrt zwangs- 
laufig dazu, ein Antiferromagnetikum als aus zwei oder mehreren ineinander- 
geschachtelten Untergittern bestehend zu behandeln, deren jedes ein durch 
die Magnetisierung des andern beeinflusstes Ferromagnetikum oder Antiferro- 
magnetikum darstellt. 

Diese auf N&EL [49], [50], [51] zuriickgehende Behandlungsweise wurde 
spater durch ANDERSON [2] und VAN VLECK [75], [76] weiter ausgebaut. In 
bezug auf die magnetischen Ionen einfach kubische und kubisch raumzentrierte 
Gitter lassen sich, wie die Figuren 2 und 3 zeigen, ohne weiteres so in zwei Unter- 
gitter A und B aufteilen, dass jeder Elementarmagnet von A lauter solche von 
B als nachste Nachbarn hat und umgekehrt. Bei kubisch flachenzentrierten 
Gittern sind, wie man sich an Hand von Figur 4 iiberzeugt, die nachsten Nach- 
barn eines herausgegriffenen Dipols zum Teil unter sich wieder nachste Nach- 
barn. Der Idealfall, in dem alle nachsten Nachbarn eines beliebigen Dipols zu 


Fig. 2 


Aufteilung des einfach kubischen Gitters in zwei antiparallele Untergitter. 
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Fig. 3 
Ordnung erster Art im kubisch raumzentrierten Gitter. 


Fig. 3a 


Fig. 4 
Aufteilung des kubisch flachenzentrierten Gitters in vier Untergitter A, B, C und D. 
Ordnung erster Art. 
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Fig. 4a 
Ordnung zweiter Art im kubisch flachenzentrierten Gitter. 


diesem antiparallel orientiert sind, lasst sich daher beim flachenzentrierten 
Gitter nicht realisieren. Zur Behandlung dieses Gittertyps ist eine Aufteilung 
in vier einfach kubische Gitter A, B, C und D notwendig. 


a) Einfach kubische und kubisch raumzentrierte Gutter 


Das innere Feld H; , am Ort eines Dipols A wird einerseits der Magnetisie- 
rung M, des Untergitters A und anderseits derjenigen des Untergitters B pro- 
portional gesetzt : 

H,,=—-«M,—yMy, (7a) 
und entsprechend 
H,z =—«M,z—yM,. (7b) 


Da gemiss der Unterteilung in Untergitter keine zwei Dipole desselben Unter- 
gitters nachste Nachbarn sind, entsprechen die Glieder —« My und —a« Mz, 
der Mitberiicksichtigung der zweitnachsten Nachbarn. 

Die Wechselwirkungskonstanten y und « hangen folgendermassen von den 
Zahlen Z4, und Z44 der nachsten bzw. iibernachsten Nachbarn sowie von den 
entsprechenden Austauschintegralen A 4, und A 44 ab: 


\ oo 2 2p Aang —2244 A444 
on Negipe ’ ~~ Nags - (8) 


N, und N B bedeuten die Zahlen der Atome der Untergitter A und B in der 
Volumeneinheit, g den Landé-Faktor und f das Bohrsche Magneton. Die Vor- 
zeichen in (7) sind so gewahlt, dass bei antiferromagnetischer Wechselwirkung 


a% und y positiv werden. Die Ausdriicke (8) kénnen sich je nach der Berech- 
nungsweise um kleine Faktoren dndern. 
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Fir hohe Temperaturen kann man Sattigungserscheinungen vernachliassi- 
gen und bei einem aussern Feld H fiir die Magnetisierung der beiden Untergitter 
entsprechend der Langevin-Theorie des Paramagnetismus das Gleichungs- 
system 


C = = au a Cc 25 ae a 
M,=5> (GoM p= M,), M, = So (H-aM,—yM,) (6) 
anschreiben, worin 


N= B29? S13 tod) 
3k = 


Ga 


(N = Zahl der Dipole je Kubikzentimeter, S = Spinquantenzahl der Elemen- 
tarmagnete, = Boltzmann-Konstante) und die Klammerausdriicke das in den 
Gitterpunkten A bzw. B wirkende totale Feld bedeuten. Beriicksichtigt man, 
dass die totale Magnetisierung M = M, + Mz, so erhalt man hieraus sofort 
fiir hohe Temperaturen eine Suszeptibilitat 


C 
Land Gag (10) 
mit einer sogenannten asymptotischen Curie-Temperatur 
i! 
O=-—Z(%+y)C. (10a) 


Lasst man andererseits das dussere Feld H verschwinden, so berechnet man aus 
der Forderung, dass das nunmehr homogene Gleichungssystem (9) eine nicht- 
triviale Lésung besitze (Determinante = 0), die eigentliche Curie-Temperatur 
T,, unterhalb der sich ein Ordnungszustand im Kristall aufbaut, zu 


T.=5 y-a)C. (11) 


Es ist nun interessant, zu bemerken, dass dem Verhaltnis der beiden 


Curie-Temperaturen 
OPS yt a 


4 a 


gewisse Grenzen gesetzt sind, indem sich zeigen lasst, dass fiir « > y/2 beim 
raumzentrierten Gitter eine Ordnung zweiter Art (Figur 3a) energetisch giin- 
stiger wird, bei welcher die Untergitter A und B in sich selbst antiferromagne- 
tisch werden. In Funktion von «/y wird deshalb @/T, den in Figur 5a gezeigten 
Verlauf nehmen, das heisst, wenn @ bekannt ist, kann man fiir 7, bestimmte 


Grenzen angeben. 
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Fig. 5 
(Nach VAN VLECK [76].) 


b) Fldchenzentrierte Gitter 


Wenn M,, Mz, Mc und Mp die Magnetisierungen der vier einzufiithrenden 
Untergitter bedeuten, so wird fiir das innere Feld H; 4 am Ort eines Dipols A 
in analoger Weise zu setzen sein: 


Hy = —~aM, —y(M, + Me + Mp) | 
und = S: oy = a (12) 
H,» = —« M, —y(M,+M,+M,) usw. | 


Die entsprechende Durchfiihrung der Rechnung ergibt fiir eine Ordnung gemiass 


Figur 4, welche ANDERSON [2] als die bei kleinen positiven Werten von « ener- 
getisch giinstigste erkannte 


Ti aid aa S t- 4) pO gy Se 
(3 ay 4 \e, I, (5 12 C; 7p por y — a/3 . (13) 


Diese Ordnung geht bei « > 3 y/4 wieder in eine solche zweiter Art iiber, bei 
der die Untergitter gemass Figur 4a in sich antiferromagnetisch werden, so 
dass O/T, den in Figur 5) aufgetragenen Verlauf nimmt. 
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Eine auf derselben Grundlage beruhende, etwas allgemeiner gefasste Theorie 
wurde jiingst von SMART [60] verdffentlicht. Die verschiedenen méglichen 
Arten von Ordnungen kénnen auf Grund einer Arbeit von LUTTINGER [44] 
auch streng rechnerisch bestimmt werden. GorTER und HAANT]Es [23] modi- 
fizieren die Néelsche Theorie zur Anwendung auf rhombische Kristalle am 
absoluten Nullpunkt. 

Li [41] betont, dass die Lage des Curie-Punktes von der Grisse des ange- 
legten Magnetfeldes abhangt. In diesem Zusammenhang zeigt GARRETT [20], 
dass Antiferromagnetismus auch beim absoluten Nullpunkt nur méglich ist, 
solange das Magnetfeld kleiner ist als 


th (14) 


(u = magnetisches Moment der Ionen). Dieser Effekt wird besonders deutlich 
bei tiefen Umwandlungstemperaturen, weil es dann mit den laboratoriums- 
massig herstellbaren Feldern méglich ist, H, zu erreichen. Experimentell ist die 
Erscheinung von Squire [63] an Mn-Verbindungen, vor allem MnSe, von 
SHALYT [55] an FeCl,, CoCl, und NiCl,, von PouLis, HANDEL, UBBINK, POULIS 
und GorRTER [53] an CuCl,, 2 H,O und von GorTER [22] und GARRETT [21] an 
Alaunen festgestellt worden. 

So leistungsfahig die Molekularfeldmethode auch scheint, so haften ihr, 
wie KRAMERS [36] zu bedenken gibt, doch wesentliche Schwachen an. Zum 
Beispiel wird nach den strengeren Berechnungen von BETHE [4] und HULTHEN 
[31] die lineare Kette nie antiferromagnetisch, wahrend die Molekularfeld- 
methode auch in diesem Falle einen Curie-Punkt liefert. 


2. Strengere Methoden 


Die Beriicksichtigung der Riickwirkung eines Dipols auf seine Umgebung 
gibt zu einer starken Komplizierung der mathematischen Behandlungsweise 
Anlass, so dass es in vielen Fallen nur méglich ist, den eindimensionalen Fall 
(lineare Kette) zu meistern. 

Nachdem BETHE [4] als erster eine strenge quantenmechanische Berechnung 
der Eigenfunktionen der linearen Atomkette versucht hatte, modifizierte HuL- 
THEN [31] diese Methode und wandte sie speziell auf den Fall der antiferroma- 
gnetischen Kette an. Er fand, dass dieses Modell keinen Curie-Punkt besitzt 
und somit auch bei tiefen Temperaturen in keinen ausgedehnten Ordnungszu- 
stand iibergeht. Spater verallgemeinerte KasTELEIJN [33] HuLttTuEns Rech- 
nungen unter Einbeziehung eines Kristall-Anisotropiegliedes. Er konnte zeigen, 
dass dieses einen Minimalwert besitzen muss, wenn sich ein Ordnungszustand 
einstellen soll. 

WANNIER [79] bewies, dass ein ebenes Dreiecknetz bei noch so tiefen Tem- 
peraturen nicht antiferromagnetisch wird. 
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Hingegen erhalten WAKEFIELD [77] mit statistischen Methoden und AN- 
DERSON [3] auf Grund einer halbklassischen Modifikation der Blochschen 
Spin-Wave-Methode [9] durch HELLER und KRAMERS [28] im raumlichen Modell 
schon fiir das einfach kubische Gitter bei geniigend tiefen Temperaturen einen 
Ordnungszustand, dessen Energie innerhalb gewisser Grenzen abgeschatzt 
werden kann. 

Die auf BETHE [5] zuriickgehende statistische Methode zur Behandlung des 
Ferromagnetismus ist neuerdings von Lr [42], in der durch WEIss [80] ausge- 
arbeiteten Form, auf Antiferromagnetika mit einfach kubischem und raum- 
zentriert kubischem Gitter angewendet worden. Obschon die Wechselwirkung 
mit den iibernichsten Nachbarn vernachlassigt wurde, was bei der Molekular- 
feldtheorie auf —@ = T, fihrt, erhalt L1 -O=1,5 T, bzw. —O = 1,25 T,. 
Diese Faktoren reichen aber nicht aus, um die grossen beobachteten Quotienten 
Q/T, bis 5 zu erklaren, so dass vorlaufig die Betrachtung der tibernachsten 
Nachbarn nicht entbehrlich scheint. 

ZIMAN [88] leitet ebenfalls die Hilfe des Bethe-Verfahrens im wesentlichen 
die gleichen Resultate ab wie GARRETT [20] mit der Molekularfeldtheorie. Eine 
weitere Bearbeitung mit den statistischen Methoden der Theorie der Uber- 
struktur in Metallen hat Sato [54] publiziert. 

Brooks und Doms [10] berechnen das zweidimensionale antiferromagne- 
tische Ising-Modell auf Grund einer von KRAMERS und WANNIER [37] vorbe- 
reiteten strengen Anwendung der Boltzmann-Statistik. 

TsuyA [71] gelang die Darstellung der Temperaturabhangigkeit der Suszep- 
tibilitat von MnSe durch Anwendung der Energiebandtheorie des Ferromagne- 
tismus von MryAHARA [47] auf den Fall negativer Austauschintegrale. 

STREET [65] soll die theoretischen Arbeiten bis 1951 resiimieren. 


3. Neutronenbeugungsexperimente 


Die experimentelle Bestatigung der schon 1936 von NEEL [49] vorausge- 
sagten Aufteilung der Dipole eines Antiferromagnetikums in zwei (oder mehrere) 
antiparallel magnetisierte Untergitter liess mangels passender Methoden lange 
auf sich warten. Die Réntgenmethoden mussten fiir einen direkten Nachweis 
des Ordnungszustandes versagen, weil der Réntgen-Streuquerschnitt von der 
Orientierung der magnetischen Dipole gegen die Einfallsrichtung der Réntgen- 
quanten unabhangig ist. : ; 

1949 fanden SHULL und Smarr [56] in der Neutronenstreuung eine Methode 
zum. direkten Nachweis magnetischer Ordnungszustande. Fiir kleine Streu- 
winkel liegt der Neutronen-Streuquerschnitt der Hiille in der gleichen Gréssen- 
ordnung wie derjenige des Atomkerns und hat die Eigenschaft, von der Orien- 
tierung des Hiillenmomentes gegen die Streuebene abhangig zu sein. 

Fir einen Ordnungszustand gemiss Figur 3a zum Beispiel wird die magne- 
tische Elementarzelle linear doppelt so gross wie die chemische, was sich im 
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Neutronenbeugungsdiagramm durch Auftreten von Uberstrukturlinien aussern 
muss. 

Dieser Effekt ist tatsichlich bei den untersuchten Antiferromagneten ge- 
funden worden. Aus den Intensitatsverhiltnissen lasst sich schliessen, dass in 
MnO, MnS, MnSe, FeO, NiO und CoO, welche alle im NaCl- -Typ kristallisieren 
und daher in bezug auf die Metallionen flachenzentriert sind, unterhalb der 
Curie-Temperatur eine Ordnung zweiter Art (vgl. Figur 4a) sich einstellt. SHULL, 
STRAUSER und WOLLAN [57] konnten ebenfalls zeigen, dass die Dipole wie bei 
den Ferromagneten im Kristall an gewisse Vorzugslagen gebunden sind, in 
diesen Fallen die 100- und aquivalenten Richtungen. 

Die Tatsache, dass eine Ordnung zweiter Art vorliegt, besagt gemass Ab- 
schnitt III/1, dass die Wechselwirkung zwischen tibernichsten Nachbarn be- 
sonders stark sein muss. Nach den in II/2 dargestellten theoretischen Uberle- 
gungen ist dies auch zu erwarten, da ja in den in Frage stehenden Stoffen die 
ubernachsten Nachbarn mit dem die indirekte Wechselwirkung vermittelnden 
Kopplungsion auf einer Geraden liegen, wahrend zwischen nachsten Nachbarn 
sowohl die direkte Wechselwirkung infolge des grossen Abstands der Jonen als 
auch die indirekte Wechselwirkung infolge des gewinkelten Kopplungsweges 
klein sind. 

Ebenfalls eine Ordnung zweiter Art fand Erixson [13], [14] bei dem in 
bezug auf die Mn-Ionen tetragonal raumzentrierten MnO,. Hingegen zeigen 
MnF, und FeF, bei gleichem Kristallgitter eine Ordnung erster Art. Die Vor- 
zugsrichtung der Spin steht senkrecht zur c-Achse. 

Bei dem in der NiAs-Struktur kristallisierenden CrSb fand SNow [62], dass 
alle Dipole in einer Netzebene senkrecht zur c-Achse unter sich parallel, jedoch 
antiparallel zu den Dipolen der benachbarten Netzebenen stehen, wodurch 
eine Schichtstruktur resultiert. Die magnetische Vorzugslage ist die c-Achse. 

Von den Ubergangsmetallen (vgl. II/1), welche nach NEEL [49] antiferro- 
magnetisch sein sollten, haben SHULL und WILKINSON [58] V, Cr, Mn, Cb, Mo 
und W nach der Methode der Neutroneninterferenzen untersucht. Sie konnten 
aber nur bei Cr mit Sicherheit einen antiferromagnetischen Ordnungszustand 
feststellen. Der Curie-Punkt scheint im Widerspruch zu den spater zu be- 
sprechenden Suszeptibilitatsmessungen von McGuire und KrigssMan [46] bei 
zirka 150°C zu liegen. Die Resultate beziiglich Mn sind noch nicht ver6ffentlicht. 


IV. Magnetische Eigenschaften 


1. Doménenstruktur, Bloch-Wdnde, H ystereseerscheinungen 


Man kann sich mit N&EL [50] fragen, ob der Ordnungszustand unterhalb 
der Curie-Temperatur sich gleichmassig tiber den ganzen Kristall ausdehnt 
oder ob in Analogie zu den Ferromagneten Weiss-Bezirke vorliegen. An sich 
stellen diese ja im ferromagnetischen Fall das Mittel dar, das die Natur benutzt, 
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um die magnetostatische Energie der spontanen Magnetisierung méglichst tief 
zu halten. Beim Antiferromagnetismus ware wegen der verschwindenden spon- 
tanen Magnetisierung aus diesem Grunde eine Aufteilung in Weiss-Bezirke 
nicht nétig. Wir wissen aber vom Ferromagnetismus her, dass infolge von in- 
neren mechanischen Spannungen einzelne, im spannungsfreien Zustand mit 
andern gleichberechtigte Magnetisierungsrichtungen energetisch bevorzugt 
werden. Andert der Spannungszustand in einem antiferromagnetischen Kristall 
rdumlich, so lasst sich vermuten, dass damit auch antiferromagnetische Weiss- 
Bezirke auftreten. Ihre Form und Grdsse sind aber im Gegensatz zum ferro- 
magnetischen Fall nicht wesentlich durch die Gestalt des Kérpers, sondern 
vielmehr durch seinen Spannungszustand, Reinheit usw. bestimmt. 

Wenn Weiss-Bezirke vorhanden sind, miissen auch Grenzen zwischen ver- 
schiedenen Bezirken auftreten. Man kann wiederum durch Vergleich zum 
Ferromagnetismus auf die Struktur dieser Bloch-Wande schliessen. Die 6rt- 
liche Anderung der Spinrichtung in diesen Ubergangsgebieten erfolgt so, dass 
die Energie, die sich aus Austauschenergie und Kristallenergie zusammensetzt, 
iiber die Wanddicke integriert, minimal wird. Bei den Antiferromagneten sind 
diese Beitrage von der gleichen Gréssenordnung wie bei den Ferromagneten. 
Man kann daher fiir die Dicke der Wande ebenfalls zirka 100 Atomabstande 
und fir ihre Energie ungefahr 1 erg/cm? einsetzen. 

Da die Bloch-Wande der Antiferromagnete gleiche Struktur wie diejenigen 
der Ferromagnete besitzen, wird ihre Beweglichkeit im Kristallgitter auch von 
der gleichen Gréssenordnung sein, das heisst, man braucht fiir eine Wandver- 
schiebung denselben «magnetischen Druck», welcher durch den Energiegewinn 
bei der Verschiebung eines Quadratzentimeters Wand um 1 cm dargestellt wird. 
Die kritische Feldstarke H;,;,, bei welcher im Antiferromagnetikum Wand- 
verschiebungen vor sich gehen, lasst sich daher aus der folgenden Gleichung 
abschatzen 


Js He = % Hiyi - é (15) 


Setzt man hierin fiir die spontane Magnetisierung J, des Ferromagnetikums 
1700 Gs, fiir seine Koerzitivkraft H, 1 Oe und fiir die Suszeptibilitat des Anti- 
ferromagnetikums ¥ = 3 - 10-4, so wird H;,,;, von der Gréssenordnung 3000 Oe, 
das heisst die Wandverschiebungen gehen erst bei wesentlich héheren Feld- 
staérken als in Ferromagneten vor sich. Entsprechend erwartet man die Dre- 
hung der Magnetisierungsrichtung, die sich bei Ferromagneten im allgemeinen 
oberhalb 100 Oe abspielt, erst bei Feldstiarken iiber 30000 Oe. 

Der Hysteresis der Magnetisierung bei den ferromagnetischen Stoffen ent- 
spricht in unserem Falle eine Hysteresis der Suszeptibilitat. 

Eine solche ist bis heute allerdings noch nicht beobachtet worden. Wohl 
aber stellten GORTER [22] an CrK(SO,).12H,O und Kurt [38] an ammoniaka- 
lischem Eisenalaun bei 0,03° K deutliche Hysterese der Magnetisierung fest. 
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2. Die Suszeptibilitat 


Die Berechnung der Suszeptibilitat oberhalb des Curie-Punktes ist in Ab- 
schnitt III/1 beschrieben. Im Anschluss an Arbeiten von N&EL [49], [50], 
BITTER [7] und VAN VLEcK [75], [76] lasst sich die Behandlung des Suszepti- 
bilitat und Magnetisierung der Untergitter unterhalb 7, folgendermassen skiz- 
zieren : 

Die Magnetisierung der Untergitter ergibt sich aus der Losung des folgenden 
Gleichungssystems: 


- ; pa aly H+ H, 
4 = Myo B; fee M; = Mz. B; ay eh (16) 


M4. und Mz, sind die Sattigungsmagnetisierungen der Untergitter A und B 
am absoluten Nullpunkt, « das magnetische Moment der Dipole und’B;,(z) die 
Brillouin-Funktion fiir den Spin 7 


Bz): = hs ce : ctgh ( es - : 2) — a ctgh et : (17) 
Fiir H; 4 und H;, miissen die Ausdriicke (7) bzw. (12) eingesetzt werden. 

Es ist klar, dass eine geschlossene Lésung von (16) fiir beliebige Tempera- 
turen nicht mdéglich ist. Die Berechnung der Suszeptibilitat 7 = (M. at M. B) /H 
ist daher unter der Annahme durchgefiihrt, dass nur die Richtung, nicht aber 
der Absolutbetrag der Magnetisierung der Untergitter durch das aussere Feld 
verandert werde. Gleichzeitig wird angenommen, dass sich die Kristallaniso- 
tropie durch die gegentiber der Kopplungsenergie der Untergitter kleine Energie 
W, = —K,cos2 # darstellen lasse, wo # den Winkel zwischen der magneti- 
schen Vorzugsrichtung A und der Magnetisierungsrichtung bedeutet. Das 
Magnetfeld H schliesse den Winkel B mit A ein. 

Unter der Wirkung des Feldes drehen sich die Magnetisierungen der Unter- 
gitter, wie aus Figur 6 ersichtlich, aus der Vorzugsrichtung J heraus und sind 
nicht mehr genau antiparallel, was zu einer resultierenden Magnetisierung M 


fiihrt. Da wegen der starken Wechselwirkung zwischen M 7, und M; 3 der Winkel 
zwischen diesen Magnetisierungen noch nahezu ein gestreckter ist, kann man 
nach wie vor von einer Richtung & der antiparallelen Magnetisierung sprechen. 
Ihr Winkel mit dem dusseren Feld H sei y. 

Aus der Bedingung fiir das Momentengleichgewicht erhalt man 


sin 2 B 
tg 2 Pests cos2 B — (H?/2 y K,) ; (18) 


Man berechnet hieraus unter der Annahme, dass sich im Mittel die Komponente 


16 HeINrIcH LABHART ZAMP 


der resultierenden Magnetisierung senkrecht zu H aufhebt, die Suszeptibilitat zu 


C= a sin? p. (19) 

Die Ausdriicke (18) und (19) enthalten die folgenden Aussagen : 
a) Fir Felder H < Ay, it = 2 y K, bleibt die Magnetisierung in der Vor- 
zugslage liegen. Oberhalb der kritischen Feldstarke H,,;, Tichtet sie sich da- 


A 


Fig. 6 Fig. 7 
Zur Berechnung der Theoretische Kurven der 
Suszeptibilitat unterhalb T,. Temperaturabhangigkeit der Suszeptibilitat. 


gegen in zunehmendem Mass senkrecht zum dusseren Feld. H,,,;; hat bei den mei- 
sten Antiferromagneten die Gréssenordnung 104 Oe (A, = 10° erg/em? y = 500). 

b) Die Suszeptibilitat ist bei schwachen Feldern abhangig vom Winkel 
zwischen H und A. Sie verschwindet fiir 8 = 0 und hat ihr Maximum ¥ nae = 1/y, 
wenn H | A. 

c) In der Nahe der kritischen Feldstarke wird sie feldabhangig, weil dann w 
eine Funktion des Feldes ist. Die Suszeptibilitat steigt mit zanehmendem 
Felde. 

d) Da bei H | A die temperaturabhangige Kristallenergie und Magneti- 
sierung der Untergitter keine Rolle spielen (die Magnetisierung bleibt bei allen 
Feldstarken in A liegen), ist unter den getroffenen Annahmen eine temperatur- 
unabhangige Suszeptibilitat y= 1/y zu erwarten. Dieser Wert schliesst, wie 
aus (10) und (11) ersichtlich, an den paramagnetischen Suszeptibilitatswert 
bei 7, ohne Sprung an. Der hieraus zu erwartende Temperaturverlauf der 
Suszeptibilitat ist in Figur 7 aufgetragen. Die Temperaturabhingigkeit von Xi 
foigt aus der Lésung des Gleichungssystems (16) fiir diese spezielle Orientierung. 
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Fiir einen Polykristall wird 


De iene 
x Lo Ki ; (20) 


Xoow = 3 


was erwarten lasst, dass die Suszeptibilitat am absoluten Nullpunkt auf zwei 
Drittel des Wertes beim Curie-Punkt T, gefallen ist. 

Diese theoretischen Resultate kénnen mit experimentellen Daten verglichen 
werden. Als Einkristalle sind vor allem FeCO,;, MnF,, MnO, und MnO durch 
BIzETTE [8] untersucht. Bei der ersten, rhomboedrisch kristallisierenden Sub- 
stanz sowie beim tetragonalen MnF, besteht eine Vorzugslage A fiir die Magne- 
tisierung. In Ubereinstimmung mit der Theorie ist die Suszeptibilitat unterhalb 
T, praktisch temperaturunabhangig, wenn H | A. Ebenso strebt sie erwar- 
tungsgemdss mit abnehmender Temperatur gegen Null, wenn H |] A. 

Beim ebenfalls tetragonalen MnO, bilden die a- und b-Achsen die Vorzugs- 
richtungen. Man findet daher bei Magnetisierung in c-Richtung temperatur- 
unabhangige Suszeptibilitat, wahrend in a- und b-Richtung immer eine Kom- 
_ ponente parallel H zeigt und deshalb gegen den absoluten Nullpunkt hin ein Ab- 
fall des Suszeptibilitat beobachtet wird. Das kubische MnO verhilt sich isotrop. 

Bei allen diesen Substanzen konnte unterhalb 7, die theoretisch verstind- 
liche Zunahme der Suszeptibilitaét mit steigender Feldstirke festgestellt werden. 

Figur 8 demonstriert an Hand von Messungen an polykristallinen Proben 
von MnF, und MnO, dass die theoretisch erwartete Temperaturabhangigkeit 
der Suszeptibilitat von Figur 7 gut bestatigt wird. Die Bedingung 7(0) = 2 ¥(T,)/3 
wird allerdings nur innerhalb sehr weiter Grenzen erfiillt. Hiertiber berichtet 
VAN VLECK [76] Naheres. 

Linpsay [43]' stellte bei MnSe eine sehr starke Temperaturhysterese der 
Suszeptibilitat fest, welche durch Annahme einer langen Einstellzeit fiir das 
thermodynamische Gleichgewicht zwischen verschiedenen Kistallmodifikatio- 
nen des MnSe erklart wird. 

Der charakteristische Temperaturverlauf der Suszeptibilitat von Figur 7 
und 8 kann als eines der hauptsdchlichsten Kennzeichen fiir das Auftreten von 
Antiferromagnetismus gewertet werden. In diesem Sinne interessant sind die 
erwaihnten Messungen von McGutrE und KriessMAn [46] an Cr, welche eine 
bis 1400° C stetig ansteigende Suszeptibilitat ergeben, was sich durch Antiferro- 
magnetismus mit sehr hohem 7, deuten liesse. Eigenartigerweise wird oberhalb 
150° C mit Neutroneninterferenzen keine Uberstruktur mehr gefunden. Ahn- 
liche Schliisse lassen die Suszeptibilitatsmessungen von HOARE und MATTHEWS 
[29] an Pt, Rh und Pd zu. Wahrend Pt bis zu 20° K hinunter einen negativen 
Temperaturkoeffizienten von x zeigt, ist derjenige von Rh im untersuchten 
Temperaturgebiet, das heisst hinauf bis zu 300° K, positiv, was wiederum auf 
Antiferromagnetismus mit relativ hohem 7, hindeutet. Infolge seines typischen 
Suszeptibilitatsmaximums bei 80° K kénnte Pd unterhalb dieser Temperatur 


ebenfalls antiferromagnetisch sein. 
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Fig. 8 
Experimentelle Kurven der Temperaturabhangigkeit der Suszeptibilitat von MnF, und MnO. 
Feldabhangigkeit von y (nach Bizerre [8)). 


3. Antiferromagnetische Resonanz 


300 T°K 


ZAMP 


Ein magnetischer Dipol besitzt in einem Magnetfeld H,;, eine Resonanz- 


frequenz 


(21) 


welche experimentell durch ein Absorptionsmaximum fiir elektromagnetische 
Wellen entsprechender Frequenz erkannt wird. Eine gewisse Schwierigkeit liegt 
in der Berechnung von H,,;, welches mit dem von aussen an einen Probek6rper | 
angelegten Feld durchaus nicht identisch ist, sondern wesentlich durch das ent- 
magnetisierende Feld Hp, durch das Austauschfeld H,; und durch das die Kri- 
stallanisotropie darstellende Feld H, mitbestimmt ist. Bei den Ferromagneten 
liegt H,,, immer in der Richtung der Magnetisierung M und spielt, wie aus der 


Bewegungsgleichung (D = Drall des Elementarmagneten) 


sty — [mM “ Hy] 


(22) 


ersichtlich ist, keine Rolle. Aus Abschnitt 1V/2, Figur 6, geht jedoch hervor, 


dass in Antiferromagneten die Richtung von M im allgemeinen nicht mit der 
Richtung von H,,, zusammenfallt, was bedeutet, dass unterhalb des Curie- 
Punktes H; nicht mehr vernachlassigt werden darf. Die Resonanzfrequenz w 
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verschiebt sich deshalb beim Abkiihlen unter den Curie-Punkt plétzlich stark. 
Da bei T, = 100° K H, in der Gréssenordnung 108 Oe liegt und im allgemeinen 
bei dusseren Feldern von einigen 103 Oe und Frequenzen um 102° Hz gearbeitet 
wird, verlagert sich die Resonanzstelle gewohnlich ausserhalb des Messbereiches 
der verwendeten Anordnungen. Hingegen konnten UsBink, Poutis, GERRIT- 
SEN und GorTER [74] im Falle der tiefen Curie-Temperatur von 4,3° K bei 
CuCl, 2H,O die Absorptionsfrequenzen infolge des entsprechenden kleinen H; 
weiter beobachten. 

Die oben skizzierte Theorie des Effektes wurde in grosser Allgemeinheit 
durch KEFFER und Kitret [34] und WanacsnEss [78] entwickelt. Eine friihere 
Bearbeitung stammt von Nacamrya [48]. UBBINK [73] behandelt den Fall 
eines rhombischen Kristalls im Hinblick auf die eben erwahnten Messungen an 
CuCl,, 2H,O. Messungen von OKAMURA, TorIzuKA und Koya [52] an MnO, 
MnS und MnSe, von Trounson, BLeIL, WANGSNESS und MAxweELt [70] an 
Cr,O3, von HUTCHISON [32] an MnF, und von BicKForp [6] an Fe,O, bei der 
Umwandlungstemperatur von 118° K zeigen alle teils schlagartiges, teils all- 
mahliches Verschwinden der Absorptionsmaxima unterhalb der Curie-Tempe- 
ratur. Bei Fe,O, ist das Auftreten von zwei neuen Absorptionsmaxima bei 
kleiner dusserer Feldstaérke beobachtet worden. 

Tsuya und IcHIKAWA [72] wie auch KEFFER und KITTEL [34] diskutieren 
die Breite der Absorptionslinien. 

MaxwELL [45] gibt einen kurzen Uberblick iiber die experimentellen Ar- 
beiten auf dem Gebiete der antiferromagnetischen Resonanz. 


V. Erscheinungen beim Ubergang vom paramagnetischen zum 
antiferromagnetischen Zustand 


Wenn der Ubergang vom paramagnetischen zum antiferromagnetischen 
Zustand wie beim Ferromagnetismus einem Ubergang von Disorder zu Order 
entspricht, so miissen beim Antiferromagnetismus die typischen, beim Ferro- 
magnetismus mit gerader Potenz von der Magnetisierung abhangigen Uber- 
gangserscheinungen auftreten, also Anomalie der spezifischen Warme, Gitter- 
verzerrung und Anomalie der thermischen Ausdehnung sowie Beeinflussung 
des Elastizitatsmoduls. 

Uber die Anomalie der spezifischen Warme liegen nur wenige Messungen vor. 
Zu erwahnen sind diejenigen von SHALyT [55] an NiCl,, CoCl,, FeCl, und CrCl,. 
Die erwarteten Maxima bei der Curie-Temperatur sind beobachtet, jedoch sind 
keine Angaben tiber ihre Héhe und ihre Flache publiziert. 

Gitterverzerrungen unterhalb des Curie-Punktes wurden in den Arbeiten von 
Tomps und Rooxssy [67], [68] sowie von GREENWALD und Smart [24], [25] 
nachgewiesen. Die oberhalb der Curie-Temperatur exakt kubischen Gitter von 
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MnS, MnO, FeO und NiO deformieren sich unterhalb 7, zu rhomboedrischen ~ 


Strukturen. Das kubische CoO wird tetragonal, wobei die c-Achse um 0,5% 
kiirzer als die a-Achse wird. Dieses unterschiedliche Verhalten ist erstaunlich, 
da nach den Shullschen Neutroneninterferenzversuchen der Ordnungszustand 
in CoO derselbe ist wie bei MnO, FeO und NiO. Cr,O3, das Korundstruktur 
besitzt, zieht sich, wenn es antiferromagnetisch wird, langs der 111-Richtung 
zusammen. 


» «10 
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Fig. 9 


Anomalie der thermischen Ausdehnung beim Curie-Punkt (nach Forx [16]). 


Solche Strukturaénderungen miissen sich auch im Verlauf des thermischen 
Ausdehnungskoeffizienten dussern, was Forx [16] (vgl. Figur 9) bei MnO, FeO 
CoO und NiO tatsichlich beobachtete. 

Eine weitere Folge der Strukturanderung ist der von STREET und LEwIs [66] 
bei NiO und CoO bei steigender Temperatur beobachtete Anstieg des Elastizi- 
tatsmoduls um 75 bzw. 150% am Curie-Punkt. Interessant ist auch die Beob- 
achtung einer wesentlich vergrésserten mechanischen inneren Dampfung im 
antiferromagnetischen Zustand, welche in ahnlicher Weise wie die grosse innere 
Dampfung der Ferromagnetika durch hysteresebehaftete Verainderung der 
Domanenstruktur zu erklairen versucht wird. 

Obschon nur lose in diesen Zusammenhang gehérend, sei noch auf die Beob- 
achtung von Forx und WucuHer [19] hingewiesen, wonach Cr,O, im amorphen 
Zustand erwartungsgemass nicht antiferromagnetisch wird. Hervorzuheben ist 
die gefundene Abnahme des magnetischen Momentes der Elementarmagnete 
im amorphen Zustand, welche méglicherweise durch eine Kompensation der 
Spinmomente durch die frei werdenden Bahnmomente verursacht wird. 

Endlich fiihren wir noch die Arbeit von CoRLISS, DELABARRE und ELLIot 
[12] an, welche sich auf Mischkristalle von MnF, mit ZnF, bezieht, bei welchen 
sich der Curie-Punkt linear mit dem ZnF,-Gehalt gegen Null verschiebt, was 
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einen Anhaltspunkt dafiir liefert, dass man in der Molekularfeldtheorie die 
Wechselwirkungskonstanten gemiass Gleichung (8) proportional der mittleren 
Zahl der nachsten magnetischen Nachbarn setzen darf. 


VI. Zusammenstellung der bis heute bekannten Antiferromagnetika 


In der folgenden Tabelle sind diejenigen Substanzen zusammengestellt, 
welche heute mit einiger Sicherheit als Antiferromagnetika gelten diirfen. Es 
ist aber sehr wahrscheinlich, dass diese Liste im Laufe der Zeit noch wesentlich 
erweitert werden muss. Einerseits ist naémlich die Untersuchung vieler Sub- 
_ stanzgruppen, zum Beispiel der intermetallischen Verbindungen, von denen das 
antiferromagnetische CrSb ein Vertreter ist, noch lange nicht abgeschlossen, 
andererseits ist sicher, dass sehr viele der bis heute als normal paramagnetisch 
bekannten Stoffe sich bei sehr tiefen Temperaturen antiferromagnetisch erwei- 
sen, wie zum Beispiel die Alaune. 


Zusammenstellung antiferromagnetischer Substanzen 


ats Vor- 
owl eel Vwi Lee Kristall- vey : : 
Substanz IES O°k ae ran eae ee Literatur 
Sree sree nee. | 420 1,8 [46], [58] 
Riau te 2 aes Ih lO 610 6 NaCl 100 | [8], [63] 
OMe sod 2 tan, il AOS 570 8 NaCl 100 | (8) 
SOO mate otek... F208 280 5,3 ~ |-Nacl 100 | [8], [40], [69] 
OGRE ws ees sft O47 2470 NaCl 100 | (8), [40], [69] 
Biase ces 165 528 6,0 | NaCl [8], [63] 
Masem te t o05 & > 1 247 740 | 19 NaCl [43], [63] 
Neate. Yas d ccm 10307 4,6 | Nids (63) 
Gesieeten. kets oe |) 700: 1,9 | NiAs [62] 
i Oh aire 92 2.8 wel util c a 
Wem. See ic. | 273 1,0 [17] 
ih ie ee hexagonal [18], [19], [30] 
WerOs — 6 vac) tl - 950 2000 hexagonal [11] 
IRXCOh sae. 3 ches Ue 57 14 48 hexagonal c [8] 
PEE gery uid 72 TAS tence Rutil a,b | [8], [26] 
SLE ASI Re RONEN 79 LF \pee22 Rutil a,b | [8] 
GO eee tt ae ame 70 53) 
C6 a aoe 24 ce 
BOCineeer Pe a: 25 9] 
ENG La Mant als 50 
MON 30 [64] 
100 | [53] 
OuCls 2H,0 oa be 4,3 rhomb. Se ae ee 
(SO,),CrK,12H,0 . . 0,004 2-108 | kubisch (22}, (38) 
(SO,).FeNHy,, 12 H,O. 0,043 108 (22), vee 
(SO,)>Mn(NH,).,6 H,O 0,12 oes 
(S0)),Cuks,GH,O . © 0,05 | eh 
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Diese sehr tiefen Curie-Temperaturen riihren méglicherweise von magneto- 
statischer Dipol-Dipol-Wechselwirkung her. 

Neben den reinen Antiferromagneten scheinen auch Ubergangssubstanzen 
zum [errimagnetismus bekannt zu sein, bei denen sich die Magnetisierungen 
der Untergitter im feldlosen Zustand nicht exakt kompensieren. Solche Ver- 
haltnisse glaubt YosrpA [81], [82], [83] in den Systemen FeS,,, CrS,, und bei 
Pr,O,, anzutreffen. 

Den Fall ausgepragten Ferrimagnetismus findet man vor allem bei den im 
Spinellgitter kristallisierenden Ferriten, welche wegen ihrer technischen Bedeu- 
tung als Ferromagnetika mit sehr schlechtem elektrischem Leitvermégen be- 
sonders eingehend untersucht sind. Zusammenfassende Darstellungen hieriiber 
publizierten FAIRWETHER et al. [15], SNOEK [61] und Lasuart [39]. 


Zusammenfassend lasst sich feststellen, dass das Studium des Antiferro- 
magnetismus in den letzten Jahren einen wesentlichen Beitrag zum Verstand- 
nis der magnetischen Eigenschaften fester Korper geliefert und dariiber hinaus 


wertvolle Gesichtspunkte und Anregungen zur Behandlung kooperativer Pha- 
nomene im allgemeinen beigesteuert hat. 


Die Entwicklung befindet sich aber noch in vollem Gange, und die Trag- 


weite der neuen Erkenntnisse fiir die Physik des festen Ko6rpers oder gar fiir die 
Technik ist noch keineswegs abgegrenzt. 
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Energicbedarf zur Verhiitung von Vereisungen 
an Freileitungen’) 


Von THOMAS BRUNNER®), Weissfluhjoch, Davos 


1. Das Problem 


Nachdem MELCHER [1]8) in seiner Arbeit iiber Vereisungserscheinungen die 
Entstehung von Reif und Eis im Windkanal eingehend untersucht hat, sollte 
abgeklart werden, wie gross der Energieaufwand ist, um Vereisungen bei den 
verschiedensten meteorologischen Bedingungen zu verhiiten. 

Es gibt grundsitzlich zwei Méglichkeiten, eine Freileitung durch Warme 
vor gefahrlichen Vereisungen zu schiitzen: Entweder wird die Leitung, sobald 
sich eine gewisse Eismenge gebildet hat, mit einem starken Stromstoss belastet, 
so dass das Eis an der Leiteroberflache schmilzt und abfallt, oder es wird die 
Leitung dauernd etwas iiberbelastet, so dass die erzeugte Warme ausreicht, die 
Eisbildung zu verhiiten. 

Die erste Methode lasst sich nur anwenden, wenn die Leitung nicht unter 
Hochspannung steht; wenn namlich die Eislast plétzlich abfallt, beginnen die — 
Leitungen zu schwingen (galloping conductors [2]), was bei gegenseitiger Be- 
rihrung zu einem Kurzschluss fiihrt. Auch besteht die Gefahr, dass eine ganze 


Reihe von Masten umgelegt wird. Das Problem schwingender Leitungen ist bei 
[2] eingehend behandelt. 


1) Ausgefithrt am Eidgenéssischen Institut fiir Schnee- und Lawinenforschung im Auftrag der 
Schweizerischen Kommission fiir Vereisungsfragen. 


? Siehe Bemerkung ZAMP 8, 460 (1952), Fussnote 1. 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 84. 
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Giinstiger scheint es deshalb, eine Leitung dauernd eisfrei zu halten. Eisfrei 
bedeutet in diesem Zusammenhang nicht vélliges Fehlen von Eis; Spuren von 
Eis schaden ja nichts, falls die Menge im Laufe der Zeit nicht zunimmt. 


2. Erzeugung des Eises 


Die Versuche wurden in einem offenen Windkanal, der in einem KAltelabor 
von —10°C aufgestellt ist, durchgefiihrt. Die Temperatur liess sich in einem 
Bereich von —15° bis —5°, die Windgeschwindigkeit zwischen 0 und 20 m/s 
regulieren. 

Erzeugung von Eis durch Kondensation von Dampf kam fiir diese Versuche 
nicht in Frage, da diese Methode keine schweren Vereisungen zu erzeugen ge- 
stattet. Die Zerstaubung von Wasser nach dem Prinzip der Spritzpistole schien 
uns nicht geeignet, weil eine genaue Justierung der Diisen und somit die Ein- 
haltung einer einheitlichen Tropfchengrésse unméglich ist. Das Hindurchpres- 
sen von Wasser durch enge Bohrungen ergibt zwar feine Wasserstrahlen, aber 
viel zu grosse Trépfchen. Die fiir die folgenden Untersuchungen benutzte Diise 
arbeitet nach dem Prinzip der Selbstzerstaubung (Herrn Prof. Dr. EICHELBERG 
von der ETH. danke ich fiir die Uberlassung einer solchen Diise). Auf einer 
Diisennadel, welche die Form eines Zylinders mit aufgesetztem Kegelstumpf 
hat (Figur 1), sind vier zirka 0,05 mm tiefe Rillen eingeritzt. Setzt man auf 
die Diisennadel einen konischen Ring (Figur 2), so sind Ring und Nadel langs 


Fig. 1 
Diisennadel mit wasserfiihrenden Rillen. 


Fig. 3 
(R), Heizwicklung (H#) und Wattefilter (F) versehenen 
Diisen. 


=, 


Anordnung der mit Reguliervorrichtung 
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des ganzen Umfanges, mit Ausnahme der vier Rillen, dicht. Presst man Wasser 
zwischen Ring und Nadel, so entstehen vier feine Strahlen, die sich etwas ausser- 
halb der Diise treffen. Durch das Aufeinanderprallen der vier Wasserstrahlen 
entstehen feine Wassertropfchen, die gleichmassig in einen kegelf6rmigen Raum 
ausgestrahlt werden. 

Um eine gleichmassige Verteilung der Troépfchen tiber den Kanalquerschnitt 
zu erreichen, wurden zwei Diisen nebeneinander angeordnet (Figur 3). Eine — 
Einrichtung gestattet, Nadel und Ring in einer bestimmten gegenseitigen Lage 
zu fixieren. Auf diese Weise wird erreicht, dass die Diisen bei allen Versuchen 
gleich arbeiten. Ein Wattefilter / halt Verunreinigungen des Wassers, die die 
Rillen verstopfen kénnten, zuriick. Eine Heizwicklung H (30-60 W) verhindert 
das Einfrieren der Diisen. Die Wicklung wurde in eine Talk-Wasserglas-Paste 
eingebettet, die zum Schutze gegen Tropfwasser mit Zaponlack lackiert wurde. 


3. Versuchsk6rper und Apparatur 


Als Versuchskorper fiir die Anreifung dient ein 27 cm langes Aluminium- 
rohr von 10 mm Aussendurchmesser, in dessen innere Bohrung von 6,8 mm ein 
Heizelement gebracht wurde. Dieses besteht aus einem Pyrostearohr, auf das 
0,2 mm dicker, oxydisch isolierter Konstantandraht eng aufgewickelt wurde. 
Der geringe Zwischenraum zwischen Drahtwicklung und Wand des Aluminium- 
rohres wurde mit Wasserglas ausgefiillt, so dass keine warmeisolierende Luft- 
schicht vorhanden ist (Figur 4). 

Die Enden des Heizelementes sind zum Schutze gegen Warmeableitung 
in Glas gefasst und zur Vermeidung der Konvektion mit Glaswolle verstopft. 

Der Stab wird mit Wechselspannung, die mittels eines Variac (Regulier- 
transformer) von 20-150 V reguliert wird, geheizt. Die Netzspannung, mit der 


Glas Al.- Rohr 
Glaswolle Const.—Oraht 
Pyrostearohr 


Fig. 4 


Detail des zylindrischen Versuchskérpers (d = 10 mm). 


der Variac betrieben wird, ist stabilisiert; die bendtigte Heizenergie wird an 
einem Wattmeter abgelesen. 

Auf den Versuchsstab sind zwei Thermoelemente aufgeklemmt, eines in der 
Mitte, eines 4 cm von einem Stabende entfernt. 

Die Windgeschwindigkeit wird mit einem Handanemometer bestimmt, die 
Lufttemperatur mittels eines Quecksilberthermometers, das in den Windkanal 
ragt, ohne von den Trépfchen getroffen zu werden. 
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4. Bestimmung der Trépfchengrésse, des Gehaltes an fliissigem Wasser 
und der Reifrate 


Eine direkte Bestimmung der Tvdp/chengrisse unter dem Mikroskop ist einer 
indirekten Methode vorzuziehen. MELCHER hat die Trépfchen auf einer Gela- 
tineschicht, die er blau einfarbte, aufgefangen und unter dem Mikroskop be- 
trachtet. Diese Methode lefert keine Absolutwerte. Die Schwierigkeit der Er- 
zeugung einer gleichmassigen Gelatineschicht und einer giinstigen Einfarbung 
fallt weg, wenn man Photoplatten verwendet, auf die ein leichter Schleier ent- 
wickelt wird. 

Die Bestimmung des Gehaltes an fliissigem Wasser wurde nicht wie von 
MELCHER einfach aus dem Verbrauch des Wassers berechnet. Tropfwasser an 
den Diisen sowie Reifbildung an den Wanden des Windkanals wiirden Fehler 
verursachen. Es wurden deshalb 5 bzw. 10 cm oberhalb des heizbaren Probe- 
stabes zwei weitere, aber ungeheizte Aluminiumstabe in den Windkanal ge- 
bracht. Die wahrend der Versuchsdauer auf ihnen abgelagerte Eismenge ist ein 
Mass fiir den Gehalt an fliissigem Wasser G}°°. Dabei nimmt man an, dass alle 
Trépfchen, die im angeblasenen Querschnitt vorhanden sind, sich auf dem 
Stab ablagern. Diese Annahme ist gemass der Theorie von LANGMUIR [3] unter 
den herrschenden Bedingungen erfiillt. 

Als Reifrate bezeichnen wir die pro Quadratzentimeter und Stunde abge- 
lagerte Eismenge in Gramm. R, ist dabei die Reifrate, wie sie an einem Stab 
ermittelt wird, dessen Querschnitt sich wahrend der Eisablagerung nicht ver- 
grossert hat, R diejenige eines Stabes, dessen Querschnitt im Laufe der Zeit 
zagenommen hat. 

Der Zusammenhang zwischen R, und R hangt stark von den Versuchsbe- 
dingungen ab. 

Bei der Bestimmung der Reifrate ist noch zu beriicksichtigen, dass die Ver- 
teilung der Wassertrépfchen tiber den Querschnitt des Windkanals nicht immer 
gleichmassig ist. Die grésseren Trépfchen haben namlich bei kleinen Wind- 
geschwindigkeiten Zeit, emen gewissen Teil der Héhe des Kanales zu durch- 
fallen, wodurch dann die Trépfchendichte oben klein und unten gross wird. 
Durch Vorversuche kann diese Verteilung fiir verschiedene Windgeschwindig- 
keiten gemessen und als Korrektur bei den eigentlichen Versuchen angebracht 
werden. 

Der Gehalt an fliissigem Wasser pro Kubikmeter nimmt — bei konstanter 
Leistung der Diisen — mit zunehmender Windgeschwindigkeit ab. Wenn der 
Wassergehalt fiir verschiedene Windgeschwindigkeiten konstant gehalten wer- 
den musste, wurden bei kleinen Windgeschwindigkeiten die Diisen wahrend 
kiirzerer oder langerer Zeit periodisch mittels eines Schiebers zugedeckt. Der 
Wassergehalt schwankt daher. Da aber das ganze Heizsystem eine gewisse 
Tragheit hat, werden dadurch die Resultate nicht verfalscht. Eine andere Re- 
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gulierung des Wassergehaltes kam nicht in Frage, denn zum Beispiel eine Er- 
héhung des Wasserdruckes oder eine Vergrésserung der Rillen in der Diise hatte 
sofort auch eine Veranderung der Troépfchengrésse und -verteilung zur Folge. 


5. Der Warmeiibergang an einem trockenen, zylindrischen Stab 


Die Heizleistung wird bestimmt, um den Stab bei verschiedenen Windge- 
schwindigkeiten auf einer bestimmten Temperatur zu halten (zum Beispiel 
15°C iiber Labortemperatur). Die Warmeiibergangszahl « ist definiert durch: 


W,=aFAT, (1) 


wo W, die pro Zeiteinheit abgefiihrte Warmemenge, F die Oberflache des zylin- 
drischen Stabes und AT die Temperaturdifferenz bedeutet. 

NusSsELT [4] bekommt fiir die Warmeiibergangszahl eines trockenen, zylin- 
drischen Stabes, der senkrecht angeblasen wird, folgenden Ausdruck: 


a = 0,067 A,, (1273 + Re)®7® cal/em? Grads (2) 
Am bedeutet die mittlere Leitfahigkeit der Grenzschicht, fiir Luft: 
5,7 + 10-* cal/s cm Grad. 


Re bedeutet die Reynoldssche Zahl. 


Vdo 
Kee on (3) 
wobei V = Windgeschwindigkeit, d = Durchmesser des Stabes, 9 = Luftdichte, 
» = Zahigkeit der Luft bedeuten. 
Tabelle 1 gibt die nach NussELT berechneten Werte der Warmeiibergangs- 
zahl, berechnet fiir @ = 0,97 kg/m*, 1 = 166+ 10-® g cm-! s~1 sowie die notige 
Heizleistung W,), um einen zylindrischen Stab von 1cm Durchmesser und 


100 cm Lange auf einer Temperatur von 10°C iiber der Lufttemperatur zu 
halten. 


Tabelle 1 


, _ _Wéarmeiibergangszahl und Warmeverlust 
berechnet fiir zylindrischen Stab (1 cm Durchmesser, 100 cm Lange, AT = 10°) 


V [m/s] 


aS eee eee 5850 8800 11700 


| nes fe —————E—EE————EEEE 


6,35 9,05 LL 14,1 


| Rees) en ee ee ee 


Se Cer ‘ 20,0 


a 
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Die graphische Darstellung (Figur 5) zeigt die berechneten und die ge- 
messenen Werte der Warmeiibergangszahl «. 


af 10° 
ro 


? Grad, 


Watt 


ie 
= 
So 


Zp. Nusselt 


on 


So 


te) 5 10 15 20 


Warmeubergangszahl « 


vy m[sek 
Fig. 5 

Warmetibergangszahl ~ eines zylindrischen Stabes von 1 cm Durchmesser in Funktion der Wind- 

geschwindigkeit v. 


Die Abweichung zwischen gemessenen und berechneten Werten mag von 
der Schwierigkeit herriihren, die Oberflachentemperatur zu bestimmen. 

Die Warmeiibergangszahl bei der Windstarke 0 m/s (Windstille) enthalt den 
Einfluss der freien Konvektion, der Strahlung und der Ableitung. 


6. Eisfreihaltung eines Stabes 


Zuerst wird der Stab so stark geheizt, dass sich auf ihm nur Tropfen und keine 
Spuren von Eis zeigen; die aufgewendete Heizleistung bedeutet eine obere 
Grenze der zur Eisfreihaltung benétigten Energie. Dann wird die Heizung gedros- 
selt, bis alles Wasser gefriert, wodurch sich eine untere Grenzenergie ergibt. 
Dann wird die Heizleistung innerhalb dieser Grenzen so lange variiert, bis sich 
ein Gleichgewicht zwischen Wasser und Eis auf dem Stab bildet, wobei nur eine 
geringe Eismenge, die sich nicht vermehren darf, zugelassen wird. Die Feststel- 
* lung, wann dieses Gleichgewicht erreicht ist, ist natiirlich in gewissem Masse 
‘eine Ermessensfrage. Immerhin scheint eine Genauigkeit von + 10% méglich. 

Im folgenden sollen nun die Einfliisse der verschiedenen meteorologischen 
Faktoren auf den Energiebedarf untersucht werden. 


a) Verdunstung 


Der zusitzliche Warmebedarf, den ein feuchter Stab gegeniiber einem 
trockenen benotigt, rithrt in der Hauptsache von der Verdunstung her. Diese 
wiederum ist abhaingig von der Dampfdruckdifferenz zwischen Stab und um- 
gebender Luft sowie von der Windgeschwindigkeit. Man kann analog zum 
Warmeiibergang einen Massentibergang definieren. Die fiir die Verdunstung 
aufzuwendende Leistung W; ist proportional der Dampfdruckdifferenz Ae und 
der Flache F, an der Verdunstung stattfindet: 


W,=BeFAe ; (4) 
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f ist dabei die «Masseniibergangszahl». ¢ gibt an, welcher Bruchteil der Flache F 
feucht ist. Bei kleinen Windgeschwindigkeiten ist der Stab nur auf der Vorder- 
seite feucht, somit ¢e ~ 0,5. Je grésser die Windstarke ist, um so mehr Trépfchen 
fliessen oben iiber den Stab auf die Riickseite ; e nimmt mit zunehmender Wind- 
starke zu. 

Der gesamte Energiebedarf des feuchten Stabes betragt unter Beriicksichti- 
gung von Gleichung (1) 


W=W,4+W,=aFAT+ BeFAe. (5) 
Nach der Ahnlichkeitstheorie von NussELT [5], die den Zusammenhang zwi- 
schen Warmeiibergang und Verdunstung angibt, ist 


pe REO 
a: Mae ke (6) 
wo k der Diffusionskoeffizient Wasserdampf/Luft (0,28 cm? s~! bei 0°C), 
y die Verdampfungswarme (596 cal g~! bei 0°C), 


zder Umrechnungsfaktor absolute Feuchte in Millimeter Dampfdruck * 


(1,05 - 10-§ g cm-* (mm Hg)? bei 0°C), 
A die Leitfahigkeit der Luft. 
C hat bei 0°C den Wert 2,8. 
Setzt man die Beziehung (6) in (5) ein, folgt 


W=aF (AT+CeAe). (7) 


Zur Priifung der Theorie wurde aus gemessenen Warmeverlusten unter Ver- 
wendung berechneter Werte fiir « und C mit Hilfe von Gleichung (7) das Pro- 
dukt ¢ Ae ermittelt und mit Dampfdruckwerten verglichen, die auf anderem 
Weg bestimmt wurden. 


Das in Figur 6 wiedergegebene Diagramm zeigt ¢ Ae als Funktion der Wind-— 


geschwindigkeit. Mit steigender Geschwindigkeit nihern sich die Funktions- 
werte asymptotisch dem Dampfdruckunterschied Ae, der sich aus der Tempe- 
ratur und Feuchtigkeit des Laboratoriums errechnet, was wohl als Bestatigung 
der Theorie gewertet werden darf. Das Abfallen der Kurven gegen kleinere 
Windgeschwindigkeiten geht zu Lasten der Benetzung, die, wie bereits erwahnt, 
bei schnellen Strémungen den Zylinder ganz umschliesst (e = 1,0), bei kleineren 
Geschwindigkeiten dagegen die Riickseite trocken lisst. Bei der Ungenauigkeit 
der Feuchtemessung und der Unsicherheit beziiglich « kann keine allzu hohe 
Genauigkeit erwartet werden, 


Vorderhand unerklirlich ist der Verlauf der Kurve fiir — 15° oberhalb einer 
Windgeschwindigkeit von 10 m/s. 


Nachdem sich, von der erwihnten Ausnahme abgesehen, gute Ubereinstim- 


mung zwischen Theorie und Messung ergeben hat, lasst sich der Energiebedarf — 


ese 
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eines feuchten Stabes in Abhangigkeit von T und Ae mit Gleichung (5) berech- 
nen. Dabei soll die relative Feuchtigkeit zu 100% angenommen werden. Fiir 
den speziellen Fall eines Stabes von 100 cm Lange, 1 cm Durchmesser und 0° 
Oberflachentemperatur sei die Leistung mit W, bezeichnet. 


150415 


o n.Nusselt berechnet 
x gemessen 


Wo 


&-fe- mmHg—— 


Werte von ¢ Ae, berechnet gemass 


- ‘Fig. 6 


Gleichung (7) in Funktion der 


Windgeschwindigkeit. 

é = Bruchteil der benetzten Flache. 
Ae = Sattigungsdefizit der anstré- 
menden Luft, bezogen auf die 
Temperatur des Versuchsstabes 


(OMG): 


vo m/sek ——> 
Fig. 7 
Energieverlust W ) eines ringsum feuchten 
Stabes von 1 cm Durchmesser und 0° C Tem- 
peratur, der von einer seitlichen Strémung 
gesattigter Luft verschiedener Temperatur 
und Windstarke angeblasen wird. 

Linke Ordinatenskala: bezogen auf eine Stab- 
lange von 100 cm; rechte Ordinatenskala: 

bezogen auf 1 cm? projizierter Flache. 


Tabelle 2 


Energiebedarf W, (W) eines feuchten, auf 0° C gehaltenen Stabes 
(Luftfeuchtigkeit 100% ; Stablange 100 cm, Durchmesser 1 cm, ¢ = 1,0) 


Windgeschwindigkeit in Meter pro Sekunde 


1 2 5 10 | 15 | 20 

= 5°C, 

Ae = 1,57 mm Hg 10,3 12,6 18,7 26,7 34,5 41,7 

=? Cy 

Ae = 2,63 mm Hg 19,0 23;3 34,6 49,4 63,8 viet 
eae 15°C, 

Ae = 3,34 mm Hg 26,6 32,6 48,6 69,3 89,4 108,1 
fa 202'G, 

Ae = 3,81 mm Hg Sey 7/ 41,0 612 87,3 | 112,6 | 136,3 


eT E 
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b) Reifrate 


Da die Nebeltrépfchen unterkiihlt sind, ist ein gewisser Energiebetrag notig, 
um sie auf 0°C zu erwarmen. Die benétigte Leistung betragt bei einer extrem 
grossen Reifrate von 4 g cm~* h~! und einer Unterkiithlung von —10°C 0,046 W 
pro Quadratzentimeter projizierter Flache oder 4,6 W fiir einen Stab von 
einem Meter Linge. Dieser Betrag fallt mit zanehmender Windgeschwindigkeit 
immer weniger ins Gewicht, da der Betrag W, im allgemeinen bedeutend 
grosser ist. 

Die gleiche Uberlegung gilt natiirlich auch in bezug auf den in der Luft 
enthaltenen Niederschlag (Schnee). Ein extrem starker Schneefall mag eine 
Wassermenge von 10 mm/h ergeben, entsprechend 1 g/cm? je Stunde. Diese 
Wassermenge betragt also nur einen Viertel der oben angenommenen Reifrate 
und hat deshalb ebenfalls keinen grossen Einfluss. 


c) Gehalt an fliissigem Wasser 


Viele Messungen haben ergeben, dass der Gehalt an fliissigem Wasser G10° 
keinen oder héchstens einen ganz geringen Einfluss auf den Energiebedarf 
hat. Sicher ist der Einfluss einer Windgeschwindigkeitsanderung um + 1 m/s 
grdésser als derjenige, der durch eine Veranderung von G3°° um + 0,5 gm-$ 
bedingt ist (der in der Natur maximal vorkommende Wassergehalt diirfte 
1 gm? nur unwesentlich iibersteigen). 


d) Tropfchengrésse 


Die Abhangigkeit des Energiebedarfes in Funktion der Tropfengrésse wurde 
nicht untersucht. Aus den Untersuchungen von SCHAFER auf Mount Washing- 
ton [6] weiss man aber, dass die Trépfchengrésse einen betrachtlichen Einfluss 
auf den Energiebedarf hat. Bei einem Trépfchendurchmesser von 12 ist — bei 


sonst ahnlichen Bedingungen — der Energiebedarf etwa doppelt so gross wie 
bei Trépfchen von 6 «4 Durchmesser. 


e) Anblasrichtung 


Bei unseren Untersuchungen stand der Stab stets senkrecht zum Wind 


(Kreuzstrom). EpER [7] untersuchte den Zusammenhang zwischen Energie- 
bedarf und Anblasrichtung #. 


Die Abhangigkeit wird wie folgt ausgedriickt : 


W(8) = W f(d) , 


-— 


; aoe en ae aes he 
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wo / folgende Werte hat: 


o f 
90° (senkrecht) . 1,00 
COs et ae eee 0,90 
SLO” wey oR aE BON 0,68 
O” Goeveewilel) 2 ¢ 0,52 


Der Warmeverlust eines parallel angestro6mten Drahtes ist also nur rund 
halb so gross wie derjenige eines senkrecht angestromten. 


7. Vergleich mit den Resultaten von Schafer 


ScHAFER hat im Jahre 1946 Vereisungserscheinungen wahrend der Stiirme 
auf Mount Washington untersucht. Der vertikal aufgestellte, heizbare Probe- 
stab war 40 cm lang und hatte einen Durchmesser von 8 mm. SCHAFER kommt 
zum Schluss, dass 10 W/cm? projizierter Flache geniigen, um diesen Stab bei 
extremen Bedingungen trocken, 3 W/cm? projizierter Flache, um ihn eisfrei 
zu halten. 

Berechnet man den Warmebedarf an Hand unserer Resultate fiir die von 
ScHAFER beobachteten meteorologischen Bedingungen, so findet man nur 
etwa die Halfte der von SCHAFER angegebenen 3 W/cm?. 

Die Erklairung dieser Diskrepanz fallt schwer. Es ist méglich, dass SCHAFER 
sozusagen kein Eis auf dem Probestab zugelassen hat, was den Warmebedarf 
etwas vergréssern wiirde. 

Wie bereits erwahnt, wird auch die Trépfchengrésse eine Rolle spielen. Bei 
unseren Untersuchungen betrug sie 0,1 bis 0,01 mm, wahrend sie bei SCHAFER 
rund eine Zehnerpotenz tiefer war und eigentlich noch viel kleinere Energie- 
verluste hitte bewirken sollen. Es ware nun denkbar, dass sich die von SCHA- 
FER festgestellte Zunahme des Energiebedarfes mit der Trépfchengrésse nicht 
monoton fortsetzt, sondern von einer gewissen Grosse an durch eine Abnahme 
abgelést wird, und dass die von ihm untersuchten Trépfchen gerade im Grdssen- 
bereich des Maximums lagen. Der Entscheid muss vorderhand offengelassen 


werden. 


8. Schlussfolgerung 


Die Untersuchungen haben ergeben, dass auch bei extremen meteorologi- 
schen Bedingungen 


T = —20°C; v = 20 m/s; 10 mm Niederschlag pro Stunde 
und einer Reifrate von 4g cm~* h~? 


ZAMP. IV/3 
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mit einer Heizleistung von 1,5 W/cm? projizierter Flache ein zylindrischer Draht 
von 1 cm Durchmesser eisfrei gehalten werden kann. Dieses Resultat lasst sich 
nicht ohne weiteres auf Drahte anderer Dicken oder andere meteorologische 
Bedingungen tibertragen, da sich die Reynoldssche Zahl andert und die Ver- 
dunstung nicht linear in den Ausdruck fiir die Heizleistung eingeht. 
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Summary 


Ice was produced in a wind tunnel at a laboratory temperature of —10° C by 
spraying water through a nozzle. An electrically heated cylindrical rod, 30 cm 
long, 1cm in diameter, is exposed to various meteorological conditions. The 
amount of energy needed to keep the rod free of ice is investigated. At a tempera- 
ture of —15°C and a wind velocity of 20 m/s, 1,1 W/cm? of projected area are 
sufficient to keep the rod clear. The result agrees with the theory of NUSSELT. 


(Eingegangen: 3. August 1952.) 
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Beitrage zur Kenntnis 
des Biorthogonalisierungs-Algorithmus von Lanczos 


Von HEINZ RUTISHAUSER, Ziirich!) 


§ 1. Problemstellung 


Unter einer Kodiagonalmatrix verstehen wir eine Matrix A mit der Eigen- 
schaft a;; = 0 fiir |¢ — j| >1. Eine solche Matrix hat also nur in der Haupt- 
diagonalen und in den beiden anliegenden « Kodiagonalen» nichtverschwindende 
Elemente: 


Vy Wy 
Us Ug a 2 
Ug Vg Ws 
(1) 


Es wird sich zeigen, dass man jede Matrix auf Kodiagonalform transfor- 
mieren kann, sogar auf unendlich viele Arten. Die praktische Durchfthrung 
einer solchen Kodiagonaltransformation, welche insbesondere fiir die Eigen- 
wertbestimmung in letzter Zeit eine grosse Bedeutung erlangt hat, erfolgt 
zweckmassig nach einem von Lanczos [1]*) angegebenen Verfahren oder mit 
dem n-Schritt-Verfahren von E. STIEFEL und M. R. HEstTENEs [2], [3]. Im 
Prinzip beruhen beide Verfahren auf der Orthogonalisation einer Vektorfolge 
x, Ax, A?x, ..., bzw. bei nichtsymmetrischen Matrizen auf der Biorthogonali- 
sation von zwei Vektorfolgen-x, Ax, A?¥%,... und y, A*y, A**y, ..., wobel 
A* die Transponierte von A bedeutet. 

Die Kodiagonalform einer Matrix ist nicht eindeutig bestimmt, sondern 
hangt bei gegebener Matrix A noch weitgehend von den bis zu einem gewissen 
Grade willkiirlich wahlbaren Vektoren x und y ab. 


1) Institut fiir angewandte Mathematik an der ETH. 


2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis anf Seite 56. 
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Diese Arbeit befasst sich nun mit der Abhangigkeit der Kodiagonalform von 
den Anfangsvektoren x, y, wortiber zur Zeit wenig bekannt ist (vgl. jedoch [4] 
fiir symmetrische Matrizen). Insbesondere wird auch untersucht, wieweit es 
méglich ist, die Anfangsvektoren x und y so zu wahlen, dass die Kodiagonal- 
elemente, das sind die uw; und w; in (1), klein werden, denn diese Eigenschaft ist 
fiir viele Anwendungen sehr erwiinscht. Fiir symmetrische Matrizen ist dies bei 
[4] bereits untersucht. Die vorliegende Arbeit befasst sich aber auch mit nicht- 
symmetrischen Matrizen und ergibt folgende Resultate: 


a) Fiir Matrizen mit ausschliesslich reellen Eigenwerten (auch bei nichtlinearen 
Elementarteilern) kann man die Anfangsvektoren x, y tatsachlich immer so 
wahlen, dass simtliche Kodiagonalelemente beliebig klein werden. 

b) Dies gilt nicht mehr fiir eine Matrix mit komplexen Eigenwerten. Man kann 
aber in diesem Fall erreichen, dass fiir jedes Paar komplex konjugierter 
Eigenwerte A -- 7 w zwei Kodiagonalelemente w;.,,, w, auftreten, so dass die 
diagonale Untermatrix 


Unis Ue 


von (1) ungefahr die Eigenwerte A -+- iu hat und ferner alle Kodiagonalele- 

mente, die nicht in einer dieser Untermatrizen enthalten sind, beliebig klein 

werden. 

Die theoretischen Untersuchungen des Verfassers erstrecken sich im Fall 
komplexer Eigenwerte nur auf lineare Elementarteiler, doch lassen numeri- 
sche Versuche vermuten, dass dieselbe Aussage auch fiir Matrizen mit nicht- 
linearen komplexen Elementarteilern richtig ist, sofern man die Eigenwerte 
mit Vielfachheit zahlt. 

Die in dieser Arbeit nur fiir endliche Matrizen durchgefiihrten Untersuchun- 
gen lassen sich sinngemass auch auf unendliche Matrizen sowie Differential- und 
Integraloperatoren iibertragen (vgl. die Originalarbeit von C. Lanczos [1)). 
Allerdings erfahren dann einige der hier erhaltenen Aussagen eine Abschwa- 
chung; beispielsweise wird man im allgemeinen nicht alle unendlich vielen 
Kodiagonalelemente £,, gleichmdssig klein machen kénnen. 


§ 2. Der Biorthogonalisierungs -Algorithmus!) 


Sei A die gegebene N-reihige Matrix und seien x und y zwei Vektoren, die 
zunachst der einzigen Bedingung (x, y) + 0 unterliegen sollen. Dann bestimmt 
man, von %,= * und y, = y ausgehend, mit nachstehender Vorschrift die Vek- 


1) Fortan kurz mit «der Algorithmus» bezeichnet. 
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ROVOML DX, Gey Ag WO 4/7, Vo, ate Ves 
Ri Ak 0g My Vo A yy Vy Mit. ay = — 
Wanner Wine ——2 no eee wiVe 


Huy = A xy — Oy Xy = Bua Miri Mie Ve OE Vi Bua Vet (3) 
mit 
(A Xp Vx) 


(Aare Sree) 
(ip Va) Pra = aay (4) 


Cp 
(%p1) Ves) 


Dieser Algorithmus liefert nach der allgemeinen Theorie!) zwei biorthogo- 
nale Vektorsysteme *;, v;,, das heisst zwei Vektorsysteme mit der Eigenschaft 


(%j, Ye) =O (fir 7 + R), (5) 


die, falls fiir alle k (xz, v,) + 0 ist, je eine Basis des Raumes bilden. Der Algo- 
rithmus kann auf Grund der Formeln (4) genau so lange durchgefiihrt werden, 
als noch (%,, V,) + 0 ist, und kommt spatestens fiir k = N zu einem Abschluss, 
da ja mehr als N Vektorpaare nicht biorthogonal sein kénnen; die Formeln (3) 
miissen dann %y,1 = Yy.1= 0 liefern. Das vorzeitige Versagen der Formeln 
(3), (4) wird in § 4 behandelt. 

Zwischen den Vektoren x;, y, und den Skalaren «,, 6;, bestehen noch eine 
Reihe weiterer Beziehungen, die man durch elementare Umformungen aus 
(2), (3) und (4) erhalt und von denen wir hier ohne Beweis einige angeben: 


(Va, A Xp1) = (Xe, Vu) = (A Xp, Vicaits (6) 

— (%p, Vx) a] 

Pro mee ee ( ) 

(A'x;, Vu) = (yj, A* ;) = 0° ftir oy — RS 1. (8) 


§ 3. Vereinfachungen 


Wenn die Matrix A symmetrisch ist, ist es zweckmiassig, die beiden Anfangs- 
vektoren einander gleichzusetzen, es ist dann immer %;, = Yz, $0 dass man tiber- 
haupt nur mit einem Vektorsystem arbeiten muss, welches dann ein Orthogo- 
nalsystem ist. 

Aber auch wenn A nicht symmetrisch ist, gelingt es oft, gewisse Figen- 
schaften auszuniitzen, so dass man ebenfalls mit einem Vektorsystem aus- 


1) Vergleiche die Originalarbeit von C. Lanczos [1]. 
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kommt. Dies ist beispielsweise der Fall, wenn eine Matrix R mit der Eigenschaft 
A* — RA R~11) gefunden werden kann; man unterwirft dann die Anfangs- 
vektoren der Bedingung y = R x. Es wird dann 

Vg = A* yy, — 0, 9,= RAR, —, y= RAY —o, Rx 


= K (Ax, — 0) = ae 
und ebenso zeigt man 
Vu = RR ae CU ee) oe eh (9) 


Man kann also auch in diesem Fall mit den Vektoren x, allein arbeiten. 


§ 4. Vorzeitiges Abbrechen des Algorithmus 


Wie bereits erwahnt, muss der Algorithmus nach spatestens N Schritten 
zum Abschluss kommen. Unter gewissen Umstanden kann er aber schon friiher 
abbrechen, indem fiir ein gewisses k < N einer der drei nachstehenden Faille 
eintritt : 

a) Man erhalt mit Formel (3) ¥444 = Yeu = 0, 
b) Es verschwindet genau einer der Vektoren *;..4, Vgsy- 
c) Es ist x41 +0, Very +0, aber (%p41, Ves) = 0. 

Im Fall a) sind die von den Vektoren 1, ..., x, bzw. Vy) +++) Vp aufgespann- 
ten Unterréiume invariant in bezug auf die Matrizen A bzw. A*. Diese Erschei- 
nung kann auf ungliickliche Wahl der Anfangsvektoren zuriickzufiihren sein: 
es gibt aber auch Matrizen A, fiir die bei jeder Wahl von x, y schon fiir k << N 
einer der drei Fille a), b), c) eintreten muss. Wenn namlich der Grad m des 
Minimalpolynoms der Matrix A nicht den Wert N erreicht, so sind die Vek- 
toren x, A x,...,A™ x linear abhangig und ebenso die y, A* y, Atty, -.. Atma 
so dass der Algorithmus notwendigerweise ~,, 41 = YVmi1 = 0 ergeben muss, wenn 
nicht schon vorher eine Stérung eintritt. 

Fall b) ist auf ungliickliche Wahl der Anfangsvektoren «, y zuriickzufiihren ; 
wiirde man nimlich x, y der Bedingung y = Rx unterwerfen, wobei R die in 
§ 3 erwihnte nichtsingulire Matrix bedeutet, so kénnte Fall b) wegen (9) 
niemals eintreten. Aber auch Fall c) kann durch geeignete Wahl der Anfangs- 
vektoren x, y vermieden werden, es gilt namlich: 

Satz 7: Hat das Minimalpolynom der Matrix A den Grad m, So ist es immer 
moglich, die Anfangsvektoren x, y so zu wihlen, dass der Algorithmus bis zu 


seinem natiirlichen Ende (k = m) durchgefiihrt werden kann und auch noch 
(Xm, Vm) + 0 ist. 


) Es existiert immer eine solche nichtsingulare (sogar symmetrische) Matrix R (vgl. etwa [5], 
S. 44), Sind alle Eigenwerte reell und die Elementarteiler linear, so ist R sogar positiv definit. 
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Beweis: Auf Grund elementarer Sitze iiber Determinanten und der For- 
meln (2), (3), (5) ergibt sich 


(x9) (asi) ae (ge A haere) 
(Ae) Pa Rae ee (Ae A Fee? 4) 
aeomeerVmed rt Ala 2) LARS A ¥ea1 y) 
(%, 41) (1, Vx) 
(%, 1) (%2, Vx) : 
=> aa] LL 3) 
alia) 9] | fel elfo) ol ielveliie: (el © (e! 6!) 0) \e) louie q=1 
(¥n, V1) (Xn» Vu) 
Somit gilt 
G, 
ky eg 10 
(x Vu) = Gin ( ) 


wobei G,, die links oben stehende Gramsche Determinante der beiden Vektor- 
Pcie Ay, aed UNG, Ve Aye, AY P44 Ist. 
Gemiass (10) geht also unsere Behauptung dahin, dass bei geeigneter Wahl 


von x und y keine der Gréssen Gy, Gy, ..., Gm verschwindet. Fiir den Beweis 
betrachten wir eine einparametrige Schar von Anfangsvektoren 
x(t) = eA? x, y(t) = 64"? yy, (11) 


die auch fiir die weitere Theorie (insbesondere § 8 und § 9) von ausschlaggeben- 
der Bedeutung ist. Es wird dann 


(At x, A*” 9) = (AMT 64? to, 9); 


das letztere ist aber gleich der (w+ v)-ten Ableitung der Funktion fQ) = (41%, Vo)- 
Somit wird G,(¢) eine Wronskische Determinante: 


aN ae ae es 


ia iim f{®) 
G,,(t) = 
Ce an poe) | 
Es geniigt offenbar, wenn man zeigen kann, dass kein G,(t) fir k= 1,0, 


identisch verschwinden kann. 
Wiirde erstmals G,,, identisch verschwinden, so misste f einer linearen 
homogenen Differentialgleichung der Ordnung s mit konstanten Koeffizienten 
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gentigen!). Da bei geeigneter Wahl von %p, Vo sicher G, = f nicht identisch 0 
ist, braucht man nur noch zu zeigen, dass / keiner solchen Differentialgleichung 
mit s < m geniigen kann. p 

In der Tat, wenn das Minimalpolynom der Matrix A die Gestalt IT (~—A,)** 


hat und wenn die beiden Vektoren %9, Vo allgemeine Lage haben, das heisst zu 
keinem Eigenvektor oder Hauptvektor von A* bzw. A orthogonal sind, so ent- 
halt f unter anderem die Summanden 


a (2 
> mit ) My, =m 
CuO uO NC ai 1 
1 ent 


und kann daher keiner Differentialgleichung der Ordnung s < m geniigen; 
was zu beweisen war. 

Wenn m der Grad des Minimalpolynoms ist, so bewirkt also der Algorithmus 
die Abspaltung eines m-dimensionalen invarianten Unterraums R;,. In bezug 
auf diesen Unterraum ist m/(x) gleichzeitig Minimalpolynom und charakteristi- 
sches Polynom von A. 

Man kann nun den Algorithmus im N-m-dimensionalen Komplementar- 
raum von k, fortsetzen [sofern der Algorithmus nicht etwa nach Fall b) oder 
c) abbricht], indem man neue Anfangsvektoren x,» 41 Vas Ves o< ee 
Vm41 L %1) %y, +++) %m Wahlt und mit diesen an Stelle von x,, y, den Algorithmus 
neu beginnt. Auf diese Weise spaltet man einen weiteren invarianten Unter- 
raum FR, ab, wobei das charakteristische und das Minimalpolynom von A in 
bezug auf R, ebenfalls iibereinstimmen usw. 

Auf Grund dieser Betrachtungen kénnen wir uns fortan auf solche N -reihige 
Matrizen beschranken, fiir die das Minimalpolynom den Grad N hat und also 
mit (x) tibereinstimmt. 


§ 5. Die Kodiagonalform 


Die durch den Algorithmus — sofern dieser tiberhaupt bis zu seinem natiir- 
lichen Ende durchgefiihrt werden kann — gelieferten Werte «,, 8, bestimmen 
nun die durch die Matrix A vermittelte lineare Abbildung z > A z vollstandig. 
Mit dem Ansatz 


Ie Do be Xy, Az= Dy, Ne X~ (= Aes A Xxx) 
erhalt man naémlich durch innere Multiplikation mit y, und Beriicksichtigung 


1) Fiir den Beweis siche etwa [Ol S.S2 7k 


. ee : : ‘ 38 
Vol. IV, 1953 Beitrage zur Kenntnis des Biorthogonalisierungs-Algorithmus von Lanczos 41 


der Beziehungen (2), (4), (5), (6) und (8) sowie der Tatsache, dass xy ,;= yy .,—0 
ist, sofort: . oy a 


=> ‘ec y) es ee & = 
My = % 6 + Br és, Wi pay GE Pies kot + Nn = yc, + Oy En 


Sh) 


Die Matrix A hat also im Koordinatensystem der x; die folgende Gestalt: 


Hy By 
| Xo Bo 
Wares ps | 
| : ; } \° (12) 
-: LMG Ae Bay } 
1 An | 


Falls man den Algorithmus in mehreren Stufen ausfithren muss, weil das Mini- 
malpolynom den(Grad m < N hat, so bewirkt dies den Zerfall der Matrix (12) 
in Hauptuntermatrizen; zum Beispiel: 


Ly By 


1 &% Bo 
[et 0 


Nachdem die Matrix A einmal auf Kodiagonalform transformiert ist, kann 
man verhialtnismassig leicht die Eigenwerte berechnen. Dank der einfachen 
Gestalt der Matrix (12) reduziert sich namlich die Berechnung des charakte- 
ristischen Polynoms (A) nach den bekannten Determinantenregeln auf folgende 
Rekursionsformeln : 

Berechne, von #5 = 0, £; = 1 Hae fire = Ly 24 og Ne 


PrialA )= (A 4) Pi lA = By. Pr _1(A) . (13) 


Dann ist 7(A) = pyi1(A). 
Man erhialt so auch die Eigenvektoren; es ist namlich 
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ein Eigenvektor, wenn / ein Eigenwert ist, denn es ist A v — A v auf Grund von 
(13) unter Beriicksichtigung von (2), (3), (6) = —7(A) xy/(xy, Vy) - 

Es bestehen aber noch andere Méglichkeiten zur Eigenwertbestimmung, die 
damit zusammenhangen, dass das Gleichungssystem (A — A £) € = 0 in der 
Kodiagonalform dreigliedrig ist und damit in bekannter Weise in eine Ketten-_ 
bruchgleichung fiir 2 umgewandelt werden kann; die Bestimmung der Eigen- 
werte erfolgt dann auf ahnliche Weise wie bei der Mathieuschen Differential- 
gleichung (vgl. hierzu [6]). 


§ 6. Zur numerischen Durchfiihrung des Algorithmus 


Im allgemeinen wird man beobachten, dass die Gréssenordnung der ite-— 
rierten Vektoren x,, y, mit fortschreitendem k rasch unbequem wird; es ist 
deshalb praktisch, den Algorithmus fiir numerische Zwecke wie folgt abzu- 
andern : . 

(4%, 44) 


So %g = Ax, — Oy %, Se ¥p= A* y,— a y,, mit Oy (x, V1) 5 ght 


“3 —_ me 3 , == *, 
Spg1 Veg = A Xy = Oy Xp — Ve1w%k-1>» Shit Vew = A* Ve— On Vg — Vr-1 Vir (3’) 


mit 
(A Xp, Vp) (A vp, Vea) 


ir (pr Ve)” face (%e-a Vea) © (*/) 

Dabei werden die frei verfiigbaren Skalare Sy, Sg, ..., Sy SO gewahlt, dass die 
Vektoren x;, y, eine numerisch vorteilhafte Gréssenordnung annehmen. Man 
kann zeigen, dass diese Skalare die Rechnung nicht wesentlich beeinflussen: 

Satz 2: Die Gréssen «, und By; = x1 S$, hangen nicht von der Wahl der 
LANIET.S5,'Sa.s nes Sy BD: 

Auf Grund dieser Invarianz stimmen die f; mit den frither definierten 
Werten ;, tiberein, die dem Spezialfall s; = 1 entsprechen. 

Beweis: Die s, beeinflussen nicht die Richtung, sondern nur die Lange der 
Vektoren x;, y,. Gerade a, ist aber gemiss (4’) von der Linge der Vektoren 
Xp» Ve unabhangig. Fiir den Beweis der Invarianz von Yr-1 8, bendtigen wir die 
Relationen, die sich auf Grund der abgednderten Formeln (2’), (3’), (4’) an 
Stelle der Relationen (6) und (7) ergeben?): 


(Vics A Xp_3) = Sz (Xp, Vr) Eee (A Vey Yu-a) , (6’) 
(Vp, Ve) 

Jy, 4 == S$}, . 4 

Ve yA Vee) i) 


1) Die Relationen (5) und (8) bleiben unverandert. 
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Aus (6’) und (7’) erhalt man sofort 


sy A 2 
Bua == eM ( (re 4 Xn) ’ (14) 


Xp Vi) (%pa» Vea) 
woraus sich die behauptete Invarianz sofort ergibt, was zu beweisen war. 


Es sei noch erwahnt, dass in diesem Fall die Matrix A im Koordinatensystem 
der x; die folgende Gestalt erhalt: 


(12’ 


Ferner hat man als Rekursionsformeln fiir die Berechnung von (A) an Stelle 
von (13): 
Sera Pes = (A — Ox) bu — Yu Pas - (13’) 


Die numerische Rechnung wird nun durch Aufrundungsfehler gest6rt; diese 
haben zur Folge, dass die Biorthogonalitatsrelationen (5) und (8) nicht genau 
erfiillt sind. Da die Biorthogonalitat fiir die Eigenwertbestimmung wesentlich 
ist, ist es empfehlenswert, wenigstens die Relation (5) durch geeignete Korrek- 
turen zwangsweise zu erfiillen. Man bildet zu diesem Zweck die innern Produkte 
des auf Grund von (3’) ermittelten Vektors »;,,, mit allen frithern y; und kor- 
rigiert damit %;,,, 4): 


_ (Fer V5) 
V+, I (45, 4';) : 
x und entsprechend fiir die y. (15) 
Nps (korr) — “k+4 — DS Yea; x; | 
j= 


Die Aufrundungsfehler haben aber auch einen Vorteil: Sie verhindern das 
vorzeitige exakte Nullwerden der Vektoren ~;, yx selbst in den Fallen, wo dies 
theoretisch zu erwarten ware, so dass man den Algorithmus in praktischen 


Fallen in der Regel bis zum Ende durchfithren kann’). 


1) Beim n-Schritt-Verfahren weiss man wesentlich mehr iiber die Wirkung der Aufrundungs- 
fehler (vgl. [3], S. 43ff.). Es wird dort auch eine grundsatzlich andere Methode zur Korrektur der 


Orthogonalitat angegeben. ; f 
2) Ein numerisches Beispiel mit einer symmetrischen Matrix A mit mehrfachen Eigenwerten 


(N = 8, m = 7), wo also der Algorithmus nach sieben Schritten beendigt sein mtisste, findet sich 
bei [7]. 
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§ 7. Die Bedeutung der Kodiagonalelemente /, 


Wir wollen den folgenden theoretischen Untersuchungen wieder den Algo- 
rithmus in seiner urspriinglichen Form (2), (3), (4) zugrunde legen, die Gréssen 
f,, stimmen dann mit den Invarianten A, = yz Sp, von § 6 tiberein. 

Die naheliegende Vermutung, dass die Gréssen f;, fiir die Abweichung der 
Xp, Vz Vor System der Haupt- und Eigenvektoren der Matrizen A bzw. A* 
massgebend seien, wird durch das folgende Beispiel widerlegt : 


Bildet man mit diesen Anfangsvektoren bei kleinem ¢ die Kodiagonalform von 
A, so erhalt man 
| l+e e+2¢3 


1 2+36 4e— 28 &? 
it 3—46 


Die f sind also alle sehr klein; trotzdem sind auch im Grenzfall ¢« +0 die 
Vektoren x;,, y; bei weitem keine Eigenvektoren: 


== (19):0 | 0)5 age e(-1)) 26) 0), aya2 oh eae 
(Le Lh) t)., Vo = Ue i eg Vo.—= (0..]. Oh sath 


Dagegen folgt bei Matrizen mit linearen Elementarteilern, wenn man die 
Anfangsvektoren der Bedingung y = R x unterwirft, aus der Kleinheit aller 
B,, auch die Ubereinstimmung der Xp, Vz mit den Eigenvektoren von A bzw. A*. | 
Ferner kann man aus der Kodiagonalform (12) unmittelbar ablesen, dass die 
%, gegen die Eigenwerte konvergieren, wenn die 6, simultan gegen Null streben, 
genauer: 

Satz 3: Sind in der Kodiagonalform einer Matrix A simtliche f-Werte 
~positiv, so gilt der folgende Einschliessungssatz!). Sind wu, v zwei beliebige 
positive Zahlen mit der Eigenschaft uv > Br-++ By (wobei hier und auch 
spater die Gréssen By und fy iiberall, wo sie auftreten, durch 0 zu ersetzen 
sind), so liegt mindestens ein Eigenwert von A im Intervall 


e—-UsSAca,t+y. 


) Dieser ist allerdings nur eine andere Formulierung bekannter Einschliessungssatze, insbeson- 
dere tritt die Grosse B,,_ + f;, in anderer Gestalt bereits bei WEINSTEIN [10] und WIELAND? (11) auf. 


E oe ae vee ne * is ae F 
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Beweis: Man kann in diesem Fall die Kodiagonalmatrix (12) durch Trans- 
formation mit einer Diagonalmatrix reell symmetrisch machen: 


Hy V By | | 


Fir reell symmetrische Matrizen gilt aber auf Grund des Entwicklungssatzes: 
Wenn fiir einen Vektor z die Form T(A, u) = (Sz—Az, Sz — w2)3) negativ 
ist, dann liegt mindestens ein Eigenwert von S zwischen / und y. Nun ist mit 
Be — é,") offenbar 


S é, = V Bus C4 Ay Cx + VB. Cha» 


alsomitA= a, —u, w=a,+0: 
A, ys (Vins Cpa FU Cy + V Br Cht4» V Bus oe ee 2 VB. Css) 
= 6,—uv+ Bus; 


was zu beweisen war. 

Auf Grund dieser Tatsachen ist es wohl klar, dass es fiir viele Anwendungen 
erwiinscht ist, dafiir zu sorgen, dass die f, klein werden. Wir geben deshalb in 
§ 8 einige allgemeine Regeln fiir die Abhangigkeit der 6, von den Anfangsvek- 
wOren! 4,4) all. 

Allerdings wird mit dem Kleinhalten der f, ein Nachteil erkauft: Aus der 
Formel (7) ergibt sich naémlich, dass dann ¥;,,, im allgemeinen viel kleiner als 
x, sein muss. Das ist aber nach Formel (3) nur méglich, wenn x; ,, als Differenz 
fast gleicher Gréssen entsteht, was hinsichtlich der Aufrundungsfehler immer 
sehr ungiinstig ist. Insbesondere leiden die Relationen (5) und (8) darunter, so 
dass man unbedingt die Korrekturformeln (15) anwenden sollte’). 


§ 8. Uber die Gréssenordnung der Kodiagonalelemente 2; 


Um die Abhingigkeit der Parameter «, 6 von den Anfangsvektoren zu un- 
‘tersuchen, betrachten wir die durch die Beziehung (11) definierte einparame- 


1) T(A, “) hangt aufs engste mit dem Templeschen Quotienten zusainmen, es ist namlich das hier 
verwendete 7(A,,¢) bis auf einen Faktor (S z)2 identisch mit dem Ausdruck [7‘(?) —A,) (H2—?#) 
bei [8], S. 174. 

2) Die e,, sind die sogenannten Einheitsvektoren; ¢, = (1 [O|...| 0); eg = (0/1|0| oe | 0) usw. 

3) Praktische Versuche des Verfassers zeigen, dass es in den Fallen, wo theoretisch kleine 
B-Werte (zirka 0,1% der Hauptdiagonalwerte) zu erwarten waren, schon fiir N = 6 fast unmoglich 
ist, den Algorithmus ohne diese Korrekturen erfolgreich zu Ende zu fuhren. 

Fir das Auflésen linearer Gleichungssysteme mit dem 1-Schritt-Verfahren ist dieser Nachteil 


allerdings unerheblich. 
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trige Schar von Anfangsvektoren x(t), y(t), wobei wir zusatzlich annehmen, — 
dass die Vektoren %», Vg und damit auch x(¢) und y(¢) in allgemeiner Lage seien. 
Gemiass (7) und (10) benétigen wir dann fiir die Abschatzung der f nur die 

Gramschen Determinanten G;,(t). 


a) Reelle und verschiedene Eigenwerte 


1, > Ay >... > Ay; die zugehGérigen Eigenvektoren von A bzw. A* seien— 
il 2 ‘N > g' 
Oi Ve; «ss, Un DZWs UT, esos, Uh eS sot aes 


N . 
(x, y) = f(t) = (e4* x, ¥) = Dee, | 
k=1 
wobei : 
4 — (v%, %) (Uz, Vo) ( 
Ki (Uz, vi) . 


bei allgemeiner Lage der Vektoren %g, y» sicher positiv ist; es ist namlich bei 
geeigneter Normierung: vf = Rv,. Auf Grund bekannter Sitze der Determi- 
nantentheorie wird dann 


Gylt) = DP Citar 4a)! V2A. a, 


BY a, 


ae (16) 


fet 


wobei tiber alle Kombinationen von k unter den N Indexwerten 1, 2, vee UN 
summieren ist, und V(A,,...,4,) die Vandermondesche Determinante der 
Werte 1, ..., A, bedeutet. 
Damit sind einmal alle G,(¢) fiir alle ¢ positiv, und man kann den Algorith-_ 
mus fiir jedes ¢ bis zum Ende durchfiihren. 
Das dominante Glied in (16), welches das Verhalten fiir ¢ co bestimmt, | 


erhalt man fiir 7, = 1, 7, = 2, ..., J, = k; es wird somit asymptotisch | 
Cr (Efe BSS a! VN Mee te ee | 
eee ; (17) | 
, x, perk o, erkt T 
(x ke» Vx) = 46 ae VA, hy oS) = 2, @ ¥ W?(A;, A, assy Ax-1) , | 
release WA, basi asa dace) 
ge FB gem ae_t YAR Ay ka 18 
Bis Fiwy el | EG SE GER) (18) 


mit der Abkiirzung W(A| a, dg, ..., ay) fiir (A— a) (A— Ay) ... (A — a,). Fiir 
¢ -> co konvergieren also in der Tat alle B, gegen Null. 

Satz 4: Dann und nur dann, wenn alle Eigenwerte reell und alle Elementar- 
teiler linear sind, kann man die Anfangsvektoren «, y!) so wahlen, dass samt- 


liche B-Werte positiv werden, und zwar gentigt es, x in allgemeiner Lage und 
y = R x zu wiahlen. 


1) Falls das Minimalpolynom nicht den Grad N 


hat, muss man natiirlich in jedem invarianten 
Unterraum geeignete Anfangsvektoren wihlen. 
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Beweis: Die Bedingung ist hinreichend, wie aus den obigen Ausfithrungen 
hervorgeht. Sie ist aber auch notwendig, denn wenn alle f;, positiv sind, kann 
man ja die Kodiagonalmatrix durch Transformation mit einer Diagonalmatrix 
auf die Gestalt (12”) bringen; was zu beweisen war. 


b) Es treten mehrfache, aber nur reelle Ergenwerte auf 


Auf Grund der Voraussetzung, dass das Minimalpolynom den Grad N 
haben soll, bedeutet dies auch nichtlineare Elementarteiler. 
Wir nehmen die Eigenwerte wieder der Grosse nach geordnet an: 


ee eee ae 


und bezeichnen die Vielfachheit von A, mit fielesuist.also. 2/4, == N)..s ist 
dann 


p 
I= DF eo Pkt) ; (19) 
k=1 
wobei die P,,(¢) Polynome vom Grad ju; — 1 sind, deren héchste Potenz bei all- 
gemeiner Lage der Vektoren xy, yy nicht den Koeffizienten 0 haben kann 
(vgl. § 4). Es ist nun G,, der k-te Hauptminor der Matrix 


PRD NLD UD wt TDNESs 

‘ f a 

(Ge. ESE ISS Pee aon ae ree (D-= Ge] 
TSN WIE Pareto MEDI tad 


und setzt sich aus Summanden der Form e(*i*4i.*****4ix)'+ Polynom (é) zu- 
sammen, wobei die 2 im Exponenten mit ihrer Vielfachheit auftreten kénnen. 
Wir unterteilen diese Matrix in #2 Untermatrizen U,,, (g, 4 = 1, 2,..., p) 
mit uw, Zeilen und mu, Spalten und machen dann die folgende Umformung}), 
gegeniiber der siimtliche Hauptminoren und |G| selbst invariant sind: 


My Ms a) 
a OF Oi PU a Oh 
Ms || Up, Use 
uk eh Se ease I 
Mg OF eee St wire! OFF, 
| | 
Wea CH Ca) eee Seas re Wr 


1) Diese Umformung ist eine elementare Determinantenumformung : Addition. des ee ae 
einer Zeile oder Spalte zu einer andern Zeile oder Spalte. Damit jedoch alle Hauptminoren invariant 
i i j , rejter rechts lie 
bleiben, wird immer nur zu einer Zeile oder Spalte addiert, die weiter unten bzw. weiter rechts liegt. 
? 
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Das in der g-ten Zeile und y-ten Spalte der Untermatrix U,, stehende 
Element 
Dit lsteet tga tlt Mat Mateo Fen tt “{f} 


wird ersetzt durch 
(D —A,)" (D —A,)*®...(D—Ajgy (BD —Ayt e— BA cd Dae 
x D—Ayght-" (Deas oe 


: 

; . 

(wobei alle Elemente der Matrix G durchlaufen werden). Auf Grund der Eigen-— 
schaft . 
emt? POG) far is < uy | 

D— hy) fem! P, ee (20) 


0 fir Ss 2 fy, 


werden dann in der Untermatrix U,,, alle Glieder mit e'', ..., e*"—*' eliminiert, 
wobei 1 die gréssere der Zahlen g und / bedeutet. | 

Wir wollen nun zum Grenzfall ¢ co itbergehen und deshalb nur immer 
die Glieder héchster Ordnung in Betracht ziehen, wobei aber nachtraglich 
gezeigt werden muss, dass nicht gerade die mit den dominanten Gliedern ge- 
bildete Determinante verschwindet. In der Untermatrix U, , dominiert offenbar 
in jedem Element ein Glied der Form e*"'+ Polynom, wobei » = Max (g, h) ist, 
also sieht die Matrix G schematisch wie folgt aus: 


ent eo?! ae ent 
| : | 
[fie oa hs 
| e t en! etn! 
The +e 
} etn! : oN! e*Nt | 


Beim Versuch, die Determinante zu entwickeln, erkennt man sofort, dass sogar 
nur die Elemente in den diagonalen Unterdeterminanten U,, entscheidend 


sind, die wegen (20) wie folgt aussehen, wenn man sich auf dominante Glieder 
beschrankt : 


Ha 28 Oy tia se eect cca neal Deca) 
AO, | Re 0 | 


| 
! | 
1 
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Die Polynome 


OMe (D = 1,)"* (D — 1,)?-.. (D — A,_,)?#9-1 {e%2! Pt} 


haben denselben Grad wie P, und denselben héchsten Koeffizienten Coens 
Fs ist nun so, dass man sich fiir ¢ + co wieder auf die héchsten Potenzen von 
t beschranken und deshalb die Polynome Q, durch c, t”“9~! ersetzen kann. 
Damit werden die G,(¢) asymptotisch gleich den Hauptminoren der Matrix 


Lees ae 
Ole: Um OF 8.0 
6) Det Pgemeae te 


re (age jeer (apa ty l 
aad 
Agt I} By ; 
1 (uy — 1)! | 
| (4-1)! 0 Te otha 


wobei U,, =; 


Diese Hauptminoren sind alle von 0 verschieden, was die Weglassung aller 
nichtdominanten Glieder rechtfertigt, es ist namlich (nattirlich immer nur 
asymptotisch)?) : 


Gy lb) aie, ea 
Galt) = — of oP (yy — 1) BE, 


allgemein : 


Gt) = Had) (fir k= 1,2; ..., ay) , wobel 

2 (21) 

H,,(t) = 1) op 1-21-3! se (R= L(g; = 1) abe 

(14; — fe te 1)1 phlua-*) eh Ait 

A VEE 

fur Se Lr PY j=1,...,p) > 
My 

6,6) = (1? ops em (4g — I} = | Ta 


i reic aber elementaren Rechnung 
1) Die nachfolgenden Resultate beruhen auf einer umfangreichen, aber elementaren Kechnung 
nit Determinanten, die hier nicht wiedergegeben werden kann. 


ZAMP Iv/4 
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Fir die Hauptminoren G,, ,; (k = 1, ..., Mg) tritt |U,,| als konstanter Faktor 
auf, also 
G 


Tuy +hk 


(ty s= | Oh) He g(t) Tat eee ee eer 
allgemein : (21a) 
Gacwent) = | Uy, | | Use | 22 | Uj-1, 9-1 Hae 

(fir 2 ch, 20% cali eee ee 


Durch Quotientenbildung gemiss (10) und (7) findet man hieraus schliess- 
lich die asymptotischen Werte fiir die f: 


Br=—k(u,—4)t? (far k=1)2,...,4,—1), 


_ ls __ 4\fr= 1 2 e—At get nk 
Pun vy Cy ( ) : (44 — Lye , 
Bisa RP (ty =k) t* (fr aed (22) 
+s C3 ie) ba—1 A Ag—Aslt fiat Ma 2 
Bran is oP) ee : (fo = ae 
usw. 


Hieraus erkennt man, dass auch in diesem Fall fiir £ + co alle 6 gegen Null kon- 
vergieren miissen, wenn auch im allgemeinen langsamer, als wenn nur ein- 
fache Nullstellen auftreten. Es lassen sich sogar an Hand der Konvergenz die 
Kastchen der Jordanschen Normalform erkennen, denn iiberall dort, wo in der 
Jordanschen Normalform in der Kodiagonalen eine Eins steht, hat man in der 
Lanczosschen Kodiagonalform einen B-Wert, der nur wie ¢-2 gegen Null strebt. 

Die Falle a) und b) zusammenfassend, kann man nun den folgenden Satz 
aufstellen : 

Satz 5: Dann und nur dann, wenn alle Eigenwerte reell sind, kann man die 
Anfangsvektoren x, y so wihlen, dass der Algorithmus beliebig kleine b-Werte 
ergibt. 

Beweis: Es bleibt offenbar nur noch zu zeigen, dass die Bedingung auch 
notwendig ist. Wir betrachten zu diesem Zweck die in §5 erwahnte Ketten- 
bruchgleichung fiir die Eigenwerte, die sich aus der Kodiagonalform der Glei- 
chung (A — AE) & = 0 ergibt: 


0O=A—a, — Py 


; fate af fe e : as far : 
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Fir einen komplexen Eigenwert 2 + ?w ist nun immer | 2 — a, | > uw. Fiir 


| Br 


< w?/4 folgt somit aus der Induktionsvoraussetzung 


Th == Cb, ae Br Ue 
| P Bas lew 


| A— Oni — 


die fiir k = N erfiillt ‘st, sofort 


| Bra fe 
pee cs, 
Ai SSR = | 
und damit | A— «,_, — - Baa F = + 
| Pige= ples k 
| 


Somit gilt sie fiir jedes k, insbesondere k = 1; damit kann aber die oben- 
stehende Kettenbruchgleichung nicht erfiillt sein; was zu beweisen war. 


c) Es treten einfache komplex-konjugierte Eigenwerte auf 


Wir wollen auch fiir diesen Fall noch quantitative Beziehungen fiir die f,, her- 
Beiten: Es sei 1g 7 = Agog =A +t, wobei-noch A,> A>... Age A> Ags go>... 
Es ist dann notwendigerweise z¢ 1 = %« 2. Es wird sich nachtraglich zeigen, 
dass die Betrachtung nicht gestért wird, wenn noch weitere Paare komplexer 
Wurzeln auftreten. 

Da die Eigenwerte verschieden sind, gilt (16) mit dem einzigen Unterschied, 
dass fiir ¢ >oo und k = K + 1 zwei komplex-konjugierte dominante Gleder 


auftreten; es wird: 


Gry 


Cs wots yee Mlemer: 

Si OES GRAY ES Ze) ee gm 

ferry? Cure aeme ers eer NAL ele Agee) etn 
Ge eff S cos (ut) + T sin (uw #) } 


a eee WAT ap AQ 


ee R21)! Vr). Ae, Ad yh 4M) oy in ae oe 


I 


PA Geter t( S24 7%), 
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Die folgenden G werden wieder normal, wobei einfach 4 + 7 und A 7 als 
zwei verschiedene Eigenwerte auftreten, die dominanten Glieder sind aber 
reell; darauf beruht die obige Bemerkung, dass noch weitere komplexe Wurzeln 
die Betrachtung nicht st6ren. Die 6-Werte ergeben sich nach (7) und (10) zu 


p= 2 e(4-4K)# S) cos (4 t) + T sin (ut) , 
a ig W* (hep Ae 
2 S74. 72 
Ben =—"" Ts costa Tanta] (23) 
rere’ [S cos (ut) + T sin (nu 2)] a0 
Bras => eK 3 A)t —p? (S?+ T?) J Spcg WV ged 5 bok on eel 


Die weiteren B-Werte werden wieder normal [Formel (18)]. Somit: 

Beim Auftreten eines Paars einfacher komplexer Wurzeln kann man mit 
einer Ausnahme alle 6-Werte durch geeignete Wahl der Anfangsvektoren be- 
liebig klein machen. Dagegen ist Bx ,, fiir hinreichend grosses ¢ immer negativ 
und zwar < — yu? [wie Formel (23) zeigt]. 

Die Kodiagonalform sieht dann fiir grosses ¢ an der betreffenden Stelle wie 
folgt aus (die e bedeuten irgendwelche kleine Werte): 


1 An 3 


Aus den vier eingerahmten Werten kann man die komplex-konjugierten Eigen- 
werte A + 7 leicht berechnen, wenn die benachbarten f-Werte wirklich klein 
sind, es ist namlich 


1 1 
Aw 2 (Cisr t+ OK 42), /- 7% (e411 — Kx 42)? — Beas - 


Beim Auftreten mehrerer Paare komplexer Eigenwerte treten mehrere 
B-Werte auf, die im Grenzfall nicht verschwinden. Eine Schwierigkeit tritt 
lediglich dann auf, wenn mehr als zwei komplexe Wurzeln denselben Realteil 
haben, dann werden sie durch den Ansatz x — e4#/2 %9 nicht «separiert». 
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§ 9. Die Formeln von Frenet fiir die Vektoren x ,, y, 


. Die durch den Ansatz (11) definierten Anfangsvektoren x(t) und y(é), die 
wir auch dem §9 zugrunde legen, geniigen den Differentialgleichungen 
x" =A x/2, y’ = A*y/2. Durch elementare Umformungen mit Hilfe der Be- 
ziehungen (2) bis (8) erhalt man daraus entsprechende Differentialgleichungen 
fir die tibrigen Vektoren x, ..., xy und y,..., yy des Biorthogonalsystems 
sowie fiir die Skalare «;,, By: 


Vf me . if: il , 
Xe = oe A X_ Bra GSinesten, | eRe oe A* Vig — ope Merete (k= 2a) s apealG) (24) 


Ore = Bx ay Bua ) Bi, — Br (p44 = oy) : (k ae is OL Secs INE 1) (25) 


Ferner: 
(Xz, Vu)" = Xp ie; Vr) ) | 
(A %4, Vx)’ = [az =i Bx =] Brss| (Xp, Var) » (26) 
(A Kay Vu)’ = Ons Ba (Xm Vu) - | 


Die Gleichungen (24) nennen wir die Frenetschen Formeln der Vektor- 
systeme x; bzw. y,, weil sie bei symmetrischer Matrix A, und wenn man die 
-Vektoren so normiert, dass die inneren Produkte (x,, y,) = 1 werden, in die 
Frenetschen Formeln des begleitenden Dreibeins der Kurve E = 4H? E im RN 
ubergehen: 


Oa =e By Xe, Xs ar —fB, DG Bo X3 ; xs a — Bs X, + Bs X, usw. 


Stabilitdt der Gleichungen (25) 


Nachdem durch einmalige Ausfiithrung des Algorithmus Anfangswerte 
ox, Bria (k = 1,..., N) bestimmt wurden, kénnen die Gleichungen (25) offenbar 
ohne Kenntnis der Vektoren x,, y, integriert werden. Hieriiber gilt: 

Satz 6: Sind fiir ein gewisses ¢ die Werte ,, f., -.., By —1 alle positiv, so kann 
man das System (25) bis ¢ = + ce oder ¢ = — oo integrieren, wobei in beiden 
Fallen die 6 gegen Null, die «, gegen die Eigenwerte konvergieren, und zwar 
werden die Eigenwerte fiir > + co in absteigender (% > a > ...), firt > co 
in aufsteigender Reihenfolge geliefert. 

Beweis: Nach Satz 4 haben wir unter den gemachten Voraussetzungen eine 
Matrix A mit nur reellen und verschiedenen Eigenwerten; es gilt also (16), 
wonach die G, fiir kein ¢ verschwinden, die fi, also fiir alle ¢ endlich sind. Nach- 
dem wir die Eigenwerte der Grésse nach numeriert haben, folgt aus (18), dass 
die f fiir t > 0e tatsachlich gegen Null konvergieren, und zwar sogar expo- 
nentiell. Es miissen deshalb nach (25), erste Formel, auch die «, konvergieren, 
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und zwar, wie man durch logarithmische Ableitung von (17) und durch Ver- 
gleich mit der ersten Formel in (26) erkennt, gegen /,. Den Fall t > — oo 
behandelt man durch die Substitution A + — A, A, > —/,; was zu beweisen war. 

Satz 7: Falls alle Eigenwerte der Matrix A reell sind, konvergieren die Vek- 
toren x,(t) bzw. v,(¢) fiir oo gegen die Eigenvektoren der Matrix A bzw. A*. 

Allerdings bezieht sich diese Konvergenzaussage nur auf die Richtung und 
nicht auf die Lange der Vektoren, und dies ist auch im Beweis so zu verstehen. 

Diese Behauptung folgt trivialerweise aus bekannten Satzen iiber das Po- 
tenzieren einer Matrix, da ja die x,(¢) und y,(¢) auch durch Biorthogonalisation 
der Vektoren e4"” x,(0) und e4*‘? y,(0) entstehen; es folgt namlich aus den 
Formeln (3) sofort, dass 


x,(t) = e4 4 x,.(0) + Linearkombination von H(t), «+, p(t) 


und entsprechend fiir die y. 


§ 10. Fast zusammenfallende Eigenwerte!) 


Falls mehrere Eigenwerte der Matrix A nahe beieinanderliegen, ist es 
schwierig, die Anfangsvektoren x, 1, fiir den Algorithmus so zu wahlen, dass 
die Kodiagonalelemente f;, alle klein werden und damit die %, praktisch die 
Eigenwerte sind. Man kann aber durch ahnliche Uberlegungen wie bei mehr- 
fachen Eigenwerten (§ 8, b) zeigen, dass — wenn die Eigenwerte gruppenweise 
nahe zusammenliegen — erreicht werden kann, dass die Kodiagonalform in fol- 
gender Weise in Kastchen zerfallt: 


Xp Be 
1 pen Bows 
i 
Se Pea 
LaDy By 
1 %a+4 


Jeder Gruppe von ¢— nahezu gleichen Eigenwerten entspricht ein solches 
Kastchen, wobei mit wachsendem { (vgl. den Ansatz 11) die zwischen diesen 
Kastchen stehenden Kodiagonalelemente By, Bq usw. viel rascher klein werden 
als die ibrigen £,. 

Es scheint mir nun besonders wichtig, dass man auch den Einschliessungs- 
satz (Satz 3) auf diesen Fall verallgemeinern kann, es gilt namlich: 


1) Nachtrag; bei der Korrektur eingegangen. 


= 2 pj =7) 5 - 4 c 1 C1 ) 
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Satz 8: Sind alle Kodiagonalelemente f, positiv und ist A ein Eigenwert 
der q-p-reihigen Untermatrix 


} Opig Parr 0) 1 

1 An +9 
A,= | I I, 

|| Ag oe | 

: 1 he 


_ so liegt mindestens ein Eigenwert der Matrix A zwischen A+ u und A — », so- 
lange uv > B,+ B, ist. 

Beweis: Ich kniipfe an den Beweis von Satz 3 an: Unter den gemachten 
Voraussetzungen ist die Matrix A zu einer reell symmetrischen Matrix S aqui- 
valent [vgl. (12”)], in der als Untermatrix eine zu A, aquivalente Untermatrix S, 


axbhiuiahee FY = i 
p41 Viboas | 
| | 
| Vieee An+9 | 
Vee 
V Bot Xa | 
q 
= bs Ck e+) 
pl 


ein Eigenvektor der Untermatrix Sj, also S, z= Az. Auf Grund der bei Satz 3 
gemachten Ausfiihrungen gentigt es nun, wenn das innere Produkt 


im 


Sei nun 


(Sz—-Azg—uz, S2—-A2+ v2) 


negativ wird, dann muss ein Eigenwert zwischen A+ wv und A — v legen. 
Nun ist aber 
S z= Sy2+ Cpit VB, Cy + Cq V By Cq+1 » 

also 


q 
Sz-Ae—a2— ra) bo coer caV Ba ast — aD Cy ees 
p+ 1 
und damit endlich, da ja die e, orthogonal sind: 


q 
(Sz—Az—uz, S2- Aztv2z) = chy ppt Bg—Uv Dh. 
; poi 


Dies ist aber in der Tat negativ, sobald wv > B, + Ba; was zu beweisen war. 


1) Vgl. Fussnote 2 auf Seite 45. 
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\ 
Summary . 


The paper begins with a review of the essential points of Lanczos’s orthogo- 
nalization procedure, which is of great importance for the determination of the 
eigenvalues of a real, but otherwise general matrix. Then several properties useful 
for numerical computation are proved: 

If the degree of the reduced characteristic polynomial of the matrix is m, then 
it is possible to choose trial vectors x, y such that the iteration may be continued | 
for exactly m steps. At that point the process must stop because the iterated 
vectors 441, Vm41 Vanish (theorem 1). 

If all the eigenvalues are real, then the co-diagonal elements 2, |see equation 
(12)] can be made arbitrarily small by proper choice of the trial vectors 7, y 
(theorem 5), which considerably simplifies the evaluation of the eigenvalues and 
eigenvectors. Should, in addition, all the B, still be positive (see theorem 4), then 
bounds for the eigenvalues may readily be given (theorem 3). Further remarks 
are made concerning matrices with complex eigenvalues. 

Vinally it is shown that by starting from a certain one-parameter family of 


trial vectors (11), the diagonal and co-diagonal elements a(t), B(t) are solutions 
of a system of differential equations (25). 


(Eingegangen: 8. September 1952.) 
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Theorie der Kombinationsseismographen 


Von Max WEBER, Ziirich!) 


Einleitung und Uberblick 


Die Theorie der Seismographen, wie sie E. WIECHERT [26]2) erstmals ein- 
wandfrei dargelegt hat (die vorangehenden Arbeiten sind in [5] zusammen- 
gefasst), untersucht, wie man mit Hilfe eines mechanischen Systems, das aus 
masselosen Federn, starren Hebeln, Gelenken, StoBstangen und Massen zu- 
sammengesetzt ist, kleine Bewegungen einer starren Masse messen kann. Diese 
Aufgabe hat WIEcHERT geldst. Jedoch haben sich, wie die Erfahrung gezeigt 
hat, seine Voraussetzungen als zu einfach und seine Lésung als zu kompliziert 
erwiesen [1], [18]. Trotzdem ist, wie die Literatur zeigt, das von WIECHERT an- 
gegebene und auf mechanische Seismographen zugeschnittene Messverfahren 
allgemein tibernommen worden. Man hat also, kurz gesagt, die mechanischen 
Seismographen durch elektronische ersetzt und damit wohl beachtliche, aber 
keine grundsatzlichen Fortschritte in der Erschiitterungsmesstechnik erzielt. 

Dieser Sachverhalt hat sich fiir die Entwicklung der Seismometrie, die 
heute auf zahlreiche Arbeitshypothesen angewiesen ist; sowie fiir die Lésung 
vieler Fragen der Technik als d4usserst hemmend erwiesen. Gerne verweist man 
in diesem Zusammenhang auf die grossen praktischen Erfolge der angewandten 
Seismik, vergisst aber dabei, dass die Refraktions- wie die Reflexionsseismik, 
entsprechend ihrem einfachen Arbeitsschema an die Seismographen zum Teil 
andere und vor allem geringere Anforderungen stellt. In klarer Erkenntnis 
- dieser heute nur wenig veranderten Sachlage schrieb WIECHERT in der erwahn- 
ten Arbeit: «Unsere Aufgabe fiir die folgenden Ausfiihrungen ist aber klar vor- 
gezeichnet: Wir werden nicht nur untersuchen miissen, wie man empfindliche 
Seismometer bauen kann, sondern auch, wie lesbare Diagramme zu erhalten 
und zu entziffern sind.» 

Versteht man unter lesbaren Diagrammen, die zu entziffern sind, Diagram- 
me, aus deren Inhalt man nicht nur auf die geologischen Strukturen, sondern 
allgemein auch auf die physikalischen Eigenschaften der Erde schliessen kann, 
so verlangt die Seismometrie in erster Linie eine exakte Erschiitterungsmessung. 

Mit diesen kurzen Hinweisen ist meines Erachtens zur Geniige dargetan, 
dass die Entwicklung einer einwandfreien und praktisch brauchbaren Erschiit- 


1) Institut fiir Geophysik, ETH. we 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis aut Seite 80. 
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terungsmesstechnik heute in vermehrtem Masse einem dringenden Bediirfnis 
entspricht. In der vorliegenden Arbeit wird nun eine verfeinerte Theorie der 
Seismometer gegeben und damit in geeigneter Verbindung mit der Elektronik 
eine neue Methode zur Messung von Erschiitterungen (mechanische Schwin- 
gungen, Normblatt DIN 1311 vom Mai 1939) entwickelt. 


§ 1. Die Aufgabe der Erschiitterungsmessung 


Die Aufgabe der Erschiitterungsmessung ist nach WIECHERT und anderen 
Autoren ((6], [7], [26]) die Messung von kleinen Bewegungen starrer Massen. 
Diese Beschrankung auf starre Massen ist weder notwendig noch vereinfacht 
sie die Aufgabe. Die Bewegungen eines elastischen K6rpers relativ zu einem 
Koordinatensystem sind vollstindig bekannt, wenn man in jedem Zeitpunkt 
die Lage jedes einzelnen Massenpunktes des elastischen K6rpers relativ zum 
Koordinatensystem kennt. 

Verwendet man das rechtwinklige Koordinatensystem S, mit den Einheits- 
vektoren 11, 1, ,, 1, als raumfestes Bezugssystem und betrachtet einen belie- 
bigen Massenpunkt M, so kann man seine Lage durch den Vektor 


t{ %1, Va, a}+ s{6é, , c} 


darstellen, wobei der Vektorr seine Ruhelage und der Vektor 5 seine Bewegung 
relativ zur Ruhelage erfassen mége. 

Fir Erschiitterungsmessungen auf seismischem Prinzip ist das Absolut- 
system S, als Bezugssystem gegeben. Hat aber die zeitliche Anderung der 
Schwerebeschleunigung g auf die Messung keinen merklichen Einfluss, so darf 
man auch ein Koordinatensystem im Erdmittelpunkt oder an der Erdober- 
flache als Bezugssystem wihlen. 

Die allgemeine Aufgabe der Erschiitterungsmessung besteht nun darin, ent- 
weder die drei Komponenten &, n,$ Oder OE/Ot, Onldt, OC/dt oder 07E/0t?, 
0? /0t?, 0?¢/0t? als Funktionen des Ortes r und der Zeit ¢ zu messen. In der 
Technik sind oft einfachere Angaben ausreichend. 

Fiir die Erschiitterungsmessung ist die Tatsache von Bedeutung, dass fiir 
viele praktische Beispiele die einzelnen Komponenten &, 4, € periodische Funk- 
tionen von ¢ sind oder sich wenigstens durch einige allerdings nicht harmonisch 
lhegende harmonische Schwingungen hinreichend genau erfassen lassen [2]. 


§ 2. Der dussere Aufbau 
der elektronischen Kombinationsseismographen 


Ein elektronischer Seismograph ist aus einem Seismometer und einem Ver- 
grosserungssystem zusammengesetzt. Ein elektronischer Kombinationsseismo- 
graph besteht aus mehreren Seismometern in Reihen- oder Differenzschaltung 
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oder in einer gemischten Schaltung und verschiedenen, voneinander unab- 
hangigen, elektronischen Vergrésserungssystemen. 

Ein Seismometer besteht aus einem mechanischen System, kurz auch als 
Schwinger bezeichnet, der so in einem starren Rahmen oder Gehause befestigt 
ist, dass er kleine Bewegungen um eine stabile Gleichgewichtslage ausfiihren 
kann, und einem elektromechanischen Umformer, der die Bewegungen des 


Schwingers relativ zum Gehiuse in elektrische Spannungsschwankungen J 


ubertragt. 


7 
Zn 


Rahmen/ 


Pigs 1 
Schematische Darstellung des Schwingers, 
auf den sich mit wenigen Ausnahmen alle 
praktisch brauchbaren Schwinger zuriick- 
fithren lassen. 


$s 


ute 


Fig. 2 


Schematische Darstellung der wichtigsten 
elektromechanischen Umformer mit Speise- 
spannung. Zn praxt werden im allgemeinen 
nur die aktiven Elemente (eingerahmt) in 


die Seismometer eingebaut. 


Den einfachsten Zusammenhang zwischen V, und der Messgrésse ergeben 
unter der Vielzahl der méglichen Schwinger diejenigen, deren Massenelemente 
sich, wenigstens naherungsweise, auf geradlinigen Bahnen bewegen und damit 
praktisch nur eine Bewegungsrichtung besitzen. Sie lassen sich mit Ausnahme 
der astasierten Pendel formal auf einen Stab, der einseitig eingespannt und 
am anderen Ende federnd befestigt und mit einer Zusatzmasse M versehen ist, 
zurtickfiihren (Figur 1). 

Die elektromechanischen Umformer, die eine Dehnung oder Verschiebung 
oder deren Anderung pro Zeiteinheit messen, lassen sich folgendermassen ein- 
teilen: 

a) Elektromechanische Umformer mit Speisespannung. In Figur 2 sind die 
Umformer mit Speisespannung schematisch dargestellt. In Figur 2a sind die 
beiden Kapazitaten C, und C,, die auch durch «Strain gages» oder Induktivi- 
taten mit beweglichem Kern ersetzt werden kénnen, die aktiven Elemente. In 
Figur 20 ist der Widerstandsstreifen S, und in Figur 2c die bewegliche Spule 
S, das aktive Element. In praxi werden im allgemeinen nur die aktiven Ele- 
mente (in Figur 2 eingerahmt) in die Seismometer eingebaut. 

b) Elektromechanische Umformer ohne Speisespannung. Die Umformer ohne 
Speisespannung sind entweder auf dem piezoelektrischen oder magnetostrik- 
tiven Effekt oder nach dem Tauchspulen- oder Reluktanzprinzip autgebaut. 
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Ein elektronisches Vergrésserungssystem besteht aus einem Verstarker und 
einem Registriergerat (siehe § 6). 


§ 3. Die Beschleunigung, mit der ein beliebiges Massenelement des 
Schwingers belegt wird 


Die Bewegung eines Seismometers, dessen Rahmen sich zusammen mit seiner 
unmittelbaren Umgebung wie eine starre Masse verhalt, lasst sich wie folgt 
beschreiben (Figur 3). Das System S, sei das Absolutsystem. Das Koordinaten- 
system S, sei mit S, starr verbunden und sein Ursprung O, liege in der Ruhe- 
lage des Messpunktes M,,. Das System S; gehe durch eine Parallelverschiebung 
aus S, (Lagekoordinaten €, 7, €) und das System S, durch eine Drehung aus 
S; bei gemeinsamem Ursprung hervor (Lagekoordinaten y, @, y). Ist das 
System S, starr mit dem Seismometer verbunden, so beschreiben &, , € seine 
Translations- und y, , 7 seine Rotationsbewegungen. 

Die Einfiihrung von y, @, x ist nur deshalb notwendig, weil es keine Seismo- 
meter gibt [9], die nur auf Translationsbewegungen empfindlich sind (ent- 
gegen einer vielfach irrtiimlich vertretenen Ansicht, zum Beispiel in [25]). 

Die Beschleunigung b eines beliebigen Massenelementes des Schwingers 
kann im vorliegenden Fall so berechnet werden, wie wenn das Element in seiner 
Ruhelage relativ zum Gehiuse, bezeichnet durch den Punkt P mit den Koor- 
dinaten x4, V4, 2,im Koordinatensystem S,, festgehalten wiirde. 

Die Berechnung von b{S,] = b{bye, Day, by} (mit b/S,,] sei die Darstellung 
der Beschleunigung b im Koordinatensystem S,, bezeichnet) ist in [8] durchge- 
fiihrt. Fiir Erschiitterungsmessungen (§ 2, [24]) ist folgende Naherung ausrei- 
chend, aber, wie gezeigt werden soll, keine Bedingung fiir die Lésung der Auf- 
gabe. Mit 0y/0t — y, ... erhalt man fiir b 


ye =—M (P+ Pp) +Uu(HX+GY+27y) : 
+a(yot+pt+2zyyp)+it+nxy—lo-eo, 
Dyy —%X4 ( : - = YP ) E54 (y* 7 y) 


. : Ete OS eee gare (1) 
PAE 0 te ARAB ae ray SiC ey ae 2 


bap = —% (ypt+ Q) + Vy(— XP+¥) 


-a(e?+p)+égp—Hytl. 


Im allgemeinen wird das Seismometer im System S, eine beliebige Lage inne- 
haben und somit das System S, nicht mehr das dem Schwinger am besten an- 
gepasste Koordinatensystem darstellen. Daher wird ein fiinftes Koordinaten- 
system S, eingefiihrt, das mit dem System S, starr verbunden ist. Beschreibt man 
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die Drehung des Systems S; in bezug auf das System S, mit Hilfe der Matrix 

COS &; cos Bi cos yj 

Y= {| COSck Coss Cosi 

COS%; COSPs, COS Vs, 
{cos «;, cos 3; ocys,; (7 = 1, 2, 3) sind die Richtungskosinus des Systems S, in 
bezug auf das System S,], so ergibt sich fiir die Beschleunigung b[S;] [23] 

BLS] = Y b[S4] . (2) 
Dabei hat man in b[S,] die Koordinaten %,, y4, %, des Punktes P mit Hilfe der 


2, 


Fig. 3 Fig. 4 
Darstellung der eingefiihrten Koordinatensysteme Die Eulerschen Koordinaten. 
Sj, Sg, S3, Sq und Lagekoordinaten EST Gx PN Pn ve 
M,, = Messpunkt. 


Gleichung (d = Vektor 0,0;: O,, O; = Ursprung von Sy, Ss; q, = Vektor OW, 
2 = Vektor O;P) 
Of %4, Va Za} = DP Mel %5, Vs, 5} + Df dae May az} (3) 
durch die entsprechenden Koordinaten *5, V5, 4, Zu ersetzen. 

Die in Figur 3 verwendeten Lagekoordinaten sind so gewahlt, dass kleine 


Bewegungen auch kleine Betrage von y, 7, % ergeben. Es ist aber zu erwarten, 
dass in konkreten Aufgaben andere Definitionen dieser Lagekoordinaten die 


Messung vereinfachen k6nnen. 
Verwendet man zum Beispiel fiir die Beschreibung der Rotationen des Systems 


S, relativ zu S; die bekannten Eulerschen Koordinaten y,;, 9, , (Figur 4, [11)), 
so ergibt sich fiir die Beschleunigung b 
Dy, = —%4 (Pi + 1)? — Va (Gi + iy) + 24 (Vr Oy) a = E+ (Yi t+ FP) » | 
= 44 (G+ Hi) — Veli + py)? + OF] — 24 (Ys BH) — Elvi +o) + i, 
Bi, = 7426, 4 49, — % 8 — HA+ 6. 

Ferner kénnen die Erschiitterungen, Zum Beispiel in Fahrzeugen, einer belie- 
bigen Bewegung iiberlagert sein. Die Berechnung der Beschleunigung b ist fir 
diesen Fall in [23] durchgefthrt. 


(4) 
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§ 4. Der Schwinger 


Figur 1 zeigt den Schwinger, der aus einem elastischen Stab der Lange J, 
der Biegesteifigkeit EJ, der Dichte 9 und des Querschnittes g, einer punktfér- 
migen Zusatzmasse M und einer Zusatzfeder mit der Federkonstanten lo 
besteht (e = Dampfungsfaktor, H — EI/o q/14). Der Rahmen des Schwingers sei 
so mit dem Koordinatensystem S, verbunden, dass die Stabachse mit der , 
x4-Achse und die Schwingungsebene mit der (x, y,)-Ebene zusammenfallt. 

Somit unterliegt der Ausschlag des Schwingers w, (w, = w 1, x, = rl) in der 
yy-Richtung folgender Differentialgleichung 


0?w _ Ow tw by 


pat Ocge + A ee a ee (5) 


die noch durch die Randbedingungen 


(0) OO sad ==, (0). 
or 
; 07w Ow 03w 
oA eee OH +e, aaa L,w— H aa == HG, (6) 
_ HMB : MI?b,  d2w 
fey) i ogl’ aaa eS a 
zu erganzen ist. 
Zur Lésung verwendet man den Ansatz 
~ 
wir, t) > Da C,,(t) wt, (7) (7) 
nel 
mit 
a*u,, A é 5 ae 
dye — My =0 (as = FF, (8) 
der, in (5), (6) eingefiihrt, mit Hilfe der Gleichungen 
du, (1) 
iv sa An Ly u,(1) ee u,(1) 
und 
oO 
a?C dC : a : 
oo Cant + eg" + Cn) ttn = 3 Fatty =F, (9) 
1m . n=1 
lo, <) 
Ly 35h ual) ec 
nel 
ergibt. 


Damit sind die Eigenfunktionen wu, Vollstandig bestimmt. Zu berechnen 
sind noch die Funktionen F, bzw. C,,. 
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Versteht man unter JS den Operator (Stieltjessches Integral) 


S flr) tu dr = fie) ty dr + Ly f(1) g(t) (10) 


so erhalt man fiir die «belastete Orthogonalitat » (14] und fiir die Normierungs- 
bedingung zwischen den Eigenfunktionen 


wud = 0 >| ceeliae cS ue GY ek 


und fiir /,,, unter Beriicksichtigung der Gleichung (10), durch Anwendung von 
S auf Gleichung (9) 


al 


Wis =| Fe, t) Un dy = AH Gy u,,(1) ¥ 


Damit ist die Integration grundsatzlich durchgeftihrt. C,, ergibt sich aus 
der Differentialgleichung 


aPC, dC, é ae 
; qe Ogg Pn Cw Tat 
TF} Setzt man zum Beispiel folgende 
n9 Storungen (allgemeinere Funktionen 
| durch Zusammensetzung nach Fov- 
08 RIER) 
a F(r, t) = rail ger a G, = 7, ei vt 
06;—— (¢ =) —1) 
oe an, so ergibt sich fiir die stationare 
O4 Loésung: 
iwet bi Pi 
- a oe 2 - te 1 é o 
02 mit 
1 
gl = ® 
WG =| A(7) u,,(7) dr ae on u,,(1) = 


= 0 


0 

FT NGL LG G8 TT = fb 49 
Fig. 5 Die Frequenzgleichung lautet 

Graphische Darstellung der Frequenzgleichung (11). (Figur 5) : 


iby 7 lg 1 + cosa cosh « a (11) 
~ Hat te ee ee Sane cos % sinh« 
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Mit den Abkiirzungen 
C,(r) = cosa,” + cosha,7, S,(r) =sing,7 + sinha, 7, 
c,(v) = cosa,” —cosha,7, S,(7) =sina, 7 — sinha, 7 


erhalt man fiir die normierten Eigenfunktionen 


ae , 
tal?) = By (= 2G) enlt) + 54(0)) (12) 
SHGQIAE (om Gab) ype (aly sinh %, COS %, + Cosh «, Sin &, 
Sy sad 
— ay (EY + 2.L1 cal) sal) 


sinh «, COS %, — COSh %» SiN & — C,(1) sp(1) 


Sol) + Ly s2(1). 


Paes 


Spezialfalle 
a) M=0, 7, = 0 (freies Ende). 
Frequenzgleichung: 
cosa cosha + 1= 0 
(%, = 1,875, a, = 4,694, a3 = 7,855, a= 10,996, [10]). 


Eigenfunktionen: 


ane Sy (7) ? 


b) fo > co (gehaltenes Ende). 
Frequenzgleichung: 
j tga = tgha 


(a = 3,927, a = 7,069, a3 = 10,210, a, = 13/352 [10]). 


Eigenfunktionen : 
1 
U,(r) = —c,(r) + tgha, Spl?) - 
c) M> eq! (Federpendel) . 
Mit sina ~ « — «9/3, cosa > 1— «2/2 usw. geht Gleichung (11) und Glei- 
chung (12) iiber in 


2 Sf Leta ; 
yp? i: + fo , Ulr) as (5 vie i 2) 
M VL, D 5 


I 


uw stellt bis auf einen konstanten Faktor die elastische Linié dar. 


§ 5. Die Eigenschaften der Seismometer 


Ein Seismometer besteht, wie in § 2 erwahnt, aus einem Schwinger (§ 4) 
und einem elektromechanischen Umformer. 
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Die verschiedenen Umformer lassen sich einheitlich erfassen. Ist J", ein 
Operator, das heisst eine Vorschrift, die jeder Funktion /(v, t) eine Funktion 
Ln f(7, #) mit ry = 79, das heisst der Variablen ¢ allein zuordnet, so ist (e,, = Emp- 
findlichkeit des Umformers) fiir kapazitive oder induktive Umformer in Briik- 
kenschaltung (m = 1) 


Ty LIN t) ee ey (f(r, t)),. To 


fiir «Strain gages» oder naherungsweise piezoelektrische Umformer (m = 2) 


Taf.) = (5, Ht, 9) 


r=1q 


und fiir induktive Umformer ohne Speisespannung (m = 3) 


Py flr,t) = es (7 fv.) 


P= To 


Damit ergibt sich fiir die allgemeine Seismometergleichung (1, = Eimheits- 
vektor in Arbeitsrichtung): 


t 
Veeeelee tire. 1 yn / F(t) v,(¢ — t) dt u,(7) , (13) 
n= 1 aS R 


1 
‘ R= (/ ste PY 4, (0) dr + (tng) wt) 


by iS ae ane 
Pni(é) ae = ice oe a | ake 


Ve — 42 


Eine Diskussion der Greenschen Funktionen ¢,,(¢) findet man in [15]. Als 
Arbeitsrichtung bezeichnet man die Bewegungsrichtung des Schwingers. , 
Zur Berechnung des Amplituden- und Phasenganges nimmt man zum Bei- 
spiel an, der Schwinger sei mit dem Koordinatensystem S, starr verbunden 
(Arbeitsrichtung parallel 44) und die Beschleunigung se1 
iot 


ban = Dy: =0, Day = 1 = at aw" é 


& 


Dies in Gleichung (13) eingesetzt, ergibt 


yf ait 


ae is @? No Fm 610! Uni?) | un(r) dy + Ly wa) a) 
. 0 


y2 — w% + Lew 


ZAMP Iv/5 
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Bezeichnet man den Effektivwert E,,(w@) der Ausgangsspannung V, fiir 
Ho = 1 als Wegempfindlichkeit, so folgt aus Gleichung (14) fiir E,,(@) 


1 
1 — pe 2 iwt tn 7 
Bolo) = | Se peta tee ( f mle drt teat] 


—wr+tiew 
V2 | w= 0 


und damit fiir V, 
Ve = V2 to Eyl) efet-w). 


, ist durch die Gleichung 


lo 6) . 
27 tw t f° 
/9 ™~ i(wt—-Ow) wo Pm e1o U», (7) 1 d ' il 
g - Se ty(v) dr + Ly u 
V2 E,,(o) e M Sr meer / n(7) 1 @,(1) 
n=l 0 
bestimmt. 


Definiert man entsprechend E,(@) und E,(@), so folgt 


Vi =V2 19 o p(w) eV, = 2 49 w® Ey(c) ef (t-P 
mit 


p= o, = d, . 


Die Empfindlichkeiten sind gegenseitig durch folgende Gleichungen ver- 
bunden: : 


E(w) = w E,(@) = w? E,(a) . 


Kir ein vorgegebenes Seismometer, charakterisiert durch die Konstanten 
é, Ly, Ls, v4, ..+, %, den Operator J",, und die Eigenfunktionen w,, beschreibt 
somit die Wegempfindlichkeit als Funktion von w (einzelne Zahlenangaben in — 
der Literatur beziehen sich im allgemeinen auf den linearen Arbeitsbereich oder 
auf die statische Empfindlichkeit E,,(0), E,(0)], den Amplituden- und die Funk- 
tion #,, den Phasengang des Seismometers als Wegmesser (in den Diagrammen 
erscheint #, im ZeitmaBstab #,,/a). Analoges gilt fiir die Funktionen E,(@), . 
Ey(@), 0,(@), 0y(@), die den Amplituden- und Phasengang des Seismometers 
als Geschwindigkeits- bzw. Beschleunigungsmesser beschreiben. 

Werden jetzt die Funktionen 


E(w) e'?, E(w) e'®?, E,(@) eb 


in der GauBschen Ebene mit @ als Parameter dargestellt, so ergibt dies ein 
gutes Bild tiber das Verhalten eines Seismometers gegeniiber Erschiitterungen 
(Figur 6). 

Die Empfindlichkeiten, zum Beispiel E,,, sind derart von w abhangig, dass 
jedes einzelne Glied in der entsprechenden Reihendarstellung proportional einer 
Funktion von w ist, die fiir alle Glieder vom gleichen Typus ist. Die méglichen 
Funktionstypen sind in Tabelle 1 zusammengestellt (2, = 1/72 — w2+7¢ @). 
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Fig. 6 
Der Amplituden- und Phasengang eines Seismometers als Weg-, Geschwindigkeits- und 
Beschleunigungsmesser (AZ = fy = 0, W = /94, €/vy= 1/2). 


£,(0,L)/£p (0,50) 


Fig. 7 
Die relative Empfindlichkeit £,(0, L,)/Ey(9; oo) fur punktférmig bei y = 1 angreifende 
Umformer( 14, 1’). 
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Ferner ist die Empfindlichkeit proportional e,, und von L, sowie L, und 
bei Verwendung von J’, auch von 7,, der aktiven Linge des Umformers, ab- 
hangig. Es sind dies keine einfachen Funktionen. Im folgenden sind zwei durch- 
gerechnete Beispiele graphisch dargestellt. 

1. Fir punktformig bei y= 1 angreifende Umformer (J), I’). 

In Figur 7 ist die relative Beschleunigungsempfindlichkeit E,(0, L,)/E,(0, co) 
in Funktion von L, (L3 = 0; w = 0; 1, ..., ¥%, = const) dargestellt. 


Tabelle 1 
Die méglhichen Funktionstypen fiir die verschiedenen Empfindlichkeiten und Umformer 


Ties lke 
2 
wo? Q . 
E i Sgt iw Q 
2 e const fiir w >», nm 
E w 2, (€ + 0) tw? Q, 
a ew const fiir wm 9, w const fiirw >», 
E, Q, ; tw Q, (€ + 0) 
~~ const fiir w < », ~ const fir w & », 


eal Lr Umformer, die durch J’, erfasst werden. 
Mit v = VE/o erhalt man fiir F,(0): 


ies. 
2V3 wo A, G, 
E,(0) = es V : nn 
(0) = 2, ue 


1 
K, = / u,(r) dr +L, u,(1), G 
0 

In Figur 8 ist die relative Empfindlichkeit Ey (0, rq, L,)/Ey (0, 74, 0) in 
Funktion von r, fiir verschiedene Werte von L, (L3 = 0) als Parameter einge- 
tragen. 

Zur Definition des neutralen Punktes geht man von der Annahme aus, der 
Schwinger sei mit dem Koordinatensystem S; mit der x;-Achse als Stabachse 
(¥;=0 bis x,=2) und der (x, y¥5)-Ebene als Schwingungsebene starr ver- 
bunden. 

Das System S;, soll aus dem System S, durch eine Parallelverschiebung in 
Richtung der —x,-Achse um den Betrag dy (dg = dy) hervorgehen. Somit ist (3) 


X4 = x5 — dy 5 y= V5 : 24 = 25 p 


- ee F : : 
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_ Die Beschleunigung sei 


Osa = O52 = 0, d5y = (%5 — do) X = (dy — %5) Yo wre”. 


vo 


- Damit erhalt man fiir V, (J d = dy) 


ae 1 
yy 4 0 di m ae “n ; 
ite wx eiwt u, (yr) (| (d—-r\ u,(r) dy + (d—1) 1, “st : 
0 


2 { , 
at a Oe D 


£4,(0; r9; Ly) / Eg(0, ra; 0) l,=G966 
1=Y 


‘a ] 


1, =O,664 
14 


4 =,219 


L,=G,lé2 


Fig. 8 
Relative Empfindlichkeit E,(0, 7g, L,)/E,(0; %, 0) in Funktion der aktiven Lange 7, des 
Umformers fiir verschiedene Parameterwerte L, (L; = 0). 


Hat die Gleichung 
V, = 0 (fiir alle t-Werte) , (15) 


die eine Bestimmungsgleichung fiir d ist, eine reelle Losung, wie gezeigt werden 
soll, so bedeutet dies, dass auf der Stabachse insofern ein ausgezeichneter 
Punkt (neutraler Punkt) existiert, als durch die Verlegung des Ursprungs des 
Systems S, in diesen Punkt (Abstand von der Einspannstelle des Stabes = dp) 
die Beschleunigungen infolge der Drehbewegungen des Systems S, eliminiert 
werden. 

Die Bedeutung dieses Punktes ergibt sich aus folgender einfacher Uber- 
legung: Der Schwinger sei mit dem Koordinatensystem S, mit der ,-Achse als 
Stabachse (x, = J, bis x, = J, + /) und der (x4, ¥4)-Ebene als Schwingungsebene 
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starr verbunden, und die Beschleunigung sei },, = (/, + lv) ¥ (0 Sr < 1) oder, 
umgeformt, 


by =lh+nx¥=lh+d—dt+ny, (i=l). 


Laut Definition des neutralen Punktes (15) ergibt die Beschleunigung 
| (—d+-+ 7) x keinen Beitrag zur Ausgangsspannung V,, oder die von 7 abhan- 


gige Beschleunigung / (/, + v) 7 kann durch diejenige des neutralen Punktes, 


namlich / (1, + d) 7, ersetzt werden. 


Setzt man analog 
1 


[r u,(v) dy + L, u,(1) 
0 
a To 
/ u,(7) dy ate I, u,(1) 
0 


, 


1 
ee 1 Ti eiat 
E,, a V2 vw —owrtiew (| un(r) dr + Ly v(t) 


0 / 


(d,, = Abstand des neutralen Punktes und E,,, = Beschleunigungsempfindlich- 
keit der n-ten Partialschwingung), so folgt 


oa V2 w® x Ad, E,, 8 (mtn) = 0, 
n=1 


Ad, = a= d, 


n? 


das heisst den Beitrag jeder Partialschwingung kann man als Beitrag eines 
Federpendels (+ Umformer) mit der Beschleunigungsempfindlichkeit E, und 
n 


oe 
der Ausgangsspannung V,, (V, = >» V,), das im Abstand Ad, vom neutralen | 
n= 1 


Punkt des Seismometers eingesetzt wird, auffassen (Schwingungsrichtung 
parallel Arbeitsrichtung). Somit stellt die Gesamtheit dieser Federpendel 
(+ Umformer) in Reihenschaltung ein einfaches Ersatzschema des Seismo- 
meters dar. 

Wie aus Gleichung (15) ersichtlich, ist die Lage des neutralen Punktes eine 
Funktion von w. In Figur 9 ist der Abstand d (L, = 0, L; = 0, I, e = 0) in 
Funktion von /», aufgetragen, 

In Figur 10 ist ebenfalls der Abstand d fiir den Fall L,=0,L;=0,I,,e=0 
in Funktion von w/y, fiir verschiedene Parameterwerte 7, dargestellt. Fiir die 


optimale Lange 7, (das heisst d nahezu konstant) ergibt sich 0,44 (lineare Inter- 
polation), 
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0 O5 hp a ae) ZF) a7 


tenis) 
Der Abstand d des neutralen Punktes von der 


Einspannstelle des Stabes in Funktion von m/v, 
(L, = 0,L, = 0,Iy,€ = 0). 


d 

740 
730 

r=G4 
720 

r=Gd 
710 

r=O2 
700 
r=G/ 

690 x 7 = 
600 \ = 


a G2 G4 G6 48 


Fig. 10 
Der Abstand d des neutralen Punktes fur 
dens Pall) 2-0), —0, 56 — Oh und 
ftir verschiedene Parameterwerte 7, 
(optimaler Wert von 7, = 0,44). 


Fig. 1 


Vier nach der vorstehenden Theorie gebaute und geeichte Scismometer. 


wy, 
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§ 6. Die elektronischen Vergrésserungssysteme 


Die Verstarker haben zwei Aufgaben. Sie dienen einerseits zur Anpassung 
der Seismometer an das Registriergerat und liefern anderseits bei Verwendung 
von geeigneten Schaltungen, als Rechengerate, einen Beitrag zur Auswertung 
der aufgenommenen Erschiitterungen. Ihr Aufbau ist somit im wesentlichen 
vorgezeichnet. Sie bestehen je aus einem Vor- und Endverstiarker und einem 
Zwischenverstarker als Rechengerat. 


Fig. 12 
Schaltschema eines RC-Bandfilters. Die Werte der Schaltelemente sind in UF bzw. MQ angegeben. 


Zur Lésung der vorliegenden Aufgabe (§ 1) bendtigt man als Rechenelemente 
differenzierende und integrierende Vierpole [20] und Bandfilter. Ihre funktio- 
nelle Bedeutung ist in den §§ 7 und 8 dargelegt. 

Die Bandfilter werden entweder zur Begrenzung des Frequenzbandes und 
damit zur Verbesserung des Stérpegels oder zur spektralen Zerlegung von 
Signalen (zum Beispiel Fourier-Analyse) verwendet. Fiir die Erschiitterungs- 
messtechnik sind, trotz ihrer relativ kleinen Kreisgiite, Bandfilter auf der Grund- 
lage von riickgekoppelten RC-Filtern sehr geeignet. Als Vorteile sind ihr relatiy 
kleines Gewicht, ihr geringer Platzbedarf sowie ihr einfacher Aufbau und ihre 
kleine Stéranfalligkeit anzufiihren. 

Eine allgemeine Theorie der RC-Bandfilter, die etwa mit der Theorie von 
W. CAUER [3] fiir verlustfreie LC-Filter vergleichbar ware, gibt es noch nicht. 
Sie diirfte auch wesentlich schwieriger sein. 

Das in Figur 12 wiedergegebene Schaltschema ist daher als ein Beispiel zu 
betrachten und wurde «halb theoretisch, halb experimentell» gefunden. Die 
Figur 13 zeigt die gemessene Frequenzkennlinie. Die erste und zweite Stufe bilden 
die untere (Hochpass) und die dritte und vierte Stufe die obere Flanke (Tiefpass). 

Auf eine eingehende Beschreibung der Vergrésserungssysteme kann in An- 
betracht der ausgedehnten Fachliteratur [17], [20] verzichtet werden. 
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Fig. 13 


Die gemessene Frequenzkennlinie des in Figur 12 dargestellten RC-Bandfilters 
(Durchlass fiir 100 Hz = 100%). 


ah aw 


Fig. 14 
Der elektronische Teil eines nach der vorstehenden Theorie gebauten 
6-Kanal-Vergrésserungssystems. 
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§ 7. Das Schema des exakten Seismographen 


Die allgemeine Seismometergleichung (13) kann, gestiitzt auf das im §5 
hergeleitete Ersatzschema, mit 


ae SE (10) 


n=1 


als Losung des Gleichungssystems 
Amin Vn = Om On. GAS, 2,3 eee ee (17) 


aufgefasst werden. b,, ist die Beschleunigung des neutralen Punktes des -ten 
Federpendels. 
Die linearen Operatoren haben zum Beispiel fiir m = 1 folgende Form: 


2 


d Cae 
A,, = dt? + & ae Tr Vy 


Fir das lineare Vergrésserungssystem gilt, wenn der Lichtzeigerausschlag 
mit @ bezeichnet wird, die Gleichung 


Pa—@PV, (18) 
mit 
U Inj ; 
D; = Aig + ay + +++ + Ayn, Se iy = 1, 2,...) 


Soll nun auf dieser Grundlage ein exakter Seismograph, das heisst ein Seis- 
mograph, der nach Erreichung des stationaéren Zustandes des Vergrésserungs- 
systems alle zur Aufzeichnung gelangenden Erschiitterungen amplituden- und 
phasenrein vergréssert (man vergleiche [6]), aufgebaut werden, so ist folgende 
Aufgabe zu lésen: Bei gegebenem Seismometersignal V, sind die Operatoren 
, und ®, so zu bestimmen und elektronisch darzustellen, dass zwischen a und 
b, — sie sind gegenseitig, wie aus dem Ersatzschema hervorgeht, durch einfache 
Operationen verbunden — folgende Gleichung zu Recht besteht (uw = Vergrésse- 
rungsfaktor des Vergrésserungssystems) 


a= —Mey b,. 


Mit Riicksicht auf die relativ kleine Bandbreite der Erschiitterungen ist 
eine Naherungslosung, deren Giiltigkeit auf die Umgebung einer Eigenfrequenz 
des Seismometers (vorwiegend Grundfrequenz) beschrankt ist, véllig ausrei- 
chend. Damit ergibt sich naherungsweise aus den Gleichungen (16), (17) und 
(18) fiir die Differentialgleichung des Seismographen 


D, nin V, > —QD, em by 
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und aus dieser Gleichung, wie man leicht erkennt, fiir die Operatoren @, und@,: 


D, = inva D, =e vale. n Ve 
oder 


Cie Vee 

Die elektronische Realisierung dieses Vergrésserungssystems ist unter den 
gemachten Voraussetzungen einfach. Es besteht zum Beispiel fiir m = 1 aus 
einem Vorverstarker, drei getrennten Zwischenverstarkern mit den Ubertra- 
gungsoperatoren 

ae pe th gn ae 
A AOR i oy 2 ae und fy “dite? 

einem Bandfilter (§ 6), einem Endverstarker und einem Registriergerat, zum 
Beispiel K. O. In praxi erzielt man die richtige Einstellung der Zwischenver- 
starker, indem man die Verstaérkung der einzelnen Verstarker variiert, bis die 
freien Seismometerschwingungen keinen Lichtzeigerausschlag mehr ergeben. 
Die Bandbreite des Filters ist durch die zulassige lineare Amplituden- und 
Phasenverzerrung [23] festgelegt. 

Fiir ein Seismometer mit einem Federpendel als Schwinger (Spezialfall c, 
§ 4) ergibt das vorgezeichnete Schema eine exakte Lésung. 


§ 8. Die Messung von &, 7, 6, &, 7, im E, 1, & mit Hilfe eines 
elektronischen Kombinationsseismographen 
Fir die Beschleunigung b MURS Dea, by} gilt allgemein [23]: 
Bay = Ry + % Ry + Va Rit % Ais » 
Day = Ry + %4 Ray + Va Roo + 2 Rog , 
byz = Rg t Hq Roy + Va Rog + 2 op 5 


4y 


X4, V4 und 2, sind jeweils die Koordinaten des neutralen Punktes. 

Die Vergrésserungssysteme seien so beschaffen, dass 

1. ein exakter Kombinationsseismograph entsteht, das heisst jede beliebige 
Kombination von Seismometer und Vergrésserungssystem einem exakten Seis- 
mographen entspricht (uw; = “4 = const) und 

2. die Bandfilter der einzelnen Vergrésserungssysteme gegenseitig liickenlos 
abgestuft sind. Die Anzahl der Vergrésserungssysteme ist daher durch die Fein- 
heit der Unterteilung der Bandbreite der Erschiitterungen gegeben. 

Mit diesen Vergrésserungssystemen und einem oder zwei oder » + ™ Seismo- 
metern kénnen nun folgende Gréssen, spektral zerlegt, gemessen werden. 
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a) Mit einem Seismometer. Man verlegt den neutralen Punkt des Seismo- 
meters in den Ursprung des Koordinatensystems S4 (%4 = V4 = 2, = 0) und er- 
halt, wenn die Arbeitsrichtung zum Beispiel parallel zur x,-Achse ist, 


1 
k= —== & 
lad em 


und entsprechend allgemein (Schema [) 


1 
a inl ayo == 12s) 


b) Mit zwet Setsmometern. Man verlegt den neutralen Punkt des einen Seis- 
mometers in den Ursprung des Koordinatensystems S, (%,= Vy = 2, = 0) und 
den des andern zum Beispiel auf die x,-Achse im Abstand x, und erhalt, wenn 
beide Arbeitsrichtungen zum Beispiel parallel zur z,-Achse und die Signalspan- 
nungen in einer Differenzschaltung vereinigt sind, 


Ry = — aE 31 
und entsprechend allgemein (Schema IT) 
Ree ore a ters (om 13 71,23) 
Rivpraits on vo ng) k= 2p 9 = 1, 2,3) 
Ruy = - 5 = aes (R= 33 7 =1,2,3) 


c) Mit n+ m Seismometern. Mit analogen Anordnungen kann man allge- 
mein folgende Grésse messen: 


m/2 


x lies Rus t+ Dy R,. 
ka 


ame 


Mit 3 x 2 Seismometern erhalt man zum Beispiel y? oder p? oder ?. 


Die Messung bzw. Bestimmung von &, , ¢, é, ..., Is mit maximal zwei Seis- 


mometern und zum Beispiel mit Hilfe der Lagekoordinaten y, p und yx. 
Es ist daher [§ 3, (1)] 


Reet tx-6o-¢o, Reetepey. ke oe 


Ry = —7? — gp? usw. 
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a) In erster Naherung. In erster Naherung setzt man y= p=7=0 (reine 
Translationsbewegungen) und misst nach SchemalI R, fiir k = 1, 2,3 und 
erhalt 

1 “ iL ae 1 


eee a ee 
1 2 i) 
em be em had . em fad 


und durch Integration 


i 1 
f= — 7 fama, n=... E=-pof fade, g=... 
” ul z 7 Cm : é / 


Die Integration kann auch elektronisch durchgefiihrt werden (siehe ferner 
Spezialfalle § 5). 

f) In zweiter Naherung. In zweiter Naherung setzt man alle quadratischen 
Glieder, wie zum Beispiel y? oder g y, gleich Null und misst nach Schema I R,, 
fiir k = 1, 2, 3 und nach Schema II R, ; zum Beispiel fiir k = 1, 7 = 3 und k = 3, 
#2 und erhalt: 


a3 1 “ ees 1 
ES re IEF Meare siti em | (19) 
1 ‘pene 1 | 
oe en we “43 om et. 


Aus diesem Gleichungssystem ergibt sich fiir é, 7 und fe: 


ie 1 oO cr ote il & i 
2 er ag i ane (a3, ] | on): 
em Le (4: ee a4 | [4 J em M a V4 ve 


und durch Integration FS aA ee . we 

Die Auflésung des Gleichungssystems (19) kann, sofern zwei Kombinations- 
seismographen gleichzeitig verwendet werden, auch vollstandig elektronisch 
durchgefiihrt werden. : : 

y) In dritter Naherung. In dritter Naherung vernachlassigt man 4 % — ¢ 9, 
fy ++ Cy gegentiber —g y, g y und misst nach Schema I AR; ftir & = 1, 2,3 
und nach Schema II R,, fiir k= 2,7 =1; R=3,7=1; Rae 34 2: hea: 
fee = 2 = 2, k= 3,47 = 3und k= 71, 4 = 2 und erhalt: 


E_go=T,, itegv=h, fv-tyt+t=h, 7-¢y=Tt, 
—yy—P=Ty, —x9+ P= In, PEC sits or ess 
p+ y= Te, y+ Gy+2pp=Tp- 
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Dieses Gleichungssystem ergibt fiir y mit 1 + pw? ~ 1 die Gleichung 


1 a + Ae a Tae |r 7 aa Tt 2. 
Ty yp + Tg, + %% Sie ee V Tz3 + Too rica 


Aus y gewinnt man @ mit Hilfe der Gleichung 


“ 1 
p = (py — Ty) yee Fog 


Peay 


Setzt man yw und @ in die Gleichungen fiir 7,, 7, und 73 ein, so ergibt sich 
a) und ¢ durch Auflésung dieses Gleichungssystems und &, 7, ..., € durch 
Integration. 

Diese Integration diirfte im allgemeinen rechnerisch nicht gelingen und ist 
daher zum Beispiel mit einem Tischintegraphen durchzufiihren. 

0) In vierter Naherung. Die vierte Naherung wird durch R, und R,; darge- 
stellt. Die in der dritten Naherung fiir y und g bestimmten Werte kénnen un- 
verandert tibernommen werden. Die Lagekoordinate y kann, da sie in den R,; 
nirgends auftritt, ohne Integration nicht ermittelt und daher die Bestimmung 


von S 7” und z nicht allgemein durchgefiihrt werden. 
Die vollstandige Lésung, das heisst die Bestimmung von &(r, ¢), ... ergibt 
sich durch systematische Abtastung des Messobjektes. 


§ 9. Fehlerrechnung 


Die zur vorstehenden Theorie fiihrenden Annahmen sind in praxi teilweise 
nur naherungsweise erfiillt oder erfiillbar. Die dadurch bedingten Abweichungen 
sollen als Fehler bezeichnet werden. 

a) Die Einzelmasse am freien Stabende des Schwingers wurde als punkt- 


formig angenommen (§ 4). Bezeichnet man die Trigheitsmomente der Einzel- 
masse mit 0,, @, und O, wobei 


—, = oy (%, ly? dm ¢ —, = ah yi dm, 0 = ~O, 1. 0, 
M 


ist, so erhalt man fiir y= 1 an Stelle der Randbedingungen (7) die Randbe- 
dingungen: 


07?w Ow Fw 03 w 02 w 
Rie Sa | | 1 = 7 Suter Als a 
Ey ot? rT Pac sd xed Time hese ae Or Ot? oye = Gy, 
OHI : l 
Ly = TE G, = TE (— ©, Ro, + Os Ry») . 


Aus diesen Gleichungen folgt: 


1, Das Seismometer ist in der Arbeitsrichtung und in der Richtung der Stab- 
achse (@, PR.) empfindlich. 


2. Ist @, klein, so muss, damit die Einzelmasse relativ zum Stab noch als 
starr angenommen werden kann, auch @, klein sein. 
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3. Ist O, und ©, bzw. O klein, so muss, konsequenterweise, in der Differential- 
gleichung (5) des Schwingers auch der Tragheitswiderstand der Drehbeschleuni- 
gung mit berticksichtigt werden. Dies ergibt asynchrone Stabschwingungen, und 
damit ware eine Auswertung der Messergebnisse praktisch unméglich. 

Ein Seismometer ist daher zur Messung in zweiter oder héherer Naherung 
nur dann geeignet, wenn 


0, ~~ 6, RY 0 
oder 
0°w 
HE. ee HG,~ 0 
aS. 


b) In § 3 wurde angenommen, die Beschleunigung b sei unabhangig von den 
kleinen Auslenkungen des Schwingers aus seiner Ruhelage relativ zum Gehduse. 
Ist zum Beispiel die Arbeitsrichtung eines Seismometers parallel zur z,-Achse 
und sind %,, 0, 0 die Koordinaten des neutralen Punktes, so lautet im Spezialfall 
(§ 4, c) die Differentialgleichung des Schwingers: 


d?w dw 


dP Fé i Pa) as Sirs = Boy TR. (20) 


Beriicksichtigt man die Relativbewegung des Schwingers, so geht die Glei- 
chung (20) tiber in 
d*w dw 


qa t oq + P+ Ry) w= — Ry — % Ry. (21) 


Allgemein ist vy? durch v2 + R,,; (7 = 1, 2, 3) zu ersetzen. Dies bedeutet, dass 
die Annahme in erster oder zweiter Naherung exakt erfiillt ist. 
Zur Messung in dritter oder vierter Naherung ist notwendig, dass 


Riy<v 


or<r 


oder, da Ry;~ w? ist, 


ist, das heisst, es sind Beschleunigungsmesser erforderlich. 

Der Einfluss des Terms R3, w entspricht einer Modulation und ist nach dem 
Verfahren der aufeinanderfolgenden Verbesserungen in [22] berechnet. 

Es sei noch beigefiigt, dass Schwinger, die durch die Differentialgleichung 
vom Typus (21) beherrscht werden, im allgemeinen instabile Bereiche haben 
k6onnen. 

Formal analoge Verhaltnisse entstehen zum Beispiel auch bei einem Horizon- 
talpendel, dessen Aufhangepunkt vertikale Schwingungen ausfihrt [13], [19] oder 
bei einem Schwinger, der senkrecht zur Arbeitsrichtung beschleunigt wird [16]. 

c) Der Schwinger wurde als linear elastisch angenommen. 

Wird bei der Herleitung der Differentialgleichung (5) die Kriimmung in 
zweiter Naherung beriicksichtigt, so ergibt sich 


0?2w Ow 04w 3 bei Oz [se () | as (I S 


Or ar Ot nae or* 2 or? | Or? Or 


Als Mass fiir nichtlineare Verzerrungen wurde in [23] der Klirrfaktor J’ ein 
gefiihrt. 
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In erster Naherung ist: 


po hia (Ga ea 
seid; 6 Vi = (3 m)2)2 + €? (3@)? 


1 1 
: i a 3 Ae 0? | 0?w /dw\? . 
y= [Fi dy So &4 u,(1), 1, = SHe-siot [ ee (a) Ju dy . 
0 


0 


Zam Schlusse méchte ich Herrn Prof. Dr. F. GASSMANN fiir das grosse 
Interesse, das er dieser Arbeit entgegenbrachte, sowie fiir viele Anregungen und ~ 
Diskussionen herzlichst danken. Ebenso bin ich den Herren R. BERGER, Fein- 
mechaniker, und M. DIETIKER, Radiomonteur, die mich beim Aufbau des Seis- 
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Summary 


The basic principles of vibration measurements and their applications are 
briefly discussed followed by the development of a new method. The essential 
features of a combination-seismograph and the operation of its several parts are 
outlined. A general seismometer equation is developed and applied to the study 
of various types of practicable seismometers. The theoretical concept of the 
construction of a precise seismograph is then worked out, and a description is 
given of its practical realization in the form of a precise combination-seismograph. 
A discussion of errors of the new method of measuring mechanical vibrations 
finishes the paper. 


(Eingegangen: 17. Juni 1952.) 
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Zum Verteilungsgesetz der Tropfengréssen bei der Zerstaubung 


Von Hans A. TROESCH, Strengelbach, und PETER GRASSMANN, BPH, Zurich 


Bei der Zerstaubung fallt immer ein Gemisch von Tropfen sehr verschiedener 
Gréssen an. Von verschiedenen Seiten wurde darauf hingewiesen, dass die Zahl 
der Tropfen in den einzelnen Gréssenklassen — das Verteilungsgesetz — wahr- 
scheinlich statistischen Ursprungs sei [MEHLING'), LEwis”)]. 


1) H. Menuine, Zur Physik der Brennstoffstrahlen in Dieselmaschinen, ATZ., H. 16, 411-421 


1934) (Auszug aus Dissertation). ae, . 
3 o (Gs Tae D. G. Epwarps, M. J. Goeita, R. I. Rice und L. W. Smiru, Atomization of 


Liquids in High Velocity Gas Streams, Ind. Eng. Chem. 40, Nr. 1, 67-74 (1948). 


ZAMP I1V/6 
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Um aber mit den Methoden der klassischen Statistik, die zum Beispiel zum 
Maxwellschen Geschwindigkeitsverteilungsgesetz fiihren, auch hier zu eindeutigen 
Ergebnissen zu gelangen, bedarf es zweier Erhaltungssatze. Sie entsprechen de- 
nen, die auch dort als Grundlage dienen: Die Satze von der Erhaltung der Masse 
und der Gesamtenergie. Der erste ist auch fiir die Zerstaubung, soweit man von 
Verdunstung absieht, unmittelbar gegeben. Natiirlich gilt fiir die Zerstaubung 
auch der Energieerhaltungssatz. Aber er niitzt zunachst wenig, denn von der 
Energie, welche der zu zerstaubenden Fliissigkeit primar zugeftihrt wird, verwan- 
delt sich bekanntlich der iiberwiegende Teil in Warme. Nur ein ganz geringer, von 
vornherein nicht angebbarer Bruchteil, erscheint schliesslich als Oberflachenenergie 
des Tropfengemisches. 

Zur Umwandlung der dem System primar zugefiihrten Energie in Warme be- 
darf es aber einer gewissen Zeit. Gerade bei hochwertigen Zerstaubertypen muss 
sich das schliesslich erhaltene Tropfengemisch in sehr kurzer Zeit ausbilden: Der 
zu zerstaubenden Fliissigkeit wird eine hohe Relativgeschwindigkeit zum um- 
gebenden Gas erteilt, wobei gréssere Tropfen durch die an ihrer Stirnflache auf- 
tretenden hohen Staudriicke aufgeblaht und zerrissen werden, ein Vorgang, der 
sich beim im Erdfeld frei fallenden Tropfen in Bruchteilen einer Sekunde abspielt. 
Bei den wesentlich hoheren Relativgeschwindigkeiten, wie sie fiir die Zerstaubung 
in Frage kommen, diirfte er aber in einem viel kiirzeren Zeitraum ablaufen, so 
zum Beispiel bei einem Wassertropfen von 70 « Durchmesser in 2,5 >* 10-* bis 
10-%s. Nach dieser Aufteilung muss das Tropfengemisch méglichst rasch in einen 
Raum kleinerer Tropfendichte iibergefiihrt werden; denn durch die Abbremsung 
am umgebenden Gas verringert sich rasch die hohe Relativgeschwindigkeit, so 
dass durch Zusammentreffen mehrerer Trépfchen wieder gréssere Tropfen ent- 
stehen kénnten. Das in einem Raum hoher Relativgeschwindigkeit ausgebildete 
Verteilungsgleichgewicht auf die einzelnen Gréssenklassen muss also so schnell 
wie méglich «eingefroren» werden. In den fiir die Ausbildung der Verteilungs- 
funktion massgebenden sehr kurzen Zeiten wird also im wesentlichen eine Um- 
wandlung von kinetischer Energie in Oberflachenenergie und umgekehrt statt- 
finden, wobei sich ein Gleichgewichtszustand zwischen diesen beiden Energie- 
betragen einspielt. Erst nach Ejinstellung des Gleichgewichts verwandelt sich 
dann die noch vorhandene kinetische Energie tiber Reibungsvorgange in Warme. 
Nach diesem Bild des Zerstaubungsvorgangs wird also ein fiir jeden Zerstaéuber- 
typ charakteristischer Bruchteil der primar dem System zugefiihrten Energie als 
gesamte Oberflachenenergie schliesslich im Tropfengemisch enthalten sein. Wie 
er aufgeteilt wird, mit anderen Worten, wie sich die Tropfen auf die einzelnen 
Grossen verteilen, ist durch die Gesetze der Statistik gegeben. Man kann also ver- 
. suchen, in die statistische Behandlung des Zerstaubungsvorgangs die Konstanz 
der gesamten Oberflachenenergie des Tropfengemisches an Stelle der Konstanz 
der Gesamtenergie des Gasgemisches der statistischen Mechanik einzufiihren. Der 
Vergleich mit den experimentellen Ergebnissen muss dann erweisen, ob dieser 
Ansatz eine giinstige Naherung an die tatsachlichen Verhiltnisse darstellt. 

Im Gegensatz zu den Betrachtungen, die zum Maxwellschen Geschwindig- 
keitsgesetz fiihren, sind aber im Falle der Zerstéubung die méglichen Energie- 
werte nach unten begrenzt. Sehr grosse Tropfen, und diese entsprechen Zustanden 
niedrigey Oberflachenenergie, sind instabil. Der Tropfen wird namlich nur durch 
die Oberflachenspannung zusammengehalten, und der dadurch bedingte Druck 
wird relativ zu den tbrigen am Tropfen angreifenden Kraften um so kleiner, je 
grésser der Tropfen ist. Durch eine derartige Stabilitatsbetrachtung lasst sich — 
wie in einer spateren Arbeit noch im einzelnen dargelegt werden soll — fiir jeden 
Zerstaubungsvorgang eine maximale Tropfengrésse ®,,,, ableiten. 
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Auch nach oben sind die méglichen Energiewerte beschrankt, da sehr kleine 
Tropfen wegen ihres héheren Dampfdrucks instabil sind. Da ausserdem kleine 
Tropfen hohen Werten der Oberflachenenergie entsprechen, ist ihre Zahl im all- 
gemeinen schon durch die Bedingung der Konstanz der Oberflaichenenergie so- 
weit erniedrigt, dass diese zusatzliche Grenze wenigstens fiir die bisher untersuch- 
ten Falle keine wesentliche Bedeutung zu haben scheint. Sie wurde deshalb in der 
folgenden Betrachtung nicht beriicksichtigt. 

Tatsachlich gelangt man, wie von TROESCH gezeigt und spater im einzelnen 
dargestellt werden soll, auf Grund dieser beiden Annahmen zu einem Verteilungs- 
gesetz der Tropfengrésse bei der Zerstaubung, das die experimentellen Ergeb- 
nisse in besserer Weise wiederzugeben vermag als das bisher hiefiir meist beniitzte 
Rosin-Rammlersche Gesetz. 

Die Zahl dN aus der Gesamtzahl N der Molekiile, die sich im Tropfen mit 
einem Durchmesser zwischen ® und ® + d@ befinden, ist gegeben durch 


dN = N q(®) d@, (1) 


wobei y(®) eine zundchst noch unbekannte Funktion des Tropfendurchmessers ist. 
Dabei sind die beiden Nebenbedingungen zu erfiillen: 


fan =N, (2) 


fedN = 1 (3) 
In ihnen ist jeweils von ® = 0 bis = @,,,,, zu integrieren. Es bedeutet E die im 
Tropfengemisch insgesamt enthaltene Oberflachenenergie, wahrend ¢ die auf das 
Molekiil eines Tropfens mit dem Durchmesser @ entfallende Oberflachenenergie 
darstellt. Ist o die Oberflachenspannung und v das Volumen eines Molekiils, so 
wird e = 600/®. Die Funktion g(®) ist dadurch gegeben, dass sie die unter Be- 
achtung der beiden Nebenbedingungen wahrscheinlichste Verteilung der Mole- 
kiile auf die verschiedenen médglichen Tropfengréssen darstellt, also diejenige 
Verteilung, die die Zahl der Komplexionen zu einem Maximum macht. Durch die 
Methoden der klassischen Statistik gelangt man dann schliesslich zum Boltzmann- 
schen Energieverteilungssatz in der Form 


dN = 9(®)d®e°*e”  (dundy Konstanten). (4) 


Wenn die Zahl dN der Molekiile bekannt ist, die Tropfen mit einem Durch- 
messer zwischen ® und ®+ d®@ bilden, so ist die Zahl dn dieser Tropfen gegeben 
durch sa 

v 

ee 5 

oe eae (5) 

Es zeigt sich, dass der beste Anschluss an die Erfahrungswerte gefunden wird, 
wenn g(®) = const gesetzt wird, also mit gleichen a-priori-Wahrscheinlichkeiten 
der verschiedenen Tropfendurchmessey gerechnet wird. Damit folgt dann als Ver- 


teilungsgesetz der insgesamt , Tropfen auf die einzelnen Grodssenklassen 


DE). ie oy ® mit Bp=—6avd. (6) 
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Da bei der Ableitung dieser Beziehung die Konstanz der Oberflachenenergie 
[Gleichung (3)] vorausgesetzt wurde, folgt, dass hierbei — im Gegensatz zur Rosin- 
Rammler-Verteilung — die Oberflache nicht unendlich wird, auch wenn man die 
Integration bis zu beliebig kleinen Tropfen erstreckt. Die Zahl dieser sehr kleinen 
Tropfen wird eben gerade auf Grund der Energiebedingung so klein, dass das 
Integral nicht mehr divergiert. 

Tragt man die nach Gleichung (6) berechnete Volumen-Haufigkeits-Kurve in 
das Rosin-Rammler-Netz ein, so ergibt sich eine Kurve nach Art der Figur 1, die 
auch tatsachlich haufig beobachtet wird’). 


Be ee ieee ee 
Se yee pa wre cote chem Orn PI TIPS SE 
RUMI OG aRMEee el Ses a 
enna NS i Ad 2 
repeal | 
a | 
y : feshah 
wae 
aor | 79% \qoai 503 49% 90% 2p% | eR 
PA 05 10 cera 
R 


Fig. 1 
Volumen-Summen-Kurve nach Gleichung (6) in Rosin-Rammler-Koordinaten fir verschiedene 
Werte von B 


x 


100 Vo I ™ sg 
a Saar aS ae we FS? + Bia) E1(—«) = / args du = Integrallogarithmus, 
pst EE iat Ess ads 
Bl eB + Ei(—B) bee 
a= B wae 
@D 


Punktierte Kurve: Rosin-Rammler-Verteilung. 


Die nach Gleichung (6) im Tropfengemisch vorhandene Oberflachenenergie 
lasst sich fiir gegebene Stoffgréssen und Anfangswerte durch eine Funktion von 
(i= B/Pmaw darstellen. Diese Funktion von # spielt bei der Zerstaubung eine 
ahnliche Rolle wie die absolute Temperatur in der statistischen Gastheorie. Tat- 
sachlich kann auch f, zumindest fiir gewisse Zerstaubertypen, als konstant an- 
gesehen werden. 

Bei einigen Zerstaubertypen kann allerdings nicht mit einem einheitlichen 
De gerechnet werden. Es spielen sich bei diesen verschiedene Zerstaubungs- 
vorgange mit entsprechend verschiedenen @®,,,,-Werten neben- und hintereinan- 


der ab. In diesem Fall muss das Verteilungsgesetz durch Uberlagerung mehrerer 
Verteilungen nach Gleichung (6) gewonnen werden. 


1) E. Girren, J. Proc. Inst. mech. Eng. 158, Nr. 2, 205/206 (1948). 
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Summary 


The distribution problem is, in the first instance, a problem of surface energy 
of the formed drops. It is shown that the distribution of drop size from atomiza- 
tion follows statistical laws. The formulae derived with methods of the classical 
statistics permit calculation of the distribution-curve with the aid of measurable 
parameters for certain types of atomizers, i.e. relative velocity between liquid 
and air, surface tension, viscosity of liquid and air, and air density. 

The formulae for the distribution of drops, surface and volume of the drops 
contain two parameters, the maximum drop diameter and the exchange-factor. 
The former can be calculated by aid of a stability criterion on a single drop, as it 
will be shown in another article. The exchange-factor, however, which is of similar 
signification as the temperature in MAXWELL’s distribution of molecul velocities 
in perfect gases, depends on the type of atomizer. 

The new formulae describe the experimental data much better than all known 
functions (e.g. Rosin-Rammler). However, for their use, which allows the know- 
ledge of distribution curves without experiments, it is necessary to know the 
formula of maximum drop diameter and to make an analysis of the operation of 
the atomizer-type, permitting the establishment of the exchange factor. 


(Eingegangen: 23. August 1952.) : 
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International Union of Pure and Applied Physics 


Under the patronage of the International Union of Pure and Applied Physics 
the following meetings will be held during 1953: 


(1) Meeting of the International Commission of Optics and Symposium on the Sight 
in Connection with Instrumental Optics 
April 1953 at Madrid (Spain), on the occasion of the 50 years of the Sociedad 
Espafiola de Fisica y Chimica. 
Further information can be obtained from Professor J. M. OTERO NAVASCUEs, 
Instituto de Optica ‘‘Daza de Valdés’’, Serrano 121, Madrid. 


(2) International Congress on Electroacoustics and International Symposium on 
the Sound Insulation of Light-Weight (Structures) 

June 1st to 6thincl., at Delft, Eindhoven and Hilversum (Netherlands). The 
Congress and Symposium are organized by the Netherlands Acoustical Society 
under the auspices of the International Commission on Acoustics of the Inter- 
national Union of Physics. 

Sections: Phonography and Recording (Disc, Film, Wire, Tape) ; Public 
Address Systems; Acoustic Measurements; Hearing Aids and Audiometers ; 
Electroacoustics in Ultrasonics; Sound Insulation of Light-Weight (Structures). 

Further information can be obtained from Mr. P. A. pe Lance, Laboratorium 
voor Technische Physica, Mijnbouwplein 11, Delft (Netherlands). 


(3) Commission and Symposium on Cosmic Rays 


July 1953 at Bagnéres-de-Bigorre, France. Detailed information can be ob- 
tained from Professor L. LEPRINCE-RINGUET, Secretary of the Commission, 17, 


rue Descartes, Paris 5°. 
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(4) Symposium on Theoretical and Fundamental Physics 


September 1953 at Kyoto (Japan) organized by the Japanese National 
Committee. Sections: Field Theory and Theory of Elementary Particles; Sta-@ 
tistical Mechanics and Low Temperature Physics; Solid State Physics and Mole- 
cular Physics. Further information from Dr. YOSHIA Fuj1okA, Secretary of the 
Conference, Science Council of Japan, Theno Park, Tokyo, Japan. 

A, PERRIER, Swiss National Committee of Physics 


Instruments and Measurements Conference Stockholm 1952 


Von der K6niglichen Schwedischen Akademie der Ingenieurwissenschaften 
(IVA.) und der schwedischen Gesellschaft fiir technische Physik (TFF.) organi- 
siert, fand vom 22. bis 25. September 1952 in Stockholm die 6. Konferenz tber— 
Instrumente und Messmethodik statt. Gleichzeitig wurde eine Ausstellung von ~ 
Messinstrumenten durchgefiihrt (22. bis 30. September 1952). 

Das Ziel dieser Tagungen ist die Darstellung und Diskussion neuer Messmetho- 
den sowie die Vorfiihrung neuer Instrumente zur Férderung der Forschung und — 
der allgemeinen technischen Entwicklung. Die diesjahrige Tagung wurde von 
rund 400 Teilnehmern besucht. In Parallelsitzungen wurden 80 Vortrage tiber 
Messtechnik auf folgenden Gebieten gehalten: 

a) Mechanik, b) physikalische Chemie, c) Elektronik, d) Kernphysik, e) Spek- 
troskopie, f) allgemeine Elektrik, g) industrielle Kontrolle. 

Die Ausziige dieser Vortrage sind in den IVA. (Tidskrift for Teknisk-Vetens- 
kaplig Forskning) Bd. 23, Heft 4 (1952), erschienen. Wie in den Vorjahren werden 
die vollstandigen Referate in den « Transactions of Instruments and Measurements 
Conference» (Stockholm 1952) ver6ffentlicht werden. 

Sehr interessant und reichhaltig war die Ausstellung von Instrumenten und 
insbesondere die Abteilung fiir Forschung. Hier hatten verschiedene Institute 
der skandinavischen Lander sowie die kernphysikalischen Zentren von England 
und Frankreich in Bildern und Modellen einen Ausschnitt aus ihrer Tatigkeit 
gezeigt. N, SCHAETTI 


3. Osterreichischer Mathematikerkongress vom 9. bis 14. September 1952 . 
in Salzburg 


Die Osterreichische Mathematische Gesellschaft gestaltet seit einiger Zeit ihre 
Jahresversammlung zu einem internationalen Mathematikertreffen aus, das 
auch dieses Jahr unter der Leitung des Vorsitzenden der Gesellschaft, Prof. 
Dr. H. INzINGER, aufs beste gelang, indem etwa 300 Delegierte aus den verschie- 
densten Laindern teilnahmen. 

In der Sektion fiir angewandte Mathematik wurden folgende Vortrage ge- 
halten: D. R. RUTHERFORD (St.Andrews): The Relaxation Method. A. ZILLER 
(Strassburg): Valeur de t, du dernier terme du développement en série de Taylor. 
E, STIEFEL (Ziirich) : Uber die Methode der konjugierten Gradienten zur Auflésung 
von Systemen linearer Gleichungen unter spezielley Beviichsichtigung des Einsatzes 
von Rechenautomaten, S. EKELGF (Goteborg): Der mechanische Diffeventialanaly- 
satoy dey Chalmersschen Technischen Hochschule und einige damit geliste Probleme. 
E, Dory (Léwen): Sur des orientations récentes de l’enseignement supérieur dans 
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le domaine des mathématiques appliquées. H. R. PANETH (London): Correlation 
Functions. H. Ricuter (Freiburg i. Br.): Zur Begriindung des Additions- und 
Multiplikationssatzes der Wahrscheinlichkeitstheorie. D. A. KAPPos (Erlangen) : 
Die Totaladditivitat der Wahrscheinlichkeit. HH. BERGSTROM (Goteborg): Uber 
Folgen von Vertetlungsfunktionen, die eine stabile Grenzverteilungsfunktion haben. 
S. Vajpa (Epsom): Algebraische Theoreme dey Streuungszerlegung. B. TEDESCHI 
(Rom): Leggi di capitalizzazione e leggi di sconto. L. Cortatz (Hannover): Rand- 
wertaufgaben monotoney Art. H. GOERTLER (Freiburg i. Br.): Neue Beitrdge zur 
Berechnung laminarey Grenzschichten. W. GROEBNER (Innsbruck): Schwingungen 
von elastischen Kreisvingen mit vechteckigem und dreieckigem Querschnitt. P, Laa- 
SONEN (Helsinki): Uber Eigenwertaufgaben bei Differentialgleichungssystemen und 
thre iterative Losung. G. GRIOLI (Padua): Sull’equilibre dei corpi elastici. R. FAURE 
(Hanoi): Intégrale premiére en mécanique ondulatoive. M. J. DE SCHWARZ (Rom): 
Zur Saint-Venantschen Theorie hohler Prismen. J. L. DESTOUCHES (Paris): Suy 
deux problémes mathématiques non rvésolus posés récemment pay les théories micro- 
physiques. M. Caz1n (Garches): Applications de la transformation de Laplace a la 
physique théorique. R. NaRDINI (Bologna): Uber zwei Eindeutigheitssdtze fiir ma- 
gnetohydrodynamische Wellen. A. PIGNEDOLI (Bologna): Recenti ricerche sui moti 
di parziale veloci. E. STIEFEL 
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Lecons sur les principes de l’électrodynamique classique. Par A. MER- 
CIER (Editions du Griffon, Neuchatel 1952.) 75 p., fr. s. 7.80. 

In dieser Schrift werden Ausschnitte der klassischen Elektrodynamik einer 
eingehenderen Betrachtung unterzogen. Dabei wird auf einen logischen Aufbau 
der Theorie Gewicht gelegt, was eine gewisse Vertrautheit mit den behandelten 
Problemen bedingt. Andrerseits wird mancher Leser nach der Lektiire erkennen, 
dass noch blasse oder schon verblasste Kenntnisse tiber den Zusammenhang 
etwa von B und E mit H und D zu neuem Leben erwacht sind. 

Nach einleitenden Worten iiber Veranlassung und Plan der Arbeit wird in 
Kapitel 3 die Dynamik, in Kapitel 4 die Energetik des elektromagnetischen 
Feldes behandelt. Die Verkniipfung mit der speziellen Relativitatstheorie bringt 
in Kiirze das Kapitel 5. Das nachste analysiert zunachst eingehend den Begriff 
der elektrischen und magnetischen Polarisation; hierauf kommen deren Eigen- 
schaften im Sinn der Kapitel 4 und 5 zur Sprache. Mit Hilfe der Cliffordschen 
Zahlen wird in Kapitel 7 zusammenfassend eine kurze Axiomatik der Elektro- 
dynamik dargestellt, im Sinne von friiheren Veroffentlichungen des Verfassers. 
Die straffe Gliederung und Beschrankung auf das Wesentliche, im Verein mit der 


bewahrten Drucktechnik des Verlages, sichern dem Werk eine gute Lesbarkeit. 
2 W. Baumgartner 


The Theory of Electromagnetic Waves. A Symposium (Interscience 
Publishers, New York 1951). 393 pp.; $6.50. 

Die in diesem Buch ver6ffentlichten Vortrage — 21 an der Zahl — behandeln 
spezielle Probleme aus der Theorie der elektromagnetischen Wellen. Dabei ist 
aber durchaus nicht immer das physikalische Moment, sondern auch dasjenige 
der rechnerischen Methodik betont; oft ist der Zusammenhang mit dem Titel 
des Buches nur lose und verfiihrt direkt dazu, das Gelesene in seiner Wirkung 
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bei anderen Aufgaben zu erproben. Es ist also seine Lektiire nicht bloss fiir den 
Spezialisten einer Ausbreitung elektrischer Wellen in Hohlleitern und in der 
Atmosphire von Interesse. Selbst der Mathematiker von Zunft diirfte etwelche 
Anregung zu eigenen Uberlegungen finden. Von Vorteil ist der zusammenfassende 
Charakter mancher Artikel, der durch ausfiihrliche Literaturzitate unterstrichen 
wird. Nicht zuletzt ist schliesslich die saubere Ausstattung des Buches bemer- 
kenswert. Einige Inhaltsangaben mégen zur Illustration dienen, wobei ihre Aus 
lese nicht etwa eine Wertung bedeuten soll; sie wiirde je nach Leser verschieden 
ausfallen: , 

H. Levine und J. ScHWINGER referieren iiber Beugung an einer Offnung im 
ebenen leitenden Schirm. R. E. LANGER berichtet iiber asymptotische Lésungen 
der Wellengleichung mit variablem Brechungsindex, H. Portrsky tiber die Ver: 
allgemeinerung der Weylschen Integraldarstellung der Kugelwelle auf Wellen be- 
liebiger Form, H. BREMMER iiber das WKB.-Verfahren, S. A. SCHELKUNOFF 
behandelt die Wellenausbreitung in geschichteten Medien, ahnlich N. MARcu- 
vitz und B. FRIEDMANN. M. Kine bespricht eine asymptotische Lésung de 
Maxwell-Gleichungen, S. O. Ricr die Reflexion an rauhen Oberflachen usw. 
W. Baumgartner 


Allgemeine mathematische Berechnungen auf Brunsviga-Doppel 
rechenmaschinen. Von F.WACHENDORF, bearbeitet und erweitert von B. SCHRA 
DER. Herausgegeben von Brunsviga-Maschinenwerke-AG., Braunschweig, im 
Selbstverlag, 1951. 


Es hat sich hauptsiachlich fiir geodatische Berechnungen als niitzlich erwiesen, 
zwei gekoppelte Handrechenmaschinen zu benutzen, die es erméglichen, Teil- 
resultate von der einen auf die andere Maschine zu transferieren und so eine ge- 
wisse Speicherkapazitat zu erreichen. Die Brunsviga-AG. baut Maschinen dieser 
Art und gibt fortlaufend Anleitungen iiber ihre Benutzungsmdéglichkeit heraus. 
Erwahnt seien: Das Rechnen mit komplexen Zahlen und Vektoren in der Ebene. 
Gleichungen und Schnitte von Kegelschnitten. Lésungen von quadratischen und 
kubischen Gleichungen, Trigonometrische Rechnungen, speziell Vorwarts- und 
Riickwartseinschneiden in der Geodisie. E, Stiefe 


The Preparation of Programs for an Electronic Digital Computer 
By M. V. WiLkes, D. J. WHEELER and S, Gitt (Addison Wesley Press Inc. 
Cambridge [Mass.], 1951), 162 pp.; $5.00. 


Nebst einer Beschreibung der mathematischen Eigenschaften der elektroni 
schen Rechenmaschine EDSAC (University of Cambridge, England) enthalt das 
vorliegende Buch eine ausfiihrliche Anleitung zum Herstellen von Rechenpro 
grammen fiir diese Maschine. Insbesondere wird im Detail behandelt, wie man 
Unterprogramme in einen Rechenplan einfiigt; ferner wird auch gezeigt, wie man 
allfallige Fehler im Rechenprogramm entdeckt und beseitigt. Es folgt sodann eine 
Aufstellung von solchen Unterprogrammen (mit detaillierten Angaben), wie sie 
in Cambridge andauernd verwendet werden. Diesen kann man entnehmen, wie 
man mit der EDSAC e®, log x, \x* und andere Funktionen berechnet, Differential- 
gleichungen integriert usw. 

Da die iiberaus reichen Erfahrungen der Verfasser mit Rechenmaschinen it 
diesem Buch ihren Niederschlag gefunden haben, wird jeder Mathematiker, der 
sich mit der Herstellung von Rechenprogrammen’befassen muss, das Buch mit 
Gewinn studieren, obschon es naturgemass auf die EDSAC zugeschnitten ist. 

HH, Rutishause 


PHYSIKALISCHE BLATTER 


In dieser Monatsschrift finden Sie die zusammenfassende Darstellung der naturwissenschaft- 
lichen Grundlagen und eine standige Unterrichtung Uber den gegenwartigen Stand der physi- 
kalischen Erkenntnisse. 

48 Seiten DIN A 5, Abb., Tab., 9. Jahrgang 1953. Bezugspreis monatlich DM 2.25 


PHYSIKALISCHE VERHANDLUNGEN 


Durch Autorreferate von den Tagungen zahlreicher wissenschaftlicher Gesellschaften vermittelt 
diese Zeitschrift eine aktuelle Kurzinformation lange vor dem Erscheinen der Originalarbeit in 
dem entsprechenden wissenschaftlichen Fachorgan. 
4. Jahrgang 1953, jahrlich 600-700 Kurzreferate, mindestens 300 Seiten. Erscheinungsweise 
spatestens 4-6 Wochen nach den Tagungen. Bezugspreis jahrlich DM 15.-. 


ATOMWAFFEN 


Unter Verwendung des offiziellen Berichts 
der US.Atomic Energy Commission versucht 
diese Veréffentlichung, von wissenschaft- 
licher Seite her das Grundsatzliche der 
Wirkung von Atombomben und die Méglich- 
keit von Schutzmassnahmen darzustellen. 
72 Seiten, DIN A 5, 30 Abb., 8 Tab., bro- 
schiert DM 4.80. 


PHYSIKER-ANEKDOTEN 


Wolfgang GAEDE zu einem Kollegen: ,,Ach, 

ich suche schon seit langem den richtigen 

Assistenten. So einen vorziiglichen Mann, 
ganz ohne eigene Gedanken." 


177 literarische Blitzlichtaufnahmen aus 
dem Leben bekannter Physiker, 56 Seiten 
DIN A 5, viele Randzeichnungen, DM 4. 80. 


PHYSIK VERLAG, MOSBACH, BADEN 


Jetzt wieder lieferbar: 


OSTROWSKI 
Vorlesungen iiber Differential- und Integralrechnung 
Band I: FUNKTIONEN EINER VARIABLEN 
in zweitem, berichtigtem Nachdruck 


Von Prof. Dr. A. Ostrowski, Professor an der Universitat Basel 
373 Seiten mit 42 Figuren. Ganzleinen Fr. 36.40, brosch, Fr. 32.25 


«Das Gesamtwerk wird ein grundlegendes Lehrbuch unserer Epoche sein, das den 
Studierenden in modernem Geist in die klassische Analysis einfiihrt, ohne ihn mit 
den modernen Begriffsbildungen und Methoden zu tberfiittern.» 


(W. Siiss, Archiv der Mathematik.) 


Band II: Differentialrechnung auf dem Gebiete mehrerer Variablen (1951) . 
Band III: Integralrechnung auf dem Gebiete mehrerer Variablen (erscheint 1953). 


Zu beziehen durch Ihre Buchhandlung 


VERLAG BIRKHAUSER : BASEL UND STUTTGART 


in jeder Techn 
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